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ABSTRACT 


A seismo-acoustic sonar concept that uses guided interface waves (Rayleigh or 
Scholte) is being developed to detect buried ordnance in the sea floor and beach 
sediments. This thesis describes the initial research conducted into the design, 
construction, and field testing of possible seismic sources that excite preferentially the 
interface waves desired for use in such a system. The theory of elasticity shows that 
seismic interface waves have elliptical particle velocity orbits in the vertical plane along 
the path of propagation. It was therefore decided that to selectively excite the desired 
interface waves, a harmonic source employed at the interface must induce elliptical 
particle motion in this plane. Several exploratory sources were developed to produce this 
type of excitation. Field tests of the discrete-mode sources developed were conducted to 
evaluate this hypothesis, but due to the non-optimum nature of the experimental sources, 
perfect discrete source excitation was not obtained. However, it was found that the 
medium itself acted as a selective filter for the interface waves after a few tens of meters 
of propagation. The experimental results obtained here suggest that the basic concept of 


discrete-mode excitation looks promising. 
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I. INTRODUCTION 


The research described in this thesis 1s a continuation of work started at the 
Applied Research Laboratories of the University of Texas at Austin, sponsored by the 
Office of Naval Research, into the development of a seismo-acoustic sonar for the 
detection of buried ordnance using guided, seismic interface waves. Their research 
identified certain key issues requiring resolution prior to further concept development. 
The primary focus of the present research 1s to address and resolve some these issues, 
specifically; the design, construction, and field testing of a discrete-mode excitation 
source for seismic interface waves. 

This thesis contains six chapters dedicated to this purpose. The introduction 
describes the necessity and motivation for this research. It also includes a brief 
exposition on current and emerging mine warfare technology as it relates to buried 
ordnance detection. The second chapter provides an introduction into general seismology 
with a detailed description of interface waves and the characteristics that make them 
ideally suited for such an application. The third chapter describes the previous research 
efforts conducted 1n this and similar areas. The following two chapters contain the major 
portion of the information on the discrete-mode source development and the experimental 
results, through the ground and surf zone testing phases. The final chapter provides 
concluding remarks and recommendations for continuing research efforts, which should 


serve to further future project development. 


A. RESEARCH OBJECTIVE 


In November of 1996, during the Mine Warfare Symposium held at the Naval 
Postgraduate School in Monterey, California, the concept of a seismo-acoustic sonar for 
the detection of buried ordnance was introduced to the mine warfare community, and was 
based on the work conducted by the Applied Research Laboratories of the University of 
Texas at Austin (ARL:UT). [Ref. 1] 

The main objective of the ARL:UT study was to demonstrate that an electro- 


mechanical transducer could be employed to excite seismic surface waves, of which there 


are many types, but especially Rayleigh waves, which could be used to ensonify a buried 
mine-like target, and which in turn would produce seismic echoes that could be received, 
processed, and displayed. Due to the generally large contrast in density and elasticity 
between buried ordnance and the surrounding soil, the target produces seismo-acoustic 
echoes that travel in all directions, in all of the possible propagating modes allowed by 
the theory of elasticity [Ref. 1]. The experimental results of the study showed that the 
concept of a seismo-acoustic sonar for the detection of buried ordnance was scientifically 
feasible. 

During the course of this first investigation into the use of a seismo-acoustic sonar 
for buried ordnance detection, certain key issues requiring resolution were identified. 
These issues need to be solved before the overall concept will have sufficient credibility 


to progress onto naval system development. These issues included [Ref. 1]: 


e Development of sources for selective, or discrete-mode, excitation of seismic 
interface waves. 

e Development of techniques for selective reception of seismic target echoes. 

e Adaptive methods for excitation and reception of optimal modes. 

e Measurement, analysis, and utilization of real mine target strengths. 

e Incorporation of new discovered cetacean sonar protocols and signal physics. 


e Re-examination of vector signal processing technology. 


The purpose of the present research is to address and attempt to resolve some of 
these key issues, specifically, to initiate development of a discrete-mode source and to 
further refine the vector signal processing techniques. The results of this research should 


bring the concept of a seismo-acoustic sonar closer to field implementation. 


Bb. RESEARCH MOTIVATION 


Because of the series of events that transpired during the Persian Gulf War, a 
great deal of attention has been directed towards “main-streaming” mine warfare in the 


US military. The use of mines in the battlefield poses a serious threat to many different 


military operations, and thus significantly impacts the United States’ ability to execute its 
national security objectives. The tactical use of mine warfare directly influences the 
shape of the battlespace and may exclude our military forces from areas in which they 
our called upon to serve. In order to accomplish its military objectives, US forces need 
effective methods to detect and counter enemy mines. As Admiral Boorda stated in his 
1995, white paper on mine warfare, “This warfare area must take its proper place among, 
as well as influence, future research and development, acquisition, and modernization 
plans and programs’ [Ref. 2]. 

In terms of cost-effectiveness and operational and psychological impact, mines 
are perhaps the single most attractive weapons available for disrupting the battlespace. 
For the foreseeable future, we must expect increases in both the lethality and number of 
mines available. Therefore, more advanced technologies must be pursued to make the 
detection, classification, and removal of mines and other buried ordnance items quicker, 


safer, and more cost-effective. 


c. CURRENT AND EMERGING MINE WARFARE TECHNOLOGY 


At the present time (1998) there is no single method currently available that can 
detect buried mines in shallow water sediments, through the surf zone, onto the beach and 
beyond. All the varying buried mine amphibious countermeasure operations can be 
effectively grouped into either breaching for amphibious assaults or post conflict 
clearance operations. Breaching operations involve the construction of safe paths 
through a minefield during combat, while clearance operations deal with eliminating 
unexploded ordnance from large areas for occupation and use. Accurate detection and 
identification are essential for both types of operations. For combat breaching it is 
important to quickly know where the minefields are located and what their density 1s, 
while ordnance clearance operations need to deal with each and every mine, on a non- 


real-time basis. 


1. Current Mine Warfare Technology 


a) Shallow Water Mine Warfare Technology 


From the high-water mark seaward Explosive Ordnance Disposal (EOD) 
divers, Naval Special Warfare (NSW) combat swimmers, Marine Corps reconnaissance 
swimmers, and the U.S. Navy marine mammals, plus a variety of towed sweep gear 
provide the bulk of the available shallow water mine countermeasure systems. 

EOD and NSW personnel serve as an extremely effective neutralization 
component, but provide only minimal detection capability. Sonar and laser detection 
systems such as the Coastal Battlefield Reconnaissance and Analysis system (COBRA), 
the Remote Minehunting System (RMS), and the Airborne Laser Mine Detection System 
(ALMDS) are limited to detection of floating mines, shallow water moored mines, and 
exposed bottom mines. Of all the sonar and laser technologies currently employed for 
shallow water detection only one system 1s truly effective in the detection of buried 
mines, namely the marine mammals. Even though great strides have been made in 
increasing the effectiveness of marine mammal systems, search rates are low and the 


mammals have difficulty working in the surf zone or in warm water. 


b) Land Warfare Technology 


From the high-water mark landward, the currently employed breaching 
and clearance technologies revolve around hand-held probes, metal detectors, ground 
penetrating radar, trained dogs, and mechanical breaching equipment. Each of these 
methods has its unique advantages and disadvantages, but they are in general slow, 
costly, labor intensive, and not optimal for combat operations. 

Manual mine countermeasures provide the most reliable and consistent 
method of ordnance location, classification, and identification. But manual mine 
clearance is slow and dangerous, with an injury rate of nearly one technician every 1,000 


mines cleared [Ref. 3]. 


“Metal detectors” are only capable of detecting ordnance containing 
ferrous components. They are generally difficult to successfully operate because they 
have high false alarm rates. The accuracy of these systems to detect, locate, and analyze 
target returns is based on operator judgment and experience. The high false alarm rates 
associated with these systems is usually due to the high ferromagnetic soil content in 
battlefield areas, as well as by the presence of man made ferromagnetic false targets such 
as hubcaps, tea pots, etc. [Ref. 4] 

Ground penetrating radar incorporates the same general technology as 
conventional radar, but on a much smaller scale. Ground penetrating radar systems fall 
into four main categories: common short pulse systems, video pulse systems, synthetic 
aperture pulse radar, and frequency modulated-continuous wave radar [Ref. 4]. They all 
use the propagation of electromagnetic energy into the soil for target detection and 
localization. Subsurface sensing radars rely on the strength of the electromagnetic signal 
reflected from the target. The strength of the reflected signal is highly dependent on 
target distance, composition, and orientation, in addition to soil characteristics such as 
layering, presence of boulders, etc. High resolution ground penetrating radar systems 
usually operate in the frequency range of 30 MHz to 2 GHz, which implies short 
operating ranges due to the high attenuation of electromagnetic energy in the ground at 
these frequencies. They also suffer from high false alarm rates. 

Trained dogs provide a high degree of detection capability without the 
requirement for prior visual location. Most explosive sniffing dogs are trained to identify 
only one or two types of explosives; thus a whole squad of dogs would be needed to 
detect all the possible explosives types used in just landmines alone. Additionally, the 
amount of time required to train a dog for buried ordnance detection is extremely high, 
while search rates are low because the dogs get tired and detection capability is reliable 
only over very short ranges. 

Lastly. breaching equipment is the brute force method of land mine 
countermeasures. The use of breaching equipment is the most expedient, but of the low- 
end types of available breaching equipment, rollers rarely achieve more that 70% 


clearance effectiveness and rotary flails or blades less than 80%. Explosive net breaching 


whereby a “fish net” containing shaped charges is rocket fired over an area is more 
effective, but much more expensive. 

Although these current methods offer the greatest assurance that an area is 
safe for use, they are also quite dangerous because they put the operator in close 


proximity to the ordnance. 


2. Emerging Mine Warfare Technology 


Emerging technologies with application to detection and clearance can be grouped 
into the following categories: infrared, microwave, photon (laser linescan) backscatter, 
nuclear or thermal neutron analysis, passive optical imaging, and “new” acoustic and 
seismic technologies. Many of these approaches suffer from large power and high data 
processing requirements, and their effectiveness is highly dependent on the propagation 


characteristics of the medium. 


JB. PROPAGATION MEDIUM CHARACTERISTICS 


To realistically address the issue of buried mine detection in the surf zone and on 
land, much emphasis must be placed on understanding the characteristics of the media in 
which the mine detection system must operate. Depth of signal penetration, resolution 
characteristics, detection capability, and other system parameters are all highly dependent 
upon the medium of propagation. These conditions will usually vary from location to 


location and from season to season. 


1. Soil Characteristics 


In general, the earth’s boundary layers consist of a variety of minerals, rocks, 
organic and inorganic compounds, water, and air. The amount of organic materials 
present is specific to particular areas and geographical regions, but most organic 
materials significant to mine warfare are only present in the first few meters of the upper 


sediment layer, and decrease with increasing depth. Changes in the water content result 


in changes of the dielectric and conduction characteristics of the medium. All these 
varying conditions have significant effects on the propagation of electromagnetic and 
acoustic waves. [Ref. 4] 

Nominally, low frequency electromagnetic waves are less attenuated and result in 
increased depths of penetration. Higher frequency electromagnetic signals are attenuated 
more, but provide better resolution characteristics. Acoustic energy is less attenuated 


than electromagnetic energy in high moisture and densely compacted soil environments. 


Zz Surf Zone Characteristics 


In offshore sediments and in the surf zone, the affect of the water and sediment 
interfaces on energy propagation are quite pronounced. Optical techniques have been 
developed to penetrate the water interface, but cannot penetrate the underlying sediment, 
so are essentially useless in buried ordnance detection. Magnetic techniques are also 
limited because of their low resolution qualities, short operating ranges, and the increased 
use of non-magnetic mines. Radar systems are ineffective in the surf zone because of the 
extremely high attenuation suffered by electromagnetic energy in seawater, penetrating 
no more than the first meter into the water. Acoustic sonar systems are capable of 
penetrating both the water and sediment boundaries, but in the surf zone extensive clouds 
of micro bubbles generated by wave action scatter and attenuate most of the acoustic 
energy. The energy which is left can penetrate the sediments, but 1s limited by Snell’s 
Law, which establishes a “critical grazing angle”, below which sound 1s totally reflected 


within the water column and does not penetrate the sediments. [Ref. 5] 


E. SEISMIC WAVE RESEARCH 


With the limitations imposed by the operating environment to current buried 
ordnance detection systems, it was obvious that another method of energy propagation 
must be identified that could more effectively utilize the soil and sediments as a 
propagating medium. With this in mind, it became necessary to extend research efforts to 


incorporate seismic wave techniques. Seismic surface waves also propagate along both 


water/sediment and air/sediment interfaces and can be easily excited in the environments 
of primary interest for buried mine detection systems. Because of their propagating 
characteristics, the use of seismic interface waves can provide a safe minimum standoff 
distance between the operator and the ordnance. 

There are no systems currently available today that can detect buried mines on 
land and in the surf zone. It is clear that development of such a system would help to 
provide the military with the kind of system needed to make the transitions to the shore, 
imp.. 2d in such documents like “...From the Sea”, physically realizable. The benefits to 
other operations, such as, humanitarian demining could also be significant. 

The most common types of seismic waves are those generated by earthquakes, 
which travel both through the body of the earth and along its surface. Seismic waves 
have also been used extensively in geological prospecting for oil and mineral deposits. 
The next chapter contains an overview of seismic wave theory and specific information 


on interface waves and their application to buried ordnance detection. 


fe SEISMIC WAVE THEORY 


This chapter contains a description of the types of seismic waves that propagate 
through an elastic medium and along its surface. Within this cornucopia of wave types 
exist a category of boundary waves known as Rayleigh and Scholte waves. The general 
characteristics exhibited by these waves that make them of practical use in buried 
ordnance detection systems will be explored as well as their relation to the other seismic 


Waves. 


A. GENERAL SEISMOLOGY 


When one studies the general theory of elasticity, two types of mechanical waves 
can be found to propagate in an unbounded, infinite, homogeneous, elastic medium. 
These waves, one compressional and one shear, are called body waves since they travel 
within the body of the propagating medium. In 1885, the English scientist, Lord 
Rayleigh (J.W. Strutt) [Ref. 6] demonstrated theoretically that there exist waves that can 
propagate over the planar interface between a semi-infinite elastic half-space and a 
rarefied gas. These waves later came to be known as Rayleigh waves. Further research 
lead to the discovery of many other boundary waves. Seismic waves can be grouped into 


one of the following two categories: body waves and surface waves. 


1. Body Waves 


There are two types of body waves. The first and fastest propagating wave is 
called the Primary wave or P-wave. A P-wave is a compressional wave, in which the 
wave motion alternately compresses and dilates the medium. The particles of the 
medium undergo displacements in a direction co-linear with the direction of propagation: 
hence the P-wave is also called the longitudinal wave. P-waves can travel through both 


solid rock and liquid materials with a (longitudinal) propagation speed given by 


4 
B+3u 
=. (2.1) 


where B is the adiabatic bulk modulus, 1 the shear modulus, and p the density of the 
medium. 

The second, and slower, body 1s wave is called the Secondary wave or S-wave. S- 
waves are shear waves, in that, as an S-wave propagates it shears the medium at right 
angles to the direction of propagation. Due to the fact that liquids do not support finite 
shear stresses, S-waves cannot propagate in liquids. S-waves travel with a speed given 


by 
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In general, when an external force deforms a solid, both P-waves and S-waves 


will normally be produced. When a wave of either type interacts with an interface 


boundary, reflected waves of both types can be produced. 


2. Surface Waves 


When a solid has a fluid surface (air or water), or when there is an interface 
boundary between two media or between layers of a single medium, surface waves are 
normally present. Surface waves are somewhat similar to gravitational ocean waves for 
which most of the wave motion is located near the surface and decays exponentially with 
increasing distance into both the elastic substrate and the overlying fluid, although the 
mechanisms are different. 

Surface waves can also be divided into many types. One type of surface wave, 
which propagates in the duct between two sedimentary boundaries, is called a Love 
wave. Love waves induce side-to-side particle motion in a horizontal plane, but at right 
angles to the direction of propagation. A second type and more prevalent surface wave is 
the Rayleigh wave. Rayleigh waves, which occur at the planar interface between a semi- 
infinite elastic half-space and a gas, have both vertical and horizontal particle motion in a 
vertical plane oriented along the direction of propagation. Rayleigh waves propagate 
with a phase velocity slightly less than the bulk shear wave velocity in the elastic 
medium. Another type of boundary wave, analogous to the Rayleigh wave, is the Scholte 
wave, which occurs at the interface between a liquid and an elastic half-space such as the 


seafloor. Scholte waves are characterized by a phase velocity that is also slightly slower 
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than the shear wave speed in the solid. Essentially, Rayleigh and Scholte waves are very 
similar, the only difference being that the liquid media component of Scholte waves is 
more pronounced than the gas media component of Rayleigh waves. The interface wave 
that propagates along the boundary created by two elastic half-spaces is called a Stoneley 
wave. Stoneley waves travel with a velocity that is slower than the shear wave velocity 


in either media. 


B. RAYLEIGH WAVES 


In an isotropic solid the stress-strain relations have the following form 

Ox, =AG + 2ye,,, Oy =AO+ 2Ne,,, 6, = AO + 2ye., 2.3) 
Oy, ~ HE, Oy, ~ HE), , O., ~ HE. - 

where A and pt are the Lame constants, 0 the dilatation, or volume strain, and o, and ¢, 

are the stress and strain, respectively, where i denotes the direction of the stress or strain 

and / defines the plane over which it is applied. Using the elastic relations from (2.3), the 


equations of motion for an isotropic solid (neglecting body forces) can be written as 





pas =(A +p) + pa, (2.4) 


where # represents the particle displacement and p the density of the medium [Ref. 7]. 
These equations of motion can be shown to lead to the propagation of the two types of 
body waves with phase velocities given by (2.1) and (2.2). 

If we consider the propagation of a plane wave through an elastic medium with a 
free boundary, such that the boundary is free from stress, we are lead to the wave 


equations 





2 2 
Fy, Ae 5) 
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where the potentials and y were introduced to separate the effects of compression and 
rotation in the medium [Ref. 7]. The simultaneous solution of (2.5) for a disturbance that 
is localized at the boundary and decays with increasing depth requires that the following 


relation be satisfied 


Ko — 8x! + (24 — 16y’)k? + (16y? — 16) = 0. (2.6) 


1] 


This relation is known as the Rayleigh equation, where « is the ratio of the Rayleigh 
wave speed to the shear wave speed, normally 0.9, and y’ is the ratio of the shear wave 


speed to the longitudinal wave speed. y’ can also be expressed solely in terms of the 


Poisson’s ratio, namely 


SG rns eee eV 
Ge a Qu 2-2Vv’ (2.7) 





where v is the Poisson’s ratio. Thus, to solve the Rayleigh equation, which is cubic in «’, 
requires that only the Poisson’s ratio for the propagating medium be known. Therefore 
the Rayleigh wave speed in an infinite isotropic solid is independent of frequency and 
depends only on the elastic constants of the medium. The Rayleigh wave corresponds to 
the root of the Rayleigh equation that lies in the (real) interval (0,1). It can be shown that 
for values of the Poisson’s ratio in the range 0 - 0.5 the Rayleigh equation has one and 
only one such root. Hence Rayleigh waves (R-waves) traveling at the boundary of an 
isotropic, homogeneous medium are dispersionless. [Ref. 7,8] 

Figure 2.1 shows the variation in the wave speed ratios as a function of the 
Poisson’s ratio for an isotropic, homogeneous solid substrate. From the figure we can see 


that Rayleigh waves propagate at almost the same speed as shear waves, varying from 
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Figure 2.1. Variation in Wave Speed Ratios as a Function of Poisson’s Ratio 
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approximately 0.87 to 0.96 of the shear wave speed. Appendix Al contains the 
MATLAB code used to generate Figure 2.1 and Appendix A2, the code used to evaluate 
just the Rayleigh wave speed ratio as a function of the Poisson’s ratio. 

Rayleigh waves can be visualized as the superposition of two waves, one 
longitudinal and one transverse, which propagate along the boundary of an elastic half- 
space. The particle motion associated with a Rayleigh wave is normally elliptical in 
nature, with the semi-major axis of the ellipse perpendicular to the boundary and the 
Semi-minor axis parallel to the direction of propagation, depending on the boundary 
conditions. As Figure 2.2 illustrates the particle motion at the surface is retrograde 
elliptical, but shifts to prograde at a finite depth. The shift in particle rotation occurs at a 
depth of approximately 0.24, for solids, and approximately 0.14, for real sediments, 


where A, 1s the Rayleigh wavelength [Ref. 8,9]. 
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Figure 2.2. Elliptical Particle Motion Induced by Rayleigh Wave 


Figure 2.3 displays the variation in the ratio of the longitudinal to vertical particle 


displacements for Rayleigh waves at the surface of an homogeneous, isotropic solid, as a 


function 


of the Poisson’s ratio, while Figure 2.4 depicts how the particle displacements 


vary with depth. 
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Figure 2.3. Ratio of Longitudinal to Vertical Particle Displacements 
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Figure 2.4. Variation of Longitudinal and Vertical Particle Displacements 


(normalized by the surface vertical particle displacement ) 
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Figure 2.5 shows the variation in the stress components that occur with depth for the 


same isotropic conditions. 


ee 





= . 
BORO CC COR OCC OO OE RETO EEE NEES OOOO OS OOnE OEE HH OO SS OSH 00000 OHH OH OH ONE HORST OOS OOH OOCOS OS OODSOOOS CEOS DOES SS 


normalized stress amplitude 


Vo. G2 0.4 0.6 0.8 1 1.2 1.4 1.6 
depth/Rayleigh wavelength 


Figure 2.5. Tensile and Shear Stresses 


From Figure 2.4 we can see that the vertical displacement first increases and 
reaches its maximum at approximately 0.1A,, then exponentially decreases with depth. 
The longitudinal displacement decreases with depth from the outset and changes 
direction at approximately 0.2A,, reaching a secondary maximum at 0.5A,. From Figure 
2.5 one can see that the tensile stress along the direction of propagation also changes 
direction at approximately 0.254, and the other stress components reach their maximums 
at 0.34,, then all exponentially decay with depth. These sets of curves serve to show the 
localization of the Rayleigh wave to a layer approximately two wavelengths thick under 
isotropic conditions. 

Appendices A3-A6 contains the MATLAB code used to generate Figures 2.2-2.5, 
respectively. Additional insight into Rayleigh wave behavior is contained in the 


computer codes and figures making up the rest of Appendix A. 


Ee: 


C. SEISMIC WAVES 


From the literature [Ref. 6,7,10], where solutions to the elastic wave equations 
have been obtained, we can form an idealized picture, as shown in Figure 2.6, of the 


wave system created by a vertical point impulse at the surface ofan ideal, isotropic, 
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Figure 2.6. Wave System from Vertical Point Impulse [Ref. 10] 


substrate medium. From the figure we can see that the wave system is composed of 
three major events corresponding chronologically to the arrivals of the P-wave, bulk S- 
wave, and R-wave. Thus a particle at the surface will first be displaced by the P-wave, 


then by the bulk S-wave, followed by a much larger oscillation due to the passing 


Rayleigh wave. This illustrates that the Rayleigh wave is the most significant disturbance 
along the surface of an isotropic half- space, and at large distances from the source may 
be the only clearly distinguishable wave. [Ref. 10] 

For a disturbance in the form of a circular footing undergoing vertical oscillations 


at the surface of an isotropic, elastic half-space the basic features of the wave field far 


away from the source are shown in Figure 2.7 [Ref. 10]. 
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Figure 2.7. Wave Field from a Circular Footing 


The two body waves propagate radially outward from the source along a hemispherical 
wavefront. The decrease in the amplitudes of the body waves within the medium is 


proportional to 1/r, where r is the radial distance from the source, while along the surface 
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of the half-space the far-field amplitudes decrease as //r’. Rayleigh waves, on the other 
hand, propagate radially outward with a cylindrical wavefront and experience a far-field 
amplitude decrease that is proportional to //7r. For a vertically oscillating point source 
at the surface of an isotropic, elastic half-space, it was determined that the distribution of 
input energy among the three wave types is 67% Rayleigh wave, 26% shear wave, and 
7% compressional wave. [Ref. 10] 

From this brief exposition on seismic wave fundamentals, the features that make 


Rayleigh waves unique (in an isotropic half-space) are that they are: 


e Dispersionless, thus propagate with a phase velocity that is independent of 
frequency. 

e Localized within a layer that is approximately two wavelengths deep. 

e The most significant disturbance along a surface boundary. 

e Easy to excite, with a particle motion that is retrograde elliptical at the 


surface. 


It is these features that make Rayleigh waves ideal for buried ordnance detection 
in the earth and along the beach. Scholte waves exhibit the same general characteristics 
as Rayleigh waves, the major difference being that Scholte waves occur along the 
interface between a liquid and an elastic half-space, thus making them ideally suited for 
surf zone applications. 

It should be remarked that in natural soil, beach, or seafloor sediments, the 
overburden weight of the surface sediments causes the shear modulus to increase with 
increasing depth. This means that the shear velocity increases with depth and that, in 
turn, means we have a condition requiring dispersion. This occurs because the lower 
frequency, long wavelength waves naturally propagate at deeper depths in the sediment, 
where the shear velocity is greater. Thus, low frequency, long wavelength signals 
propagate faster than high frequency, short wavelength signals. 

If an electromechanical source could be built that could selectively excite 
Rayleigh waves, more of the energy used in signal generation could be put on target, and 


less of the projected energy would exist to create unwanted reverberation carried in non- 


Rayleigh wave modes. However, at the target, even pure Rayleigh waves get distributed 
into other non-Rayleigh modes. With the addition of special signal processing techniques 
that exploit the known features of Rayleigh waves from target echoes, we could capitalize 
on discrete-mode sonar methods. The next chapter discusses the previous research 
efforts made in buried ordnance detection using scismic waves and the experimental 


limitations imposed by the systems used. 
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HI. PREVIOUS RESEARCH EFFORTS 


This chapter briefly reviews two previous research projects that explored the 
propagation of energy along an air/sediment interface to detect buried objects. The first 
project was an experiment conducted at the Naval Postgraduate School, investigating the 
applicability of array theory to target localization. The second project, done at the 
Applied Research Laboratories of the University of Texas at Austin (ARL:UT), for which 
this thesis is a direct continuation, developed the theory and signal processing techniques 
for the first phase of the seismo-acoustic sonar system. Other significant prior research 
consists of the dry-land experimental work done at the firm of Bolt, Beranek, and 
Newman (BBN) and theoretical work done at the Coastal Systems Station. These works 


will be cited subsequently. 


A. NAVAL POSTGRADUATE SCHOOL (NPS) 


In 1995, NPS student LT William Stewart conducted research into the feasibility 
of using a phased array to generate surface waves to detect buried ordnance [Ref.11]. His 
research was the continuation of a study conducted by BBN Systems and Technology 
from 1987 to 1992 [Ref. 12]. The source utilized in the NPS study was a four-inch, large 
cone excursion, electrodynamic speaker, whose frequency response for the radiation of 
sound was relatively constant over the range 50 — 7000 Hz. The speaker cone was 
modified by bonding to it a 25-gram mass, cone-tipped plunger, which was used to 
improve ground coupling. See Figure 3.1. 

Three single-axis seismometers were used as ground vibration receivers. Testing 
was conducted in a twelve-foot diameter, four-foot deep redwood tank filled with 
approximately 17 tons of medium grain sand. With a single source and the three 
receivers, LT Stewart was able to successfully demonstrate target localization using 
surface waves. Though no data was reported on the results for target localization using a 
phased array, his study did show that beamforming with a three element linear array of 


surface sources was possible. [Ref. 11] 


Speaker Cone 





Cone-tipped Plunger __§_» 


Figure 3.1. NPS Source Design 


B. ARL: UT 


An investigation into the use of Rayleigh waves for buried ordnance detection 
was conducted by ARL:UT on the beach of Mustang Island, near Corpus Christi, Texas, 
in 1996 [Ref. 1]. Their experimental setup consisted of a source which employed a six- 
inch by eight-inch exciter foot, from which protruded 43, %-inch nails and a receiving 
array composed of three seismometers, as shown in Figure 3.2. 

Their raw, field recorded data did not show the target. But with the use some 
sophisticated vecior signal processing techniques and some coherent subtraction methods, 
operating on the reverberation, the target was detected at its appropriate range and 
bearing. A seismo-acoustic detection model for elastic targets buried in elastic sediments 


was also developed to support the measurements [Ref. 13]. 
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Figure 3.2, ARL:UT Seismic Sonar Experiment 
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Cc COASTAL SYSTEMS STATION 


Some of the benchmark research efforts recently produced by the laboratories of 
the Coastal Systems Station concern the theoretical solutions and numerical modeling of 
acoustic wave scattering from elastic objects buried in underwater sediments and plane- 
stratified poroelastic media [Ref. 14,15]. Objects located in fluid mediums produce 
resonant scattering behaviors. These same objects when buried in ocean sediments 
display a resonant scattering behavior which is modified by the proximity to interfaces, 
sediment attenuation, and sediment loading, thus information obtained from the scattering 
field can be used for identification as well as detection. Further experimental verification 
of these solutions could have significant impact on the capabilities of seismo-acoustic 


sonar systems. 


D. EXPERIMENTAL LIMITATIONS 


What 1s common to the aforementioned experiments is the use of vertical source 
excitation, which generates not only Rayleigh waves, but also undesired P- and S-waves. 
Thus, the Rayleigh wave signals of interest were contaminated by body wave signals. To 
satisfactorily localize the target, in the ARL:UT experiments, this unwanted reverberation 
had to be removed during the signal processing by coherent subtraction of the reverberant 
background field in the absence of the target from the field with the target present. 
Vector signal processing was also used for the first time, in the ARL:UT experiments, to 
bring the target echoes above the reverberation. This technique capitalizes on the known, 
90-degree phase difference between vertical and horizontal components in propagating 
Rayleigh waves. But since coherent subtraction is not a feasible signal processing option 
in a operationally fielded system, another solution must be identified. The research team 
at ARL:UT theorized that considerable reverberation suppression could be achieved by 
utilizing an excitation source that could selectively excite Rayleigh waves [Ref. 1]. They 
also experimented with a two-degree of freedom source, which preferentially excited 
Rayleigh waves that propagated to tens of meters in range. Unfortunately, this source 
produced a signal too long in time duration for sonar application. 

The use of ultrasonic Rayleigh waves in non-destructive testing is a well- 
documented science, but selective excitation of low frequency Rayleigh waves in 
terrestrial media has not been previously attempted. The remainder of this thesis 
describes the development and testing of a discrete-mode, electromechanical, Rayleigh 


wave source for a Selsmo-acoustic sonar. 
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IV. DISCRETE-MODE SOURCE DEVELOPMENT AND TESTING 
(GROUND PHASE) 


This chapter documents the initial design, construction, and testing of the first 
four generations of the discrete-mode, Rayleigh wave source. The associated 
experimentation involved testing in both a laboratory environment and various field 
environments. This portion of the research, known as the Ground Phase, was used for 
proof-of-concept testing and qualitative determination of the relations between input 


signal parameters and the resulting surface wave field. 


A. DESIGN CRITERIA 


The first step in constructing a discrete-mode, electromechanical, Rayleigh wave 
source for a seismo-acoustic sonar system was to establish a set of guiding design 
criteria. Based on the known features of Rayleigh waves and the previous research work 
completed, it was decided that minimum requirements that a source needed to meet were 


that it: 


e Have sufficient degrees of freedom (2) to drive the ground in elliptical orbits. 

e Have good ground coupling in both the horizontal and vertical directions. 

e Generate Rayleigh waves with sufficient amplitudes to allow detection at a 
range of at least 10 m. 

e Provide a repeatable waveform capable of high repetition rates over long 
periods of time. 


e Be able to operate over a range of low frequencies in the 20 — 200 Hz band. 


The most important of these requirements is that the transducer oscillate the 
ground in elliptical orbits, so that as much as possible of the input energy be converted to 
discretely excite Rayleigh waves. 

All seismic transducers consist of two basic components: a vibration source and a 


source-to-ground coupling mechanism. The vibration source transfers the input electrical 


energy into output mechanical energy. The source-to-ground coupling mechanism 
provides the means by which the output mechanical energy is transmitted into the 
propagating medium. A constraint on the source design was that the ground coupling 
mechanism be as light weight as possible, while at the same time being as stiff as 
possible. This was desired so most of the mechanical force generated by the source 


would be used to move the ground and not just the source. 


B. CONCEPT SOURCE 


With the design criteria established for a suitable source, the next step was to 
design and construct laboratory prototypes that could be tested to examine the effects of 
varying source parameters such as frequency and driving voltage, on the resulting surface 
wave field. This was deemed a necessary first step because the results obtained would 
help in selecting an appropriate vibration source and designing the ground coupling 
mechanism for a field test source. It was decided to construct the concept source out of 
on-hand laboratory equipment and to conduct the required testing in the twelve-foot 


diameter, four-foot deep redwood, sand filled tank used in the previous NPS experiment. 


1. Vibration Source 


The Model F5B Electromagnetic Shaker (or Mini-Shaker for short) produced by 
Wilcoxon Research was chosen as the vibration source for proof-of-concept testing. The 
Mini-Shaker is an electrically driven force generator, with a constant (blocked) force 
output of 0.2 lbs (0.89 N) over the frequency range of 100 — 1000 Hz (See Appendix B1 
for manufacturer specifications). The design configuration of the Mini-Shaker allows it 
to be stud-mounted in any position without external support or shaft alignment 


considerations. 
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2 Ground Coupling Mechanism 


The source-to-ground coupling mechanism, as shown in Figure 4.1, was 
constructed from two pieces of aluminum. Aluminum was chosen because of its high 


tensile strength and low mass. Appendix C1 contains the actual design drawings. 





Figure 4.1. Picture of Concept Rayleigh Wave Source (scale in inches) 


The bottom piece was constructed from a hollowed out aluminum block. When 
buried in the sand, the inside surfaces of the block provide the required vertical and 
horizontal coupling to the medium. The top piece was manufactured to hold two Mini- 
Shakers and a tri-axial accelerometer (Accelerometer specifications contained in 
Appendix B2). Mounting studs were placed in the top piece so as to align the Mini- 


Shakers at a 45 degree angle with the horizontal surface plane of the sand. This specific 


placement of the Mini-Shakers was chosen so that any two-dimensional oscillatory 


motion could be excited by individually controlling the phase and amplitude of each 


shaker. 


3. Experimental Results 


The test equipment setup as described in Appendix D1 was used for the concept 
source testing. The source and seismometer receivers were positioned in the test tank as 
shown in Figure 4.2 and Appendix DI. This setup allowed for the greatest separation 
between the source and the tank wall and the farthest receiver and the tank wall. The 
seismometers were tri-axial, moving coil geophones, which measure the velocity of 


ground movement. 





Figure 4.2 Concept Source Testing Experimental Setup 


The first experiment was intended to qualitatively determine the effects of phase 
differences between the Mini-Shakers on the axial surface wave field, as a function of 
frequency. These results were characterized by the difference in the seismometer outputs 
for the axial and vertical directions. The transverse direction was not evaluated because 
the Rayleigh wave motion of interest exists in the vertical plane along the direction of 
propagation. Continuous sinusoids of frequencies 100, 150, 200. 250, 300, and 500 Hz 
were chosen as the test values. Since the left most Mini-Shaker in Figure 4.2 was used as 
the reference shaker, phase differences were chosen to correspond to the following six 


Cases: 


e Case 1 - Right shaker lagging by 25 degrees. 
e Case 2 - Both shakers in phase. 

e Case 3 - Right shaker leading by 25 degrees. 
e Case 4 - Right shaker leading by 45 degrees. 
e Case 5 - Right shaker leading by 90 degrees. 
e Case 6 - Right shaker leading by 135 degrees. 


In al] test cases except the second, which produced exclusively vertical motion, the 
motion of the source was elliptical as was seen by the lissajous pattern of the radial and 
vertical components of the accelerometer. Appendix E contains the qualitative 
experimental results for all phases and frequencies evaluated. The results of this first 
experiment showed that there was a great deal of variation in the axial surface wave field 
with both source frequency and relative phase. 

A second experiment was conducted to assess the effects of not only varying the 
frequency and phase, but also the driving voltages of the Mini-Shakers. This second 
experiment further served to show that the axial surface wave field, even in the nearly 
homogeneous, isotropic conditions of the sand tank, was a highly complex outcome of all 
variable source parameters. 

Even though surface waves were generated in all test cases, the specific 
combination of source parameters needed to produce Rayleigh waves of sufficient 


amplitude and duration at the desired ranges for buried ordnance detection had yet to be 


Zo 


determined. It was evident from the initial testing that the appropriate combination 
would not be found in the laboratory. This is because, not only is the correct 
combination highly dependent on the individual properties of the vibration source and the 
ground coupling mechanism, but also on the local properties of the propagating medium 
at the time of testing. Furthermore, the limitations imposed by the tank, such as 
reflections from the tank walls and short ranges, made it impossible to quantify the 
source. To proceed any further in developing a discrete-mode, Rayleigh wave source 


would require designing a source that could be tested under actual field conditions. 


C. FIELD TEST SOURCE —- MOD I 


To begin constructing a discrete-mode, Rayleigh wave, source that could be used 
in the field required that a suitable vibration source be identified. Once the vibration 
source was chosen, then the source-to-ground coupling mechanism could be 
appropriately designed and constructed. The vibration source needed to be rugged 
enough to withstand the beach and surf zone environment, yet inexpensive enough so that 
multiple pairs could be acquired, if arrays of paired vertical and horizontal sources could 


eventually be utilized. 


ie Vibration Source 


A test source that satisfied most of the necessary prerequisites was found 
manufactured by Aura Systems Inc., under the nomenclature of a Bass Shaker (See 
Appendix B3 for manufacturer specifications) [Ref. 16]. The Bass Shaker is a moving 
magnet transducer, intended to be used as a low frequency vibration source for a car 
stereo system. They are normally mounted under car seats or in door panels to provide 
listeners with the sensation of sound, without moving air like an ordinary speaker. The 
Bass Shaker generates a nominal force output of 10 Ibf (44.5 N), independent of its 


mounting orientation. 
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De Ground Coupling Mechanism 


Since the Bass Shakers were considerably different in size and shape from the 
Mini-Shakers used previously, the source-to-ground coupling mechanism had to be 
redesigned. It was decided that the Bass Shakers should still be positioned at a 45- 
degree incline to the horizontal plane to provide the necessary degrees of freedom to 


produce all variations of elliptical motion. As shown in Figure 4.3, the ground 





of ig 2 at 


Figure 4.3. Picture of Field Test Source — Mod I 


coupling mechanism was constructed from two pieces of aluminum. The top piece 
serving as a support brace for the upper part of the Bass Shakers and as a mounting 


surface for the electrical (BNC) connectors. The bottom piece serves as the attachment 


point for the lower part of the Bass Shakers and provides the coupling mechanism to the 
medium through two rows of machine screws. Appendix C2 contains the Mod I 
technical drawings. A tri-axial accelerometer was positioned in the center of the base to 
monitor the source motion, and there was also an accelerometer mounted on the base of 


the support brace. Refer to Appendix B2 for accelerometer specifications. 


SF Experimental Results 


Before commencing any field experimentation, the new Rayleigh wave source 
was put through some operator familiarization and initial shakedown testing in the sand 
tank laboratory. Using the test equipment setup as described in Appendix D2, the same 
experiments were conducted as with the Mini-Shakers. These tests provided further 
verification that the axial surface wave field was a strong function of the input source 
parameters. 

With the knowledge gained from the laboratory, it was clear that the first goal of 
any field experimentation s!_ .Id be to assess the Mod I source’s ability to generate 
surface waves. To do this required finding a nearby location that was relatively free from 
any underground obstructions that would scatter the outgoing surface waves and make 

detection more difficult. The choice was made to use the NPS softball field as the first 
| site for the field experiments. Progressing onto field experimentation from an equipment 
point of view proved to be quite a challenge at first, but as the research progressed, 
equipment difficulties were quickly overcome. 

Once at the softball field, the lab equipment was setup as described in Appendix 
D3. Figure 4.4 illustrates the initial experimental layout. The initial distances were 
chosen based on expected Rayleigh wave speeds on the order of 100 m/sec and the 
maximum length of electrical cable available. Due to the hardness of the ground, 
trenches for the two rows of machine screws had to be dug in order to get the base of the 
source flush with the ground. After all the equipment was setup, attempts were made to 
detect the source ‘notion above the background noise at the two seismometers located at 
ranges of 15.2 m (50 ft) and 30.5 m (100 ft). After many unsuccessful attempts to detect 


Rayleigh waves at these ranges. the decision was made to move back to ground zero. 


The two seismometers were repositioned at distances of 1 m and 5 m from the source. 
From the variation in source parameters attempted, it was clear from the accelerometer 
output that the Mod I coupling mechanism had good horizontal coupling with the ground, 


but poor vertical coupling. 
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Figure 4.4. Field Source Testing Experimental Layout 


At a distance of 1 m, the surface waves generated by the source motion were 
clearly visible, but at 5 m could not be separated from the background noise. This was 
clearly unacceptable, so to increase the chance of detection, improvement in the signal to 
noise ratio was required. First attempts were made by increasing the duty cycle, from 5 
cycles to 20 cycles in increments of five. The increase in signal duration was clearly 
visible at 1 m, but at 5 m the signal was still buried within the background noise. The 
next effort was to increase the drive voltages of the two Bass Shakers. Starting at 5 


Vrms, the drive voltage was increased in increments of five. until the current limit of the 


o>) 
3 


power amplifiers was reached at approximately 39 Vrms. Even at 39 Vrms, the surface 
wave motion generated by the source was barely discernable at 5 m. 

As a result of the first field experiment the following conclusions were drawn. 
First, to generate Rayleigh waves of sufficient amplitude to be detected at ranges of at 
least 10 m would require that the Bass Shakers be driven at the upper end of their linear 
response voltage limit. Second, that the ground coupling mechanism needed to be 
redesigned so as to enhance the vertical coupling. Furthermore, since there were so many 
channels of data that needed to be observed and ultimately processed, a method of 


recording all the field data was also required. 


D. FIELD TEST SOURCE — MOD II 


Using the lessons learned from the first field experiment, the next major 


milestones were to: 


e Find the operating limits of the Bass Shaker. 
e Modify the ground coupling mechanism to increase the vertical coupling. 
e Incorporate into the field test equipment a method of obtaining multiple data 


channels. 


1. Vibration Source 


The manufacturer specifications that came with the Bass Shakers did not indicate 
any electrical operating limits. The installation guide (Appendix B3) did provide an 
advised maximum continuous power rating of 25 W per shaker. With an input 
impedance of 4 92 and a 25 W power rating, the maximum continuous drive voltage is 
only 10 Vrms. These ratings apply only to its intended use as a car stereo component and 
not for use as a vibration source for a seismic transducer. Consultation with the 
manufacturer provided no further information. 

Without any formal documentation available, the only way to adequately assess 


the operating limits of the Bass Shakers was to see under what operating conditions 
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would failure occur. Sets of failure tests were conducted using the test equipment setup 
as described in Appendix D4. The goal of the failure testing was to estimate the amount 
of signal distortion as a function of frequency and drive voltage. 

The HP35665A Signal Analyzer was used to simultaneously monitor the time 
series and total harmonic distortion of the Bass Shaker using the output signal of a uni- 
axial accelerometer mounted on the shaker case (See Appendix B4 for accelerometer 
specifications). The time series of the accelerometer output was used to visually detect 
any signs of failure or non-linearity as evidenced by changes in signal shape and/or 
duration. An estimate of the signal distortion over the first ten harmonics of the 
fundamental frequency, i.e. the drive frequency, was obtained using a built-in machine 
function. 

The first set of failure tests were conducted by driving the shaker with 20 cycles 
of a sine wave, with a 2.00 sec repetition period, for fundamental frequencies of 100. 
150, 200, 500, and 1000 Hz. While holding the frequency fixed, drive voltages of 10.1, 
20.1, 30.1, and 38.6 Vrms were used, with 38.6 Vrms being the maximum voltage 
attainable without exceeding the current limit of the power amplifier. Figure 4.5 is a 
graphical representation of the signal distortion as a function of frequency and drive 
voltage for a 2.00 sec repetition period. Even though the overall signal distortion was 
low, the accelerometer time series for the 1000 Hz drive frequency showed evidence of 
non-linear behavior. The leading and trailing edges of the accelerometer signal were not 
as detailed as compared to the output signal from the function generator. There was also 
evidence of pulse elongation in the signal. The non-linear behavior was attributed to 
internal heating effects, as evidenced by noticeable case heating observed, but the shaker 
did not fail. 

To stress the Bass Shaker further toward failure required that a more powerful 
amplifier be utilized. To hopefully alleviate the distortion caused by internal heating the 
repetition period was increased to 5.00 sec. A second set of failure tests were conducted 
with the new repetition period using a higher rated amplifier. During this second round 
of testing, the test shaker failed. The shaker was disassembled to locate the cause of the 


failure, but the internal design of the Bass Shaker only made the explanation speculative. 
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Figure 4.5. Signal Distortion for 2.00 sec Repetition Period 


Since the cause of the failure of the first Bass Shaker was unknown, a second 


shaker was setup to complete the evaluation with the new amplifier. Figure 4.6 displays 


the results of the second set of tests. 
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Figure 4.6. Signal Distortion for 5.00 sec Repetition Period 


Using a repetition period of 5.00 sec seemed to effectively reduce the amount of 
internal heating and resulting distortion. As can be seen from Figure 4.6, the optimal 
drive frequency with respect to minimal signal distortion was 250 Hz. Again since the 
current limit of the power amplifier was exceeded before the shaker failed. it is safe to 
assume that with the current test equipment setup, shaker failure due to input drive 


voltage 1s not a major concern. 


ae Ground Coupling Mechanism 


To enhance the vertical coupling of the source to ground coupling mechanism, 
two end attachments were mated to the original base as shown in Figure 4.7. The 


attachments were constructed to hold three machine screws each, making the same 45- 


ce 
BEI Det, WIG apn te psy 


ie are tree gl Co 
Ate vis sa eS sib a. me cha eso 


> <a y = 
> st 5 
+ b 


{ 


oy eRe: ges an ea xk 


ay ee ree ane 

il es 

- 0" es cp - o- es - rH 

y te © aio pa a ONO gram yn Paves 2. &. 

ww & ma : wth “ y 9 a a ae neal 

: 3 “lew ed ick Zee cr em 

pee a ei ay t TER secs ogee eee 
Pa. oh £ e = @ © igs: wee en rs peed Zs a. ee aie we 7 
PO eee ew eitlen ee: waa Fa Oo a OO A a mn a) er Oe WT Sina ee 


* 
, e . 
» ¢ r 
p v 





Figure 4.7. Picture of Mod II Base Modifications 


degree angle with the horizontal surface plane as the Bass Shakers. It was hoped that 
these modifications to the base of the ground coupling mechanism would increase both 
the horizontal and vertical source coupling. 

The upper support brace for the Bass Shakers was also modified. First, it was 
reduced in size to position the shakers closer to the desired 45-degree angle. Second, 
new BNC mounts were incorporated, making the Mod II source more rugged. Figure 4.8 


is a picture of the redesigned Rayleigh wave source. Appendix C3 contains the technical 


drawings for the modifications. 





Figure 4.8. Picture of Field Test Source — Mod II 


Addition of the two end attachments meant that emplacement of the source would 


first require placing the base into the ground in a similar fashion as the Mod I source, 


then by using a cordless drill boring the six remaining angled machine screws into the 


ground. 


Sy. Experimental Results 


To resolve the data gathering issue, a SD390 Dynamic Signal Analyzer, made by 
Signal Processing Systems, was added to the test equipment setup described in Appendix 
D3. The SD390 is capable of receiving, displaying, and storing eight channels of 
information simultaneously, in real time. The eight channels provided data gathering 


capabilities for the: 


e Radial and vertical components of the base-mounted accelerometer. 
e Longitudinal, transverse, and vertical components of the two tri-axial 


seismometers. 


Field testing of the Mod II source was also conducted at the softball field, using 
the equipment setup as shown in Appendix D5, but the source experimentation did not 
progress much farther than equipment setup. It was clear from the output signals of the 
accelerometer that the source was not exciting horizontal motion equally as well as it was 
exciting vertical motion. Thus generating the required elliptical motion to selectively 
stimulate Rayleigh waves would be nearly impossible. Attempts to replace the source 
were made, but did not significantly improve the operating conditions. 

The results of the second field experiment prompted the need for further 
refinements to the design of the ground coupling mechanism and to rethink the utility of 


the softball field as the test location. 


E. FIELD TEST SOURCE — MOD Ill 


It is, of course, important to have good horizontal and vertical coupling to the 
eround for Rayleigh wave stimulation. Generating elliptical particle motion at the source 


requires moving the medium through a combination of both radial and vertical 


excitations, such that, the semi-major axis of the ellipse 1s normal to the surface boundary 
and the semi-minor axis parallel to the interface in the direction of propagation. 
Therefore, having the correct design for the ground coupling mechanism is a major 


portion of discrete-mode source development. 


1. Ground Coupling Mechanism 


From the second field experiment it was noted that boring the machine screws 
into the ground with the cordless drill seemed to lock the source in place. The problems 
with the ineffective coupling arose from having to dig holes for the pre-positioned 
vertical machine screws. Hence, it only seemed logical that the source base be modified, 
so that all the screws could be positioned by screwing them into the ground. The Mod III 


base as shown in Figure 4.9 is the result. 
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Figure 4.9. Picture of Mod HI Base Modifications 
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The Mod III base modifications consisted of the addition of the two larger end 
attachments to the original base and the switch from 3 % inch machine screws to 6 inch 
machine screws. Reuse of the original base only served to reduce the turn-around time 


between field tests, since a whole new base did not have to be machined. 


pie Experimental Results 


Due to the equipment difficulties that arose from using the softball field as the test 
location, it was decided that any following ground phase experimentation be conducted 
from the lab in Spanagel Hall. Appendix D6 contains the new test location layout. The 
same equipment setup as described in Appendix D5 was used. 

With the necessary modifications to the ground coupling mechanism completed, 
the next field experiment was used to ensure that the Mod III base provided the required 
horizontal and vertical coupling for Rayleigh wave stimulation. This phase of the source 
testing consisted of verifying, that with the correct phase settings for each of the Bass 
Shakers, both radial and vertical source motion could be individually stimulated. The 
stiffness of the ground to the source motion was quite noticeable, it was only through the 
manual manipulation of both phase and amplitude, that individual excitation of radial and 
vertical motion was possible. 

The fourth, and final, ground phase field experiment was used to investigate the 
axial surface wave field effects for various drive voltages and relative phases of the Bass 
Shakers. Refer to Appendices D5 and D6 for the corresponding equipment and 
experimental setups. The effects to the surface wave field were characterized by 
measuring the coherence between the radial and vertical signal components at the furthest 
seismometer. Coherence on the order of unity would indicate a high degree of reliability 
in the magnitude and phase information of the received signals. The coherence was 
determined using a built-in processing option of the SD390 Signal Analyzer, using the 
average of ten received signals from a 250 Hz, 20 cycle sine wave burst, with a 5.00 sec 
repetition period. Driving voltages of 20, 40, 60, and 86 Vrms and relative phases of 0. 


30, 60, 90, 120, 150, 180 degrees were chosen as the test values. 
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excitations, such that, the semi-major axis of the ellipse is normal to the surface boundary 
and the semi-minor axis parallel to the interface in the direction of propagation. 
Therefore, having the correct design for the ground coupling mechanism is a major 


portion of discrete-mode source development. 


| Ground Coupling Mechanism 


From the second field experiment tt was noted that boring the machine screws 
into the ground with the cordless drill seemed to lock the source in place. The problems 
with the ineffective coupling arose from having to dig holes for the pre-positioned 
vertical machine screws. Hence, it only seemed logical that the source base be modified, 
so that all the screws could be positioned by screwing them into the ground. The Mod III 


base as shown in Figure 4.9 is the result. 
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Figure 4.9. Picture of Mod III Base Modifications 
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The Mod III base modifications consisted of the addition of the two larger end 
attachments to the original base and the switch from 3 % inch machine screws to 6 inch 
machine screws. Reuse of the original base only served to reduce the turn-around time 


between field tests, since a whole new base did not have to be machined. 


ae Experimental Results 


Due to the equipment difficulties that arose from using the softball field as the test 
location, it was decided that any following ground phase experimentation be conducted 
from the lab in Spanagel Hall. Appendix D6 contains the new test location layout. The 
same equipment setup as described in Appendix D5 was used. 

With the necessary modifications to the ground coupling mechanism completed, 
the next field experiment was used to ensure that the Mod III base provided the required 
horizontal and vertical coupling for Rayleigh wave stimulation. This phase of the source 
testing consisted of verifying, that with the correct phase settings for each of the Bass 
Shakers, both radial and vertical source motion could be individually stimulated. The 
stiffness of the ground to the source motion was quite noticeable, it was only through the 
manual manipulation of both phase and amplitude, that individual excitation of radial and 
vertical motion was possible. 

The fourth, and final, ground phase field experiment was used to investigate the 
axial surface wave field effects for various drive voltages and relative phases of the Bass 
Shakers. Refer to Appendices D5 and D6 for the corresponding equipment and 
experimental setups. The effects to the surface wave field were characterized by 
measuring the coherence between the radial and vertical signal components at the furthest 
seismometer. Coherence on the order of unity would indicate a high degree of reliability 
in the magnitude and phase information of the received signals. The coherence was 
determined using a built-in processing option of the SD390 Signal Analyzer, using the 
average of ten received signals from a 250 Hz, 20 cycle sine wave burst, with a 5.00 sec 
repetition period. Driving voltages of 20, 40, 60, and 86 Vrms and relative phases of 0. 
30, 60, 90, 120, 150, 180 degrees were chosen as the test values. 
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The quantitative results of the experiment are shown in Figure 4.10. As can be 
seen from the figure, the coherence between the radial and vertical components of the 
received signal increases as drive voltage increases, being nearly unity over most phase 
angles for a drive voltage of 86 Vrms. This final ground phase experiment illustrates 
that, the greater the input energy into the medium, the better the stimulation of reliable, 
repeatable, surface wave fields. The reason for the low coherence at 120 degrees is 


unknown. 
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Figure 4.10. Coherence Measurements as a Function of Voltage and Phase 


The time and effort spent in conducting the Ground Phase research provided a 
great amount of detail about the behavior of the field test sources and the associated 
surface wave fields. The choice of the Bass Shakers as the vibration source proved to be 
quite effective in generating the kinds of motion desired for discrete stimulation of 


Rayleigh waves. 


42 


Since the development of the seismo-acoustic sonar concept, the ultimate goal of 
the combined research efforts was to develop a mine detection system that was capable of 
locating buried ordnance, both in the surf zone and on land. Clearly, the surf zone is the 
more complicated environment to operate in, with regards to both the physics and the 
equipment required. The next chapter describes the development and experimental 


research conducted in this new environment. This portion of the research was identified 


as the Beach Phase. 





V. | DISCRETE-MODE SOURCE DEVELOPMENT AND TESTING 
(BEACH PHASE) 


This chapter covers the beach and surf zone development and experimental 
testing of the discrete-mode, seismic wave source. Included within the chapter are 
individual sections documenting the different phases of the experimentation. The first 
section briefly addresses the final modifications made to the source-to-ground coupling 
mechanism. The next section covers the ambient background noise analysis and ground 
truthing experiments conducted at the test site location. The remaining sections describe 
the research efforts made to assess the ability of the test source to selectively excite 
Rayleigh waves. 

During the ground phase experimentation, research efforts were mainly directed 
at investigating the operating characteristics of the field test sources, in order to find the 
most effective source-to-ground coupling design. Secondly, it was also concerned with 
qualitatively assessing the effects of various source parameters on the axial surface wave 
field. The ground phase research provided the necessary background information that 
was needed to successfully interpret the beach phase experimental results. Progressing 
onto beach and surf zone field experimentation was a major milestone in furthering the 
development of a discrete mode selsmo-acoustic source and towards a greater 
understanding of the complex physics involved. 

Using the lessons learned from the ground phase experimentation with respect to 
equipment transportation and logistics, it was evident that working out the bed of a 
pickup truck would not suffice for the beach phase experimentation. Thus, it was 
necessary to find a research platform that was large enough to transport all the equipment 
to the test site and that could provide the necessary electrical power. The mobile home 
owned by the Oceanography Department was the logical and most readily available 


solution. 


A. FIELD TEST SOURCE — MOD IV 


From the results obtained during the last set of ground phase experiments, the 
Mod III source-to-ground coupling mechanism, as shown in Figure 4.9, seemed to 
provide the necessary coupling in both the horizontal and vertical directions. As a result 
of its successful testing it was decided that this design be used as a model for the 
fabrication of the final ground coupling mechanism to be used for the beach phase 


experimentation. Figure 5.1 is a picture of the Mod IV ground coupling mechanism. 


Appendix C4 contains the technical drawings for the Mod IV base. 
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Figure 5.1. Picture of Mod IV Ground Coupling Mechanism 
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To help protect the sources from the harsh beach environment, all joints and 
access holes were filled with silicone. Furthermore, a waterproof accelerometer was also 


added (See Appendix B5 for manufacture specifications). 


B. BACKGROUND NOISE AND GROUND TRUTHING EXPERIMENTS 


Because of the highly complicated nature of the surface wave field, it was 
necessary to obtain a sufficient knowledge base of the ambient acoustic properties of the 
test medium. To obtain the required information meant that ambient background noise 
analyses and ground truthing experiments be conducted at the proposed source testing 
site. The test site chosen was a stretch of beach on the Monterey Bay coastline located 
behind the NPS Beach Research Facilities. The section of beach front area used, as 
shown in Figure 5.2, was approximately 80 feet wide from the surf zone to the sand 


dunes at low tide with a longitudinal extent far exceeding that required for the source 













testing. 
=~ ~ tal. / | ; 
Soa ; Wap. 7 
ab $ cccennne N ; 
1 y , 
Se Monterey Bay iP 
RK Ps 
a MONTEREY 
FISHERMAN'S 
WHARF . 
- \ 
Beach Phase 
eS € Testing Sit 
, cusrom “—~s>=~ esting Site 
HOUSE PLAZA = 


a) MONTEREY > 
2 


STATE BEACH 


1 
~ \ 
as lies 
lg | 
0 
Pd 


oo 


ee, 
43 


Figure 5.2. Beach Phase Testing Site 


Map copyrighted 1994 by the California State Automobile Association. Reproduced by permission. 
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PF: Background Noise Analysis 


In order to make source testing and data interpretation as simple as possible it was 
desired that the source be operated in a frequency range for which competition with the 
background noise would be minimal. Background noise samples were taken at the test 
site using the seismometers as the receivers. Spectral analysis of the noise samples taken 
showed that the strongest noise components resided between 5-20 Hz as shown in Figure 
5.3. Since these frequencies were well below the envisioned operating frequencies for 
the Bass Shakers filtering out the undesired background noise from the input signals 
would be relatively simple. Appendix F contains the MATLAB programs used to 


conduct the spectral analysis. 
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Figure 5.3. Ambient Background Noise 


jas Ground Truthing 


To obtain an accurate picture of the ambient acoustic properties of the test site a 
Series of ground truthing experiments were conducted to estimate the group velocity of 
the Rayleigh wave. As was shown in section II C, Figure 2.6, a vertical point impulse 
generates a wave system comprised of P-, S- and R-waves, as well as other wave types 
not shown. 

To simulate a vertical impulsive source a bowling ball with a self-contained uni- 
axial accelerometer was dropped from an average height of 45 inches, as shown 1n Figure 


5.4. A bowling ball was used because it is spherically symmetric and impacts the ground 


48 


the same way on each drop. When the bowling ball impacts the ground several 
radially outward propagating wave systems are created. The impact of the bowling ball 
is sensed by the accelerometer, which triggers a multi-channel data acquisition system to 
record the response of a pre-positioned seismometer. Appendix B6 contains the 
manufacturer specifications for the uni-axial accelerometer and Appendix D7 the wiring 


diagram for the equipment setup. 
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Figure 5.4. Picture of Initial Ground Truthing Experiment 


Since a seismometer is a velocity-sensing instrument, the particle motion induced 
by the propagating wave system is converted to an output voltage. which is proportional 
to the particle velocity, which can be used to obtain the wave velocities. The particle 
velocity data obtained from the seismometers were processed using a time versus 
frequency decomposition method. This method involves dividing the seismometer data 
records into equal length blocks with a pre-defined overlap. The midpoint of each block 


can be associated with an arrival time based on the separation between the bowling ball 
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and the seismometer, which in turn can be associated with a wave velocity. The power 
spectral density of each block is calculated and stored in a matrix, which associates 
frequency with velocity. The results of the data processing are then displayed in the form 
of a Gabor plot, named after 1ts inventor Professor Gabor of the University of Delft, the 
Netherlands. The Gabor plot depicts the decomposition of the wave field as a contour 
plot of slowness (inverse speed) as a function of frequency. From the Gabor plots, 
estimates of the propagating velocity of various wave modes can be made. 

Furthermore, since the particle displacement of the beach sediments is 
proportional to the particle velocity, the received seismometer signals can also be used to 
graphically display the average particle motion observed at the seismometer. The 
resulting particle motion can be displayed in the form of a Hankel plot or hodogram. By 
plotting the longitudinal component of motion versus the vertical component as a 
function of time, which is known as a Hankel plot, the motion can be used to determine 
the wave type, 1.e. P, S, or R. A hodogram 1s the two-dimensional form of a Hankel plot. 
Appendix G contains the MATLAB programs used to analyze the ground truthing data. 

The first set of ground truthing measurements were conducted as shown in Layout 
I of Appendix D7. Two series of bowling ball drops were made. The first series was 
conducted approximately 45 feet from the surf zone, in dry sand. The bowling ball was 
dropped at distances of 6, 12, 18, and 24 feet from the seismometer. The second series 
of drops were conducted at the edge of the surf zone using distances of 6, 12, 18, 24, 36, 
and 48 feet. 

During the analysis of the seismometer data it became evident that even a source- 
receiver separation of 48 feet did not allow for the desired time separation of the 
propagating modes in the wavetrain. Thus, an accurate decomposition of the different 
wave types could not be accomplished. Figures 5.5 and 5.6 show the time series of the 
longitudinal and vertical velocity signals received at the seismometer for distances of 24 
feet in dry sand and 48 feet in wet sand, respectively. The time series depicted have been 
adjusted so that 0 sec corresponds with the impact of the bowling ball. Figures 5.7 and 
5.8 are Hankel plots of the particle motion at the seismometer for the associated time 


series. 


50 


Recalling that Rayleigh waves induce elliptical particle motion, its clear from the 
Hankel plots that the motion shown at these distances is due to the presence of more than 
just the Rayleigh wave. Thus to adequately estimate the various wave speeds would 
require allowing greater time for wave separation to occur. This could only be 
accomplished by using larger separation distances between the seismic source and the 
seismometer. 

To promote wave separation, the same experiment was conducted a second time 
using much larger distances. Figure 5.9 and Layout II in Appendix D7 show the 
experimental layout for the second phase of ground truthing experiments. Separation 
distances of 50, 70, 90, 110, and 130 feet were selected. Also, a second seismometer was 
placed at a distance of 10 feet behind the first to allow for the collection of five additional 
data sets at intermediate distances. Furthermore, it was decided that the ground truthing 
efforts be concentrated in the drier areas of the beach since the Bass Shakers were not 


waterproof. 
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Figure 5.5. Time Series of Seismometer Velocity Signals for a 24 ft. 
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Figure 5.6. Time Series of Seismometer Velocity Signals for a 48 ft. 
Source-Recelver Separation Distance 

0.05 
a 
m 
p 
i 
; 0 
u 
d 
e 
[V} 

-0.05 

O 
0.1 
0.15 ae eh 
time [sec] 


Figure 5.7. Hankel Plot of Seismometer Velocity Signals for a 24 ft. 
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Figure 5.8. Hankel Plot of Seismometer Velocity Signals for a 48 ft. 
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Figure 5.9. Picture of Second Ground Truthing Experiment 
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Increasing the separation distance between the drop location and the first 
seismometer greatly increased the ability to decompose the wave field. It was not until 
the data files for distances of 70 feet and greater were analyzed that modal separation in 
the Gabor plots became clear. Figure 5.10 is the Gabor plot for a separation distance of 
70 ft (21.3 m). Figure 5.10 shows the existence of at least two wave modes, the slower 
mode being the Rayleigh wave as determined by analysis of the hodogram of the time 
series velocity signals. The group velocity of the Rayleigh wave at 40 Hz is 
approximately 88 m/sec (0.0114 sec/m) and the second mode is another wave with a 
group velocity of approximately 114 m/sec (0.0088 sec/m) at 60 Hz.. The identity of the 
second wave is unknown since its speed is not consistent with that expected for a shear 


wave in this sediment. 
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Figure 5.10. Gabor Plot for a 70 ft. Source-Receiver Separation Distance 


(1dB contour separation) 
The forces generated by the bowling ball impact with the beach sediments seems 


to generate larger amplitude body waves, as compared to the point source discussed in 


section II C. This is clearly evident in the time series for the seismometer longitudinal 
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section II C. This is clearly evident in the time series for the seismometer longitudinal 
velocity signal displayed in Figure 5.11. Even at a distance of 130 feet, the longitudinal 
particle motion caused by the body waves is comparable to that of the Rayleigh wave. 
The displacements in the vertical direction are almost entirely dominated by the Rayleigh 
wave. 

Figures 5.12 and 5.13 are the Hankel plot and hodogram, respectively, for the 
time series of Figure 5.11. These figures serve to illustrate how much more prevalent is 
the retrograde elliptical motion of the Rayleigh wave at these larger distances. Figure 
5.14 is the Gabor plot of the wave decomposition at a range of 130 ft (39.6 m). As can be 
seen from the figure, two propagating wave modes are present in the velocity range 60- 
250 m/sec. From the Gabor plot, the Rayleigh wave speed is estimated at approximately 
89 m/sec (0.0112 sec/m) and the second wave has a speed of approximately 109 m/sec 
(0.0092 sec/m). The ratio of the Rayleigh and shear speeds obtained from this plot are 


within 10% of the theoretical ratio predicted by the theory of elasticity. 
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Figure 5.11. Time Series of Seismometer Velocity Signals for a 130 ft. 
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Figure 5.13. Hodogram of Seismometer Velocity Signals for a 130 ft. 
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Figure 5.14. Gabor Plot for 130 ft. Source-Receiver Separation Distance 


(1 dB contour intervals) 


Figure 5.15 shows an aggregate display of the accelerometer trigger signal (lowest 


trace) and the seismometer vertical velocity versus time recorded for the bowling ball 
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Figure 5.15. Rayleigh Wave Progression as a Function of Space and Time 


oe 


drops between 50 ft (15.2 m) and 130 ft (39.6 m) range. The behavior of the velocity 
amplitude with range should not be considered indicative of what would be observed for 
an end-on array of five seismometers receiving a single signal, as site conditions probably 
changed over the duration of the data collection. The large vertical amplitude of the 
Rayleigh wave is easily distinguishable in the more distant traces. The extracted group 
velocity 1s approximately 90 m/sec (see straight line in Figure 5.15), which agrees well 


with speeds estimated from the Gabor plots. 


c. DISCRETE-MODE EXCITATION 


With the necessary ground truthing experiments completed, the next phase of the 


research was to begin testing the Mod IV field test source (See Figure 5.16). The goal of 
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Figure 5.16. Picture of Mod IV Field Test Source 
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the research conducted from this point on was to prove that by inducing the appropriate 
elliptical particle motion into the beach sediments at the source, Rayleigh waves might be 
preferentially excited. Based on the theory for surface waves in a homogeneous, 
isotropic, elastic half-space, the particle motion induced by Rayleigh waves at the 
boundary surface is retrograde elliptical, with the vertical motion approximately 1-% 


times the longitudinal motion. 


IIe Propagation Characteristics 


Before moving onto accessing the Mod IV source’s ability to selectively excite 
Rayleigh waves, some brief experimentation was conducted to see if the beach site itself 
had a preferential direction for energy propagation. To investigate this, the Mod IV field 
test source was positioned with two seismometers at equal distances, in the longitudinal 


direction, on both sides of the source, as shown in Figure 5.17 and layout I of Appendix 
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Figure 5.17. Picture of Wave Propagation Experiment 
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D8 (Also see Appendix D8 for associated test equipment setup). 

The source as shown in Figures 5.16 and 5.17 was positioned so that, facing 
seaward, the port Bass Shaker faced in southerly direction along the beach and the 
starboard Bass Shaker to the north (See Figure 5.2 for orientation). Analysis of the 
shaker motion indicated that for an equal amplitude, positive input voltage, a negative 
phase for the port shaker and a positive phase for the starboard would generate retrograde 
elliptical source motion with respect to the southern positioned geophone. Reversing the 
phase would then induce a prograde elliptical motion. 

For a homogeneous, isotropic, elastic medium the ratio of the longitudinal to 
vertical particle displacements for the Rayleigh wave at the boundary surface is a well- 
defined function of the Poisson’s ratio. But since the Poisson’s ratio for the beach 
sediments was not known, the specific ratio to be used in driving the test source was also 
unknown. As a result, it was decided that initial source testing be done using a 90 degree 
phase difference between the two shakers. This was accomplished by setting the port 
shaker to —-45 degrees and the starboard shaker to +45 degrees, or vice versa. To 
distinguish between the two cases the following convention will be used: a —45 degree 
phase for the port shaker will be denoted by -45P and a +45 degree phase by +45P. 

Since the Bass Shakers suffered the least distortion at 250 Hz, (See section 
IV.D.1, Figure 4.6), this frequency was chosen as the initial drive frequency. A 50 V 
peak, 20 cycle, 250 Hz, narrow-band filtered tone burst was used to drive the shakers. 
Even though this frequency was optimal for the shakers, it became immediately evident 
during data collection that this drive frequency was not optimal for energy propagation 
along the beach. As a result, it was decided to reduce the drive frequency to 100 Hz. It 
has been measures that, for this frequency, the harmonic distortion is less than 1% for a 
50 Vrms drive voltage. Figures 5.18 and 5.19 are the Hankel plots for the seismometer 
signals at a range of 30 ft using a 50 V peak, 20 cycle, 100 Hz, narrow-band filtered tone 
burst and a —45P phase configuration. Figure 5.18 depicts the particle motion as 
observed at the southern seismometer and Figure 5.19 for the northern seismometer. In 
both figures the elliptical particle motion is prograde with respect to the seismic source 
and clearly shows the presence of the Rayleigh wave, but the northerly propagating wave 


is much more uniform. It can also see from the figures that the amplitude of the northerly 
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propagating wave is larger by a factor of two. Appendix H contains the MATLAB 
programs used to analyze the data gathered during the beach phase experimentation with 


the Mod IV source and to generate the accompanying figures. 
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Figure 5.18. Hankel Plot of Southern Seismometer Velocity Signals for a 30 ft. 
Source-Receiver Separation Distance (-45P) 
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Figure 5.19. Hankel Plot of Northern Seismometer Velocity Signals for a 30 ft. 
Source-Receiver Separation Distance (-45P) 


Figures 5.20 and 5.21 are the Hankel plots for the same seismometer orientation 
and range, but for a +45P phase setting, which also produced an elliptical particle motion 
that was prograde with respect to the field test source. Again we see that the Hankel plot 
for the southern seismometer (Figure 5.20) is non-uniform and has a smaller amplitude as 
compared to the plot for the northern seismometer (Figure 5.21). Based on these figures, 
it would seem that the beach sediments either have a preferential direction of energy 


propagation or that the force generated by the motion of test source is different in the two 
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Figure 5.20. Hankel Plot of Southern Seismometer Velocity Signals for a 30 ft. 
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Figure 5.21. Hankel Plot of Northern Seismometer Velocity Signals for a 30 ft. 


Source-Receiver Separation Distance (+45P) 


The next set of the beach experiments was used to investigate the typical 
propagation distances that could be achieved with the discrete mode source. Since the 
results of the previous source tests were not available at the time of conducting this phase 
of the beach experimentation, the choice for follow-on seismometer position was 
completely arbitrary. 

The two seismometers were initially positioned at distances of 30 and 45 ft and 
leapfrogged out to final distances of 90 and 105 ft, as shown in Layout II of Appendix 
D8. Using the same 50 V peak, 20 cycle, 100 Hz tone burst, sufficient Rayleigh wave 
signal strength was observed out to a distance of 90 ft. which far exceeded any of our 
initial expectations. Propagation distances of this order are just one of the features that 
make seismic sonar research so extremely important to mine warfare applications, since 


the maximum propagating distances for the best ground penetrating radar and 


magnetometers 1s only 25 ft. Figure 5.22 shows an aggregate display of the source 
vertical acceleration (lowest trace) and the seismometer vertical velocity versus time at 
ranges from 30 ft (9.14 m) to 90 ft (27.43 m). Each seismometer signal trace has been 
normalized by the source vertical acceleration for that trace, as measured by the base 
mounted accelerometer. The behavior of the velocity amplitude with range should not be 
considered indicative of what would be observed for an end-on array of five 
seismometers receiving a single signal, as site conditions (e.g. tide, surface moisture, etc.) 
were evolving during the course of the data collection. The Rayleigh wave is clearly 
visible as the large amplitude wave train in each trace. The extracted group velocity is 
approximately 86 m/sec (see straight line in Figure 5.22), which agrees very well with 


that observed during ground truthing measurements. 
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Figure 5.22. Rayleigh Wave Progression as a Function of Space and Time 


Another interesting observation from this phase of the beach experimentation was 
that analysis of the particle motion observed at the seismometers showed prograde 


elliptical motion with respect to the source for all ranges, seismometer orientations, and 
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phase configurations. This result was unexpected, but can be explained. Assuming the 
beachfront to be a nearly homogeneous, isotropic, elastic medium, then the Rayleigh 
wave speed is independent of frequency (Refer to section II B). For a Rayleigh wave 
speed of 86 m/sec, a 100 Hz drive frequency corresponds to a Rayleigh wavelength of 
0.86 m. Experimentation shows that the shift from retrograde to prograde elliptical 
motion for sediments occurs at approximately 0.12, or 8.6 cm (3.4 in) [Ref. 9]. Since the 
selsmometers nominally sense the average particle motion and are nearly 25.4 cm (10 in) 
in length, most of their bulk resides below this depth. Thus, the observed motion is 
consistent with theory. 

In order to verify this observation, another set of experiments were conducted 
using a 50 Hz drive frequency. At 50 Hz, assuming the same Rayleigh wave speed of 86 
m/sec, the Rayleigh wavelength is 0.34 m. For this wavelength the corresponding turning 
point occurs at a depth of 17.2 cm (6.8 inches). To help eliminate any chance of 
misinterpreting the data the seismometers were only buried to a depth of 12.7 cm (5 in). 
Figures 5.23 and 5.24 are the Hankel plots of the seismometer velocity signals at 60 feet 
for drive frequencies of 100 Hz and 50 Hz, respectively. 

As can be seen from the figures the sense of the particle motion is clearly 
different. Analysis of the Hankel plots showed that the particle motion as shown in 
Figure 5.23 was prograde elliptical and that of Figure 5.24 was retrograde elliptical. This 
illustration serves to show that accounting for the vector nature of the Rayleigh wave is 
extremely important, not only for discrete-mode propagation, but also for selective 


reception. 
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‘a Selective Excitation 


One of the goals of this thesis was to attempt to estimate how much of the input 
mechanical energy was being converted into Rayleigh waves. But because of the 
unpredictable nature of the source motion, this was not possible. It was assumed with the 
source design chosen and the information obtained on the operating characteristics of the 
Bass Shakers that it would be relatively easy to generate elliptical source motion. During 
the beach phase experimentation it was discovered that 1t was nearly impossible just to 
create circular motion at the source. Using equal amplitude drive voltages for the shakers 
and a -45P phase configuration should have produced retrograde circular motion and a 
+45P, prograde circular motion. But, as can be seen from the hodograms of the base 
mounted accelerometer for the two cases, Figures 5.25 and 5.26, this was clearly not the 


result. 
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Figure 5.25. Hodogram of Accelerometer Signal (100 Hz, -45P) 


67 


0.05 





vertical amplitude [V} 


o 
O 
Oo 


-0.1 -0.05 0 0.05 OF 
longitudinal amplitude [V] 


Figure 5.26. Hodogram of Accelerometer Signal (100 Hz, +45P) 


Analysis of the hodograms from all the data files showed that the motion 
produced at the source depended explicitly on voltage, frequency, and _ phase. 
Furthermore, there is also a strong dependence on the elastic properties of the 
propagating medium. Examination of the accelerometer hodograms showed _ that 
longitudinal stimulation of the beach sediments was much more effective than vertical 
stimulation. 

Follow-on experiments are planned in which it will be attempted to produce more 
optimum source motion. Even though the source motion was clearly not optimum, this 
did not seem to play a major part in stimulating the desired Rayleigh wave, at least at 100 
Hz. Figure 5.27 is the Hankel plot of the seismometer velocity signals at 60 ft for a 20 V 
peak, 20 cycle, 100 Hz, tone burst with no phase difference between the shakers. Figure 
5.28 displays the portion of the time record corresponding to the Rayleigh wave and the 


ratio of the longitudinal to vertical amplitude over the same time period. 
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Figure 5.27. Hankel Plot of Seismometer Velocity Signals for a 60 ft. 
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Figure 5.28. Ratio of Longitudinal and Vertical Rayleigh Wave Components 
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The uniformity of the wave is remarkable. Similar calculations were done for the 
seismometer velocity signals at 60 feet when phase differences of —-45P, +45P, -32P, and 
+32P were used. In all four cases, similar looking Rayleigh waves were excited and the 
ratio of the wave components was nearly identical. 

From analysis of the data files gathered during the beach phase experimentation it 
was obvious that the Rayleigh wave was the most significant disturbance along the 
interface boundary. At distances of 45 ft and smaller it appeared that no matter what 
phase configuration was used, all wave types were present, with the Rayleigh wave and 


vertically polarized shear wave the most dominant. But since body waves along the 
surface decay like £ and Rayleigh waves decay like £ , the influence of the body waves 


was negligible at distances over 60 ft. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This final chapter briefly synopsizes the noteworthy observations made during the 
conduct of the research described in this thesis and the conclusions that can be drawn 
from these observations. It also includes recommendations for follow-on research. 

The major goal of this research project was to evaluate the concept of a discrete- 
mode interface (Rayleigh or Scholte) wave source. By selectively exciting the desired 
interface waves, it was theorized that significant reverberation suppression could be 
achieved, thus increasing the signal-to-noise ratio of a target echo in a seismic sonar 
detection system. 

Based on the precepts from the theory of elasticity and the known characteristics 
of the desired interface waves, two different source concepts were designed and tested. 
The first discrete-mode source, 1.e., the concept source as shown in Figure 4.1, was used 
to qualitatively probe in a laboratory environment the surface wave field created under 
various source motions. Experimentation with the concept source showed that a complex 
relationship existed between input source parameters, such as, frequency, phase, and 
driving voltage, and the propagating characteristics of the medium, and the resulting 
surface wave field. The limitations of the concept source and those imposed by the test 
tank made development of a field test source a necessity if any real progress was to be 
made. 

To design a source that could be used under varying field conditions first required 
identifying a suitable vibration source and then constructing a matching source-to-ground 
coupling mechanism. Figure 4.3 shows the first version of the field test source. Again 
the beginning phases of the field experimentation were directed at evaluating the effects 
of varying source motions on the surface wave field. Results obtained during the various 
field tests lead to several modifications of the source-to-ground coupling mechanism until 
arriving at what was thought to be the best design (See Figure 5.16). It was with this 
final source design that all the beach phase testing was conducted. 

Based on the known features of Rayleigh waves, 1t was believed that by 
stimulating the sediments at the source in an elliptical motion, Rayleigh waves could be 


selectively excited. To test this theory a series of experiments were conducted in which 
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the field test source was driven under several different modes of excitation. During these 
experiments it was observed that multiple types of interface waves and body waves were 
generated no matter what input source parameters were used. At source-to-receiver 
ranges greater than approximately 15 m (45 ft), the Rayleigh wave was clearly the largest 
interface disturbance, while the effects of the undesired waves quickly became negligible 
and were virtually undetectable at ranges greater than approximately 25 m (75 ft). 
Evaluation of the base-mounted accelerometer on the source showed that the field test 
source did not produce the elliptical motion expected for the applied input parameters. 
Because of this, it is difficult to determine whether the excitation of the undesired waves 
was due to the non-optimum nature of the source or that these undesired waves will 
always be present in practice. However, it was found that the medium itself acted as a 
selective filter for the interface waves after a few tens of meters of propagation. 

Based on the experimental results obtained there are two major areas in which 
future discrete-mode excitation research efforts should be concentrated. The first area 
deals with correcting the deficiencies noted with the Bass Shakers and the source design. 
The Bass Shakers were not designed to be operated in the fashion utilized during the 
source testing. Thus, higher quality, greater power-handling exciters (e.g. shakers, 
actuators) must be incorporated into a new source. New sources must be identified that 
can be driven harder, but experience less distortion. Furthermore, complex source 
geometries such as the 45-degree orientation of the shakers used in the field test sources 
should be avoided. Such configurations simultaneously generate both longitudinal and 
vertical seismic motion. Future source designs should utilize simple vertical and 
horizontal excitation. This would serve two purposes, the first being that such an 
orientation still allows for the generation of elliptical motion at the source, and second 
that each component of the motion can be excited individually to study the resulting 
effects. Work has already begun in designing and constructing a new source made of 
orthogonal, servo-controlled, 20 Ibf, linear actuators. 

The second area deals with re-examining the experimental setups in order to 
better sample the various seismic wave types. This would serve to increase the 
researcher's knowledge of the different excitation potentials of the possible wave types 


and their associated excitation depths. Along these lines more experiments are needed to 
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further explore source control for optimal wave generation as a function of frequency and 
excitation depth. To better accomplish this task efforts should be made to computerize 
the source control so that the source motion can be monitored on a real-time basis. 

Even though the results of this thesis were not as informative as initially desired 
the basic concept of discrete-mode excitation is promising. |§An_ important 
accomplishment of this thesis was to get the research out of the laboratory and onto the 


beach where the follow-on efforts can begin. 





APPENDIX A. THEORETICAL PROPERTIES OF RAYLEIGH WAVES IN A 
HOMOGENEOUS, ISOTROPIC, ELASTIC HALF-SPACE 


This appendix contains the various MATLAB programs that were used to 
reproduce the theoretical graphs contained within the reference texts on Rayleigh waves, 
some of which appear as figures within the body of the thesis. It also contains other 
programs written and their output, which provide further visual aid into the behavior of 


Rayleigh waves. 


‘ps 


APPENDIX Al. VARIATIONS IN WAVE SPEED RATIOS AS A FUNCTION OF 


oP WP oO o}'0 


oP 


AO WP WW oO 


oo 


o\O 


POISSON’S RATIO 


Filename: wavespeeds.m 

Written by: F.E. Gaghan 

Date Last Modified: 22 December 1997 

Purpose: This program numerically solves the Rayleigh wave equation 
in an isotropic elastic half space. The wave equation 
solutions are expressed as the ratio of the Rayleigh wave 
speed to the shear wave speed as a function of the 
Poisson's ratio. The ratio of the three wave types to the 
shear wave speed is displayed graphically for Poisson's 
ratios in the interval [0,1/2]. 


Initialize required storage vectors 

clear all; 

SPecaurat1esqa Toor l—| |7 sspeca tae tos dedmnocr mn, 
Becedurorl@sdeaGockes—| i, 

longsqd ratiol=[]; longsqd ratio2=[]; longsqd ratio3=[]; 
VeriEsqd ratiol=(]|; vertsqd ratio7=||;,¥verrsad ratio, 
Pay Foveheqd yy yoce— |r, 


For loop used to run throuch the desired FPolssom Ss gari1cs 
j=l; 
formrcnu—/O. O40 0 -, 

geanmasad (a j= 4 la 2a (eee ee 

inv_gammasqd (j)=1/gammasqd(}j); 


ovo 


Calculate the coefficients of the Rayleigh wave equation which is 
cubic in the wave speed ratio 

p=-8; 

q=24-(16*gammasqd(j)); 

r=(l16*cgammasaqd(3))—loe; 

a= (35S ae oe 2: 


fo)\e4 


Ba 2" (pS i= (3p rd) ten eet 2: 
f=sqrt(((5°2) 74) 44 (a> 327 
G=atby Zick a; 

== (7 2) i, 


fe)\9) 


Perform necessary manipulations of the coefficients to acquire 
values in the correct quadrants when taking cube roots of complex 
quantities 
Le (ih SOy ee (O50) 

BS aes) (lea) ) >> 1k; 
else 

A=g* (1/3); 
end; *enad-1f (£20) £ (e<0) 


fo) \e4 


fo)\e4 


LE VEE Ss Ojees nO} 

B=(abs (h)*(1/3))*-1; 
else 

Ban (173). 
end;%end-if (f>0) & (h<0O) 


76 


relied 


oe 


oo 


oo 


oo 


Calculate and store the three roots of the square of the wave 
speed ratio 


re)\e) 


X1=A+B; 

ee eee 2) ( ( (AB) 72) “sarees 1) 
eres] 2) (A-B)./2) “sqme G3) i. 
yl=x1-(p/3); 

Peas leis) 3); 

y3=x3- (p/3); 


speed Fratiosgqd rootl=(speed rataosqd rootl; yl); 
Spec meee sO Or? — | Speed mae Hesecmuooee yi, 
Specdsraelosqd ©Oot3—|Spced Tariceqd soot. 73); 


oe 


Calculare and store the ratio Geslongitudinal attenuation 
coefficient squared to the Rayleigh wave number squared 
longsgqdttatiol—(longsqd ratio fi canmasqa)) 7 yl), 
fong=s@qa, matto2—| longsqd ravie7, a gamlmasgaci js ye); 
lenges tat loo=|longsqd vanes, 1 (ganmasqd( 7) 72) |; 


o\0 


fe\\2y 


Calculate and store the ratio of vertical attenuation coefficient 
Squared to the Rayleigh wave number squared 

Wo ee scuedta © P| Vere sGd pear Lola ble 
VeLesdamedielo2—(Vertsqd tabive,, | 721, 

Tepes iemedetoo— | Vert ogeeral Ou,’ al, 


fe\\2 


J=3+1; 
end;%end-for nu 


Store the vectors into matrices 

Speed ratrosqa  rOors—|speced ratllosdd woor sel ed) ral eagem oor 2, 
Spcede pat losgqa soo 4), 
Hongsgatwaetos—|| longsgqd ratio! tongsqdararioOc moncgsddrraeto |), 
VELresqeprat "Os—(voursqd ratiol, ver tsqdyrdtlO7 7 vere cgamtde Le. |] 


Perform required logical comparisons to keep only the roots that 
satisfy the Rayleigh wave criteria 
Suitable=(longsqd ratios>=0) & (vertsqd ratios>=0); 
ror m0 j=—1) 
for n=1:3 
if suitable(m,n) 
Rayletghsqd root=(Rayleigqhsgqd root; speed ratilosqd roous im, 2); 
end;%end-if suitable 
end;%send-for n 
end; send-for m 


Rayleigh root=sqrt (Rayleighsqd root); 
Du | Os enOl 025 )5 
Sicateopeea=ones (1, lenge (nu)); 


Plot the ratio of Rayleigh wave speed to shear wave speed versus 
Poisson's ratio 

figure (1); 

plot (nu,shear speed,nu, Rayleigh root,nu,inv_gammasqd) ; 

aki See. 5,077); 

Grid on;-axls square; 

xlabel('Poisson’s ratio'); ylabel('wave speed/shear wave speed'); 
POSEN CLI e7  obear JG etext. ly. 05, Raylenga )5 

text (O30 5 322, “bong Led ina ley; 
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title('Variation in Wave Speed Ratios as a Function of Poisson’s 
Rae eo. 2). 


Variation in Wave Speed Ratios as a Function of Poisson's Ratio 
i 





: : 
Orr rere ire rrr rrr eri terrier ett rrr rrr re! <errrrrrr rr rrr rr rr rrr rrr rrr irr 
: : . 
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Pree er Cr rrr yy 
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‘Longitudinal 


ry e 
Cre ee erry es Siri Si iii iii 
ry 
. 


wave speed/shear wave speed 
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APPENDIX A2. VARIATION IN RAYLEIGH WAVE SPEED AS A FUNCTION 


fo) 


AX of oP oP oP OF O10 


fe )\°) 


o\? 


OF POISSON’S RATIO 


Filename: rayleighspeed.m 

Writrenapy 8.8. Gaghan 

Date Last Modified: 02 December 1997 

Purpose: This program numerically solves the Rayleigh wave equation 
in an isotropic elastic half space. The wave equation 
solutions are expressed as the ratio of the Rayleigh wave 
speed to the shear wave speed as a function of the 
Poisson's ratio. The wave speed ratio is displayed 
Graphically for Poisson Sebaelosc Meee Tintervale|-1,1/ 2]. 


Initialize required storage vectors 

clear all; 

Speed matrosqd rootl=|)|>; speed, raticsqdyrcor2=| |, 

Spec rab Losgd Trooes—| |; 

longsqd_ratiol=[]; longsqd_ ratio2=[]; longsqd_ ratio3=[]; 
vertsqd ratiol=[]; vertsqd ratio2=[]; vertsqd ratio3=[]; 
Fay rer ghoqa soo |||; 


For loop used to run through the desired Poisson's ratios 
aoe 
Poteet (hO ros 0. 2 

Gammasqd=(1=(2*nu))/(2-(2*nu) ); 


oo 


Calculate the coefficients of the Rayleigh wave equation which is 
cubic in the wave speed ratio 

p=—-o; 

gq=24-(16*gammasqqd) ; 

r=(16*gammasqd) -16; 

See ol 2) is, 


oe 


Sea (e 3) =e pea rte) sa 
f=sormr(( (bb 2) /74yet (ars) 727 
G=2(6/ 2 jai; 

hay 2) 


fe)\e 


Perform necessary manipulations of the coefficients to acquire 
values in the correct quadrants when taking cube roots of complex 
quantities 
Le eee) OS <0 | 

Petes 4e) (a) a, 
else 

A=g* (1/3); 
end;%end-if (f>0) & (g<Q) 


oo oO 


Te: Side Oy ea. <0) 
B=(abs (nh) (173).) 751; 
else 
Behe sla 
ena send—it Khoo) —& (hao) 


fie, 


oo 


Q\o 


oe 


oe 


oe 


Calculate and store the three roots of the square of the wave 
speed ratio 


role 


x1=A+B; 

x2= (= (AFB) 72 CAB 2) Ssqret—s) jy 
23= (> (ABZ) oe eee) “Sar jar, 
yl=x1-(p/3); 

V2qs2 ae 3), 

Vy a= oo oy al 


speed ratiosqd rootl=(speed ratiosqd nooEl; yy!) 
speed ratiosqd root2=([speed ratiosqd root2;y2]; 
Speed eset osqd oot s— (speed rabvocca mee ces a2 


oo 


Calculate and store the ratio of longitudinal attenuation 
coefficient squared to the Rayleigh wave number squared 
Vongsqdweeatiol=(longsqd ratecl,; v-(cammdsad sa, 

longsqd ratio2=([longsqd ratio2;1—(gammasqa- /2))); 

tongsqd ratie3—(|Vongseqd ratios, h (anmiasge. Ve), 


oe 


ao 


Calculate and store the ratio of vertical attenuation coefficient 
sguared to the Rayleigh wave number squared 

VELrESgd Tatilel=(veresqdireti Olay. 

Ventsqa raere7—|veresogdyreac1 O24, Tai 
VOLesGdaracO3—||Verrodd mae los, Haein, 


fo)\9) 


el aue 
end; senda-for nu 


Store the vectors into matrices 

Speed ratiosqd roots—|Speed ratiosqd roct |) speccyratlooqdyroorZ, 
speed ratiosqd root3]; 

longsqd retios—[longsqd ratiol, longsqd at to7, Vongsqderarta=), 
VELtsqd rabies |vertsqd ratio, vertsed rat 167, vereacemrea-os i, 


Perform required logical comparisons to keep only the roots that 
satisfy the Rayleigh wave criteria 
suitable=(longsqd ratios>=0) & (vertsqd ratios>=0); 
FOr Amalie (yan) 
fOr un la5 
if suitable(m,n) 
Rayleighsqd root=(Rayleighsgqd root;speed ratiosqd roots (m, mn) |Z 
end;%end-if suitable 
end; send-Lfor 1 
end; tsend-for m 
Rayletgh root—sqrt (Raylerchsica rooe): 
hu (O30 3054 025]; 


Plot the ratio of Rayleigh wave speed to shear wave speed versus 
POLSSOn "Ss ratio 


figure(i)+; 

POE mu, Ray loro root), 

e/iaplie ane) a er 

xlabel('Poisson’s ratio'); ylabel('Rayleigh wave speed/shear wave 
speed'); 


title('Variation in Rayleigh Wave Speed Ratio as a Function of 
Poisson s Ratio'); 
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Rayleigh wave speed/shear wave speed 


Variation in Rayleigh Wave Speed Ratio as a Function of Poisson's Ratio 
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APPENDIX A3. ELLIPTICAL PARTCILE MOTION INDUCED BY RAYLEIGH 


ao oP of AO 


oo 


oe o\0 


oe 


ov? 


oo 


oo 


WAVE 


Filename: rayleighellipses.m 

Written by: F.E. Gaghan 

Date Last Modified: 02 December 1997 

Purpose: This program is used to graphically display the varying 
particle motions induced by propagating Rayleigh waves in an 
elastic medium. Specifically it is used to show how the 
particle motion begins as retrograde elliptical and 
progresses to prograde elliptical with depth. 


Initialize required space and wave parameters 
clear all; 

k= (22 eae 1001 300; 

w=k* 300; 

t= (Oar. 7) 7/200; 

wt=w*t; 

wt _t0=0; 

C=0.2; 


Set special plot commands 
axes (box, eon); holdton; grid son; ~ax1s "equal, 


Outer for loop used to step through five various depths along the 
Z-axis within the medium 
IBopes 9.65 —3 (Oe Oh SO ee oe leash a | 


¢ Inner loop used to step through the longitudinal ranges in the 
$ direction of propagation along the x-axis 
fomesol =| OV OG ly.) 2 pee lL: 


fo)\9) 


Algebraic manipulations required to calculate the longitudinal 
particle displacement 

al=exe( 0.84 /5"k* 3) 025773" expiUrooe se kas 
el=sSim(we-(k*x1)); 

eae 1) 

=i x, 


oe 


oe 


Algebraic manipulations required to calculate the vertical 
particle displacement 

AS=-Ume4 Ss * exp (0, 641s xo 40] Sexes ke xo 
C2=cositwe— (kee 

Ws=—6 ase es); 

X3uUs=UStxs- 


oe 


6 Plot elliptical particle motion 
EMGuUre.) | 
pod Moyea cl NW ergo) BiG iy io ira vg) ee: 
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are 


ae 


Qo 


fo\\%) 


oP 


Algebraic manipulations required to calculate the wave motion 
Mecho rongi tudinal direcrion 

wee Ge 0505. 2.51; 

ese religeit xlat0)] 

Poco xs, ones (alast) > 

Steel -scam(we CO$ (kk. eto): 

ple Cey ladle cl tO), 

wit Ot ett) LO 


Algebraic manipulations required to calculate the wave motion 
in the vertical direction 

Gol eesti tO-( kis Ea 

eeu er (ao. Came): 

ee OU Sak Otis, 

Plott CO, xous TOR b> ye 


Plot Rayleigh wave motion 
quaveritx li eO (lO: best), xs0n0 (re lace) a eeOn LO. 1ae eee 
Cee TO tast), U-ieo, Ka 


end;%end-for xl 
end;%end-for x3 


xlabel('longitudinal range [m]'); ylabel('vertical depth [m]'); 
title('Elliptical Particle Motion Induced by Rayleigh Wave'); 
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APPENDIX A4. LONGITUDINAL AND VERTICAL PARTICLE 
DISPLACEMENTS 


Filename: rayleighdisplacements.m 

Written by: F.b2> Gaghan 

Date Last Modified: 02 December 1997 

Purpose: This program is used to graphically display the horizontal 
and vertical particle displacements induced by a propagating 
Rayleigh wave for a specified value of Poisson's ratio. The 
curves displayed are the particle displacement amplitudes 
normalized by the vertical particle displacement at the 
surface as a function of depth normalized by the Rayleigh 
wave number. 


Initialize required storage vectors 
eleareall; 

Mmu=O7. 25; 
gammasqd=(1-(2*nu))/(2-(2*nu) ); 


Calculate the coefficients required to solve the Rayleigh wave 
equation which is cubic in the wave speed ratio 

p=-8; 

gq=24-(16*gammasqd) ; 

r=(16*gammasqd) -16; 

Bia S le = Goma i osie 


b=( (2% (p73) =909" po) ae 
PesSqre(( (D2) (aja 3) 727 
G=— (iy 2 

h==(b/2)-f; 


Perform necessary manipulations of the coefficients to acquire values 
in the correct quadrants when taking cube roots of complex quantities 
af  (fe0) a (g<0) 
P= (alos eye eee, 
else 
A=g*(1/3); 
end; $end=12 (f£>0)) 6] (9<0) 
Lt CES Oi aren 0) 
B=(abs (h)*(1/73)) *-1L- 
else 
Boao, 
end;%end=-1f (f£>0) & (h<0Q) 


Calculate and store the three roots of the wave speed ratio squared 
x1=A+B; 


x2=(= (APE 2)4 (( (A-B) 2) “Ssqmrn— 3) 1: 
RS= (= (AFBI ZO By 72 sqrec (ai 
‘yl=xl-(p/3); 
y2=x2- (p/3); 
V3=x3-(p/ 3); 


SPeeG  rabresqd rocel—yi; 
SOCeO Sra tOSddsLCoEl— 74, 
SPeeG lr art Osqe nOcrs—yo, 
Speco vral leosqGreeous= ly, v2, voi 
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ov”? 


fo)\°) 


fe\\2) 


fo)\od 


fo)\0) 


oe 


oo 
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foe) 


Calculate and store the ratio of longitudinal attenuation coefficient 
Squared to the Rayleigh wave number squared 
tongsdqcd.7ae10l—l1-—{gammasqd”yl); 

Dongs ta sbaetoe— =k Gammasqa”~ y2).; 

Henge 9qa ratno5—l—(gammnasqd*~ ys); 

longsqd_ratios=[longsqd_ratiol, longsqd_ratio2,longsqd_ ratio3]; 


Calculate and store the ratio of vertical attenuation coefficient 
Squared to the Rayleigh wave number squared 

vertesqd) retiol iyi; 

vertsqa ratioz—l-y2; 

Vertsddarebic 3—l ayo, 

VertsSqd rabwos. (venesqd motel Venmecgemaallogmevertsqd ratios]; 


Perform required logical comparisons to keep only the root that 
Satisfies the Rayleigh wave criteria 
suitable=(longsqd ratios>=0) & (vertsqd ratios>=0) ; 
for m=1.3 
if suitable (m) 
heave igmeqd 6Ooe—Specd raliosdd  re6ors (mM); 
index=m; 
end;%end-if suitable 
end; senad—-Lor mM 


Calculate the values of the Rayleigh wave root and the attenuation 
ratios 

Ra yleorgn reoe—Sqre( Ray lei ghsqd Toce), 

Rayleigh long ratio=sqrt(longsqda ratios (index) ); 

Kay lei@n, VeremrattO—sqrt( ver tsqdmrarlos( index) ) 7 


Calculate the normalized horizontal and vertical particle 
displacements 

q=Rayleigh long ratio; 
Ss=Rayleigh vert ratio; 
IGMGRGenStaizeg- so) (se 2) i): 
Vieltleenst—2/((s 2) +1); 
depen— 70: 0-01 21.6); 

Surtagce  lavyer—2eves (1) lengen (depen) ); 

KY Z=deeehn~ (2* pi); 

US—(€.6 4a" kez (vere (Conse-expiac ies) ) )y, 
u30—(l="verG Const }; 

u3_u30=u3/u30; 

Wed So ee Ong eons exp Ss hrZ) 5 
ul USO uly (=q-us0i) 7 


Plot the normalized horizontal and vertical particle displacements 
fener) y 

plot (depth, surface layer,depth,u3 u30,depth,ul u30); 
xlabel('depth/Rayleigh wave number') ; 

ylabel('normalized displacement amplitude'); 

Ces. 22 000 . londieudinal’ )>stext(02.62,0.5;. vertical”); 
title('Longitudinal and Vertical Particle Displacements'); 

Ge raon: 
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normalized displacement amplitude 


Longitudinal and Vertical Particle Displacements 


? longitudinal 


0.2 0.4 0.6 0.8 1 LZ 
depth/Rayleigh wavelength 
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APPENDIX A5. TENSILE AND SHEAR STRESSES 


Filename: rayleighstresses.m 

Written by: F.E. Gaghan 

Date Last Modified: 22 December 1997 

Purpose: ThiS program is used to graphically display the horizontal 
and vertical particle stresses induced by a propagating 
Rayleigh wave for a specified value of Poisson's ratio. The 
curves displayed are the particle stress amplitudes 
normalized by the horizontal particle stress at the surface 
as a function of depth normalized by the Rayleigh wave 
meu eis:, 


Initialize required storage vectors 
clear all; 

Mu Ol 5; 
gammasqd=(1-(2*nu))/(2-(2*nu)); 
gammasqdinv=1/gammasqd; 
lame=(2* nu) (l= (2 nw) ); 


Calculate the coefficients required to solve the Rayleigh wave 
equation which is cubic in the wave speed ratio 

Pao) 

q=24-(16*gammasqd) ; 

r=(16*gammasqd) -16; 

eee (2) et) Ss es er 


Dante (pose CG) ater 2a: 
fosquni(ii (se 27a (la Ssyy2)\)> 
G=— (5/7 23a; 

hae 2) ce 


Perform necessary manipulations of the coefficients to acquire values 
in the correct quadrants when taking cube roots of complex quantities 
i Vee) ee aa!) 
Beams Gg) (ay ojo 1s 
else 
A=g% (1/3); 
end;%end-if (f>0) & (g<0O) 
ie Ae 0) ee ae) 
Bolas hel eel, 
else 
B= (17 3)7 
end;%end-if (f>0) & (h<Q) 


Calculate and store the three roots of the wave speed ratio squared 
x1l=A+B; 


ee (ey 2) eB) 2 Sgr tts) ); 
So (= (ht By AS ARK By yy 7) ~seqri(—3))> 
Vlei Sie 


Va 7S) 

Vo—o- (p73) 
Speedin oc roOsdder erly i, 

speed ratiosqd root2Z=y2; 

Speed Kat losqeix oor 3—y 3; 

Speed ratiosqd roots=(yl,y2,y31¢ 
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Calculate and store the ratio of longitudinal attenuation coctiticiten: 
squared to the Rayleigh wave number squared 

Vengsdd Laetoli gatas qn. te 

longsqd ratioz—l-(gamnasgqa v2), 

Voungsqa Yaric5—l— (Ganiaisqa yoy), 

longsqd_ratios=[longsqd_ratiol, longsqd_ratio2,longsqd_ ratio3]; 


Calculate and store the ratio of vertical attenuation coefficient 
Squared to the Rayleigh wave number squared 

VErLSqdpeactol lav; 

Verlsqdeearite 7 yz, 

VEresqd evar le3—i— 73; 

Vertesqd yr abico—|Vvertsqd 2obilel, (ere sdcmedete7 ie meccam i cielool—, 


Perform required logical comparisons to keep only the root that 
satisfies the Rayleigh wave criteria 
suitable=(longsqd ratios>=0) & (vertsqd ratios>=0) ; 
Porem=— ile 3 
if suitable (m) 
Pay lergnsqd root —Sspeecdirdt tes 4qcmu peer ayn 
index=m; 
end;%end-if suitable 
end; %end-for m 


Calculate the values of the Rayleigh wave root and the attenuation 
ratios 

hewe Oh rOGt—sqnr (hay lerghsqcemraccr)., 

Rayteigh long ratioe—sqrt(longsqdiracte st imdesany 

Ray lergme vere haktoO-sort(veresogd Edutocmmicde <) 


Calculate the normalized horizontal and vertical particle stresses 
g—-Rayleigh long ratio; 

s=Rayleigh vert ratio; 

long = cense—tarq os) Glee) 1) 

Vert peoust—27 (5-2) 41); 

Cepen—(02 07011. .6); 

SUtBaCe) layver—Zer Osi |, Tenguneaepcaia:, 

krz=depen* (27 pi): 

Sigmaxx=(gammasqdinv*exp (-q*krz))- 


(lame *(q° 2) "exp(>G" krz) \=(2- tendmecnist “exes kez) 
sigmaxxO=gammasqdinv=(lame* (q°2))-(2*long comsc) ; 
sigmaxx xx0=sigmaxx/sigmaxx0; 
Sigmazz=(-gammasqdinv* (q*2) *exp (-q*krz))+ 


(Vame*exe(—G" kr2))+(2* lonGmeomst = po kee 
sigmazz xx0=sigmazz/sigmaxx0; 
SPOMAXZ=(=2" 4 exp (-O8 hem) Joc Oo exept —sS kez 
sigmaxz xx0=sigmaxz/sigmaxx0; 
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ce) 


eo Opeene noOsamalized horizontal and vertical particle stresses 
Pacuise (je 

DOE eee cert ace liver, dépiijavomaxx *xx0, depth, sigmazz xx, 
depth, sigmaxz xXx0); : 
xlabel('depth/Rayleigh wave number'); 

ylabel('normalized stress amplitude'); 


omen OnO omen oye LONG ILUCInd )” \o.seewe (Umea 03>. Shear’ \> 
bextiG.42 90.16, vertical’): 
text (Os0G, Us 4 stress); text 0263-0. stress): 


texb(O0742 702 stress): 
title('Tensile and Shear Stresses'); 
Grid» ong: 


Tensile and Shear Stresses 


: shear 
: stress 


normalized stress amplitude 


: stress 





0 0.2 0.4 0.6 0.8 1 12 qe 1.6 
depth/Rayleigh wavelength 
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APPENDIX A6. VARIATION IN RATIO OF LONGITUDINAL TO VERTICAL 
PARTICLE DISPLACEMENT AS A FUNCTION OF POISSON’S RATIO 


Filename: rayleighaxesratio.m 

Written by: F.E. Gaghan 

Date Last Modified: 10 February 1998 

Purpose: This program is used to graphically display the ratio of the 


longitudinal to vertical particles displacements at the 
surface during Rayleigh wave propagation as a function of 
the Poisson's ratio. The calculations are for a homogeneous, 
isotropic, elastic half-space. 


Establish desired Poisson's ratio values 
Clear all; 
Mo Ons On Od Oo fF 
L@rer— length (nu) 
gammasqa— (l= (47 nu (1 3) 02a 


fo)\ed 


oo 


oo oO 


fo) 


oo 


oe 


Calculate the coefficients required to solve the Rayleigh wave 
equation which is cubic in the wave speed ratio 
p=-8; 
q=24-(16*gammasadqd) ; 
r—(16-Gammasqd)—16; 
a ee, 
2s eC eee 
P—sane( 2)7 4) (esi) 2 ae 
G=—(b/2)41; 
h=—(b/ 2) er; 
Perform necessary manipulations of the coefficients to acquire 
values in the correct quadrants when taking cube roots of complex 
quantities 
tee (te Oa <0.) 
A= (abs (g)* (1/3) )*-1; 
else 
Pag 
end;%end-if (f£>0) & (g<0Q) 
ek ile Ome 0) 
Ba\aesti: (ye) ) 1, 
else 
B=h* (1/3); 
ene, senad-i1f (&>0) 6 (h<0) 


Calculate and store the three roots of the wave speed ratio 


squared 

X1=A+B; 

X2=(— (AB) 2) CAB 2) sGgriei— 2s 
Koa > CAE? 2) > CASE 2 sore (=e), 
yl=x1-(p/3); 

V¥2=x2— (0725 )i; 

y3=x3-(p/3); 


speed ratiosqd rootl=yl; 
speed ratiosqd root2=y2; 
speed ratiosgqd root3=y3; 
speed ratiosqd roots=(yi, y2,y3]; 
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Saleuiauc and Stoxe the ratio of longitudinal attenuation 
coefficient squared to the Rayleigh wave number squared 
longsqd ratiol=1- (gammasqd*yl) ; 
TongedemcaelO2—1—(gammasqd* v2); 

Wong aeleOl= |S \ gammasqd y 3); 
Longsgdtatvos—|longsqd ratiol, longsqd ratio2, longsqd ratio3]; 


oe 


fold 


Calecitarewondastore Lhe ratio Of vertical attenuation coefficient 
squared to the Rayleigh wave number squared 

vertsqd ratiol=l1-yl; 

Vertsed ratvo2—l-y2; 

verts@d ratioj-l=y3; 

VEYESRC culos = iver sqd Faltol, ete GaetO2-VeLrEsqq. ratio); 


oe 


oo 


Perform required logical comparisons to keep only the root that 
satisfies the Rayleigh wave criteria 
suitable=(longsqd ratios>=0) & (vertsqd ratios>=0) ; 
for m=1:3 
if suitable (m) 
Boavleltgisqd 2reot—speed rabtesqameosts (im), 
index=m; 
end;%end-if suitable 
end; tend-for m 


fe)\%) 


fe\\e) 


Calculate the values of the Rayleigh wave root and the 
attenuation ratios 

Ray Pewon poem —Sqru\ hay lelgnsqdl cece ji, 

Ray eagle eleng pearle—sqrk (Pongsqd rakes tende «iyi, 
Rayleigh vert ratilo=sqri (vertsqd ratios {index) }; 


fe\\e) 


fee) 


Calculate the normalized horizontal and vertical particle 
displacements and the ratio of horizontal to vertical particle 
displacements 

q=Rayleigh long ratio; 

Sahay lelgh vere Paco; 

lengmeenst—(2-q°s)/{(s°2)41)4 

Vert. const=2/((s*2)+1); 

U30 vere const, 

Wit —l—long const; 

axesratio(1i)=ul0/ (-q*u30); 

end; %end-for i 


fe\\e) 


fe) \o4 


Plot the ratio of horizontal to vertical particle displacements as a 
FUNCELeGH Of etne Poisson's ratio 

ore lb ware ade ee 

plotinu,;axesratic) ; 

xlabel('Poisson’s ratio'); 

ylabel('longitudinal displacement/vertical displacement"); 

title (‘Longitudinal and Vertical Particle Displacements'); 

axis ((nu(1),nu(i),min(axesratio),max(axesratio)]); 

Gi ba on, 


7 


longitudinal displacement/vertical displacement 


0.75 


Oy 


0.65 


Variation in Ratio of Longitudinal to Vertical Particle Displacements 


Re nt Ee) ahs ce So, ne Ee ee NC ee ice ee te 2g, Ue aa ee Ie es eg ee eg es 9 SO ee ee 
CO Cee Coe OH HS ONES Se mm occenca=sceeeeeeen er rr ee rrr rit? Siiriiriritiriiirrrrrrrrrrirrririrririririrrirrriiririirr Titi tt) 
ry . ° . 


Rane ee Se eee me MOE NR Re EOEN 7e hm ae ETS AL DR eS AM ter $1 pee Re RRC RCO et Tre ee des Oa ute T cm tr at em ee Dale a a ei) AS tee ak a ilie eh Peet 
Perr ee eri ey 0000 00 0000 00 00 2am oo O08 08 000 00085 280008 OF OO OOF OED OO SOROS SES ED HEESEOEEO ORD CS CSCO OO SEES OO SOSH SEES ROCESS OOES SERS COOEEEEeS 
CT aa ee a aa aa aa ST RN cea ETE ET hatha aha adnate A Sy 


Oe ee i er ey ih WT ree rrr irri iii ri iii iri iri iy) 
* * * = 
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APPENDIX A7. VARIATION IN RAYLEIGH WAVES AS A FUNCTION OF 
DEPTH 


Filename: rayleighwaves3d.m 

Written by: F.E. Gaghan 

Date Last Modified: 22 December 1997 

Purpose: This program produces a 3D graphical representation of 
Rayleigh waves in an isotropic elastic half space from which 
can be seen the exponential decay with depth into the 
medium. 


Initialize required space and wave parameters 
clear all; 

xi= (002221 

k=(2* 01 7100) 300. 

Keak 

w=k* 300¢ 

t=(0:0.22ey; 

wt=w*t; 

C=0.25; 

[Ac =meshgrid(xly; 


Loop through various depths in the medium to calculate the horizontal 
and vertical particle displacements. 
TOF 
if Oe She Oe oy = Oe HO. oe 
$ Calculate the horizontal particle displacement 
et=exp (060475 k* x3) =(0,. 57 73 exe 2393 ae 8 
cl=sin(wt-kx1); 
Wd=€> (al re, 
Dla | 
$ Calculate the vertical particle displacement 
ajs(—-0. 64/5 exo (064 /57k* x3) 4G So texp 023935" k xsi 
GS=C05 (Wer-Kxl 
UNo=Ce tages). 
$ Conditional statement used to display surface layer 


if i== 
[LX, LY] =meshgrid(xl,u3) ; 
[LZ] =0*LX; 


Prgure (Lie 
Suis (ian pia eae 
end;%end-iff 1 


o\o 


Plot the waves for each of the four different depths within the 
medium 

[D270 Liamesigr rd (iS); 

heldon ; 

mesh( Gl .<Ud 03)? 

xlabel('longitudinal axis [m]'); ylabel('transverse axis [m]'‘); 
zlabel('depth axis [m]'); 

title('Variation of Rayleigh Waves as a Function of Depth (3D)'); 
view(20,20); 

ee L 

end;%end-for x3 


o\o 


Variation of Rayleigh Waves as a Function of Depth GD) 





[uw] sixe uidap 


transverse axis [m] 


longitudinal axis [m] 


gh Waves as a Function of Depth (2D) 


Variation of Rayle 


£18. 6)6! 0) vie Ol weaver ele: ik) 161/68) eo) wlele! 0 4p. (6) (0 (S018) 61161 m (e186: 10 36:16 q 60 060) (8:10) ses! 0 (816) 69 (G9 O91 Se ene SS eKee ©. = Ge ee 


° . 
coe eww es eee eee eee e eee 
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longitudinal axis [m] 
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APPENDIX B. MANUFACTURER SPECIFICATION SHEETS 


The appendix contains the specification sheets from the manufacturer catalogs for 


the various equipment components used during the experimental testing phases. 
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APPENDIX BI. MANUFACTURER SPECIFICATIONS FOR ELECTRO- 
MAGNETIC SHAKERS (MINI-SHAKERS) 


WILCOXON 
RESEARCH 





The Mode: F52 Vibrauor Generaicr is a rcactice 
woo snaxergeserating dynamic forces for ele tie 
excitation in v:bration research anc testing The 
reaction princisie of cocration, hott v. eight and 
compact coniguration allow this sene-ate: to be 
stud mountes in any sosition, arroct!y 1 Siruchires 
WENO extena support of critical sratt aligament 
proble rs. 


The Modal FSB Electromagnetic Snake: is a cyt 
cneat permarert magnet shaser Tre magnet is in 
rgic contact with the outer case, A mowng co.) 
wound on ar alum rum bobb:n surrounds the mag- 
re “his con and pobbia ws suspendes from two fat. 
c.roular retail saregs ensurnng pute axia mohan A 
low cecterct cravity TITUZES (OLALON Oxsilalies Sy 
the sraser, The added dynam.c weight :¢ low sinee 
fre suscerded weigh coes not effec: ictationai a 
axat nena adove is resonarce. “hus orevents 
neors:stencecs Offen erceuntored wrerever tre 
rmationa. :‘mpedarce of struciuses s low comonres 
to the axia imsedance. 


The Mccel FSB is desiqned {fcr operat es WEry 
w SO range of aucic frequencies Itcan be sugshes 
with sensing vansducer conta: mag an ind 70m 
cteranda loree gage iMlogel 211 Ir pedance Head}. 
“ne sna<er drives ihe testod Structuré through Pre 
mpadance mead 


The Madel 211 Impedartce eas s a cyundrica’ 
structure costaiing 4 mezoolecic accelero:neter 
aA a D:evoPlectnic ioree gage The transduce: basc 
caribe usec to measuro apo od force and structure 
mcton From these measurements mechanica: im- 
pedarse cin Se obtained 


The high impecance charge signals fro-r the pie- 
rofieclnc ‘ore Sage and acce eromeie’ snculd bc 
corditanad using crarse ampiuiiersa or Charan cor: 
vere (OO751), and rah mpedarce vaitage an. 
pifiars FARAS) 


Tee Moce: 211 Impedance hirad has 4 scean en 
contac! diarneier of 6 25 incres preventing Oxces 
s:¥0 5: ffomes cy impedance head asachmeant Ibe 
“ery ow mass below the force gaye (1 grams makes 
 peoss'biec to take Measuroments on relatively igh: 
S7uGtuIes. wet a> ulsamnas, Modes ans |ghi 
MaACwINe’Y. 


Model F5B/Z11 


Kiectromagnetic Shaker Svstem 


Usarie Frequency Ranae .... 
Bloexed Force Outay 

Maximum Corunuous Current .. 
Nominal Electrical Imcedance 
OS Evcetrnen! Resistarce 
Resonance Frequercy d'oc<ed 
Cernecter Wi. 4 ©... 


7 eo 18 


SPECIFICATIONS OF MODEL 211 / IMPEDANCE HEAD 


Accelerometer Nominal Values 


Cnarge Sensitivity ww ee re 
Voltage Sensitivey’ —....... A. ce 

CABAGWANCE sine odes eee) Megs eee 
F*agency Range. OS. dB cc... een, ee tices 
CoGiGaccceicn. sac eeeee 


wen ee ems 


er eeee 


Force Gage Nominal Values 


Charge Sensitwity ow. 

Voltage Sensitray’ .. 

Canactance: FeO sen ct 
Frequency Range, +3.d8...... 0. 

Connecter oo... 

GoeaRth ees se, haddian seas 


Mass Belov: Farce Gage 
fforge gage vtached ‘co speciren}. 
Wass Belov: Force Gage 
lacceleromete’ attached te secciuren: 
Ettecive Suffress 
(force gage attached to soec.men} ....... 
Et‘ective Stiffness 
(accelerometer attached te s 
Ciameter of Mounting Surface 
Younting Stud, imeg’a! 
ee Ranse.. aa 
Maximum Screw Down T TOPDUC aes 
Baca Matenal . . 
We OlOP Se iivvcsessacsiace 
VA Rue" [7 | a7 aR ere 
VALAUNV EGR we ss ice 2 asters , 


Deayony ... 


mechanical impedance 
Refers 10 connector at the ered of cabie 


py 





SPECIFICATIONS OF MODEL F5B ELECTROMAGNETIC VIBRATION GENERATOR 


10+ 10.000 He 
S0¢ grank 

o.1 amp ms 
115 42 


3.5 9C'q 

14 MIG 

&0o pF 

°C - 10,000 He 
Mcradet 10- Se 
434 12-6 


3D aa a a) 

440 MVD 

BOG 

10+ $0.000 Hz 
Microdoi 13 - 32 
Ff43-J2-6 


1 qram 


10 grams 


1 
10 qrams 
179 grams 





NOTES: * Blocked force outou! refers te the force autpu! aga rs? @ mass of infinite 


Re'e’s to corditon atthe one of the & foo, fow-roise cable (120 9*) 


ACCESSORIES SUPPLIED: Output cab'es 
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ACCESSORIES AVAILABLE: Armpi fens sgna’ cocaivoness 


APPENDIX B2. MANUFACTURER SPECIFICATIONS FOR TRI-AXIAL 
ACCELEROMETER 


Dimens. ons shown on waches excent where sotod 








3) 10-32 Thd Scker 
ennne!) Conn. Adaptors 





398801 & 339B10 354817 S54B33 356A27 






Se ee ee 

Frequency Range (%5%) [te | te 2000 | ream | res000 | iezonh | "Suz 000 
10%) | at 4000 | 0705000] 0707000 | 07400 | 2500 

Resonant Frequs [win [2m | 20 | 75 | 22 |e | 

ema y SE 
/“Resowtion (broadband) | gpx_| 0.008 | 001 [001 | 0.003 | 00008 
[Mechanical Soc Limis | zgpx | eaoo | so [foo | 3000 [sooo 
[Temperature Range *4d| SF _| 0510-260 | 0500 1250 | ts +250 | S500) | swe 
7 i > a De ee es 
Amplitude Linearity sit ttt | att SS Se 
“Base Strain Sensitviy | pe | <0000e | s0.0003 | soos | som | cocoon 
“Excitation Vollage | voc | 20030 | 1890 | 18090 | 20ws0 | 20030 
[Gonstant Current Supply [mA | 220 | 220 | 220 | 2020 | 2020 
[Output impedance ‘| ohm | «100 [| s100 | 100 | 10 | <250 | 
[Output Bas Volage —=S~~=~‘iCCOC YC dCi d|~SCi SC 
[Discharge Time Constant | second] 20s [205 | 20s | 0s | >02~~*d 
[Ground isolation (Z) | ohm | omtonal | _aplonal__| ro -| ~~ eto8 | 0 
[Element Configuration | sirvcture | shear [shear | shear | ahwer | shows 
[Sealing SOS tye] weldectrenmete | weldedihermetic | welsnerm | weigherm | woidedtermeti | 
[~Mounnng Thread ——s—~=~‘~‘stCSCte dS | toe |S «dose | vere 
[Housing | _materiat_| Wann | tuanium | ttanum | tianan | stainiees eroet 
a 2S I 
SOPRIED ACUESSOMEGAT Apap ude bo libra pital ied Peete 
ac i mm Wl cr 
[Mounting Stud (5) |] eete0s_ | ee180s | 981605 | ostase | ostase 
[Mounting Base (6) || ceoara canara | oan [aa 
soa Sao =< nL ROLY WSEREEA ERRSLUEEEN PUPAE SAU S EEASECAO, PESOLSEES OE 
 iagraie Mewbng BaP | or [par aa |e | 
Sa HS SERRE REDSIE FETOTE, Teh anti egeg sees 
Mating Cable Connectors | toe | Y_|v _ | Na | aAnKw | wn 
[Recommended Cable | series | 010] 10 | Na] om _——*+d Sts 


NOTES: 
4. NIST traceadte caldraton certificate from 10 Hz to +5% Irequaency port 


Models 354817 and 339811 obtaw: ground isolation when used with suppled —§._ Mode! 061B05 mounting stud is 10-32 to 10-32. Mode! 061A27 mounting stud 















ibs 

















~~ 


























A> — 


mounting base. Models 356A27 and 354833 aro case Kolated.  ihweaded 5-40 to 5-40. Model 081A52 is a 0-80 threaded cap screw Mode! 
3 See cables‘accassones section beginning on page 78 for a complete Oa1AS9 ie 0 4-40 Hreseoed cap screw. Models 081AS+4 and OO1ASS ¢Fe 
Cascription of the supplied and optonal accessones. Twreaded cap screws. 


6. Models 080A33 and 080461 are electrically insulated mounting washer pads. 
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APPENDIX B1. MANUFACTURER SPECIFICATIONS FOR ELECTRO- 


MAGNETIC SHAKERS (MINI-SHAKERS) 


WILCDOXON 
RESEARCH 





The Mode: F582 Vibrauor Generatcr is a reactic 
wre snaker gencrating dynamic forces for structural 
exeitalion in vibration research anc iestine The 
rHaAcion prncate of operation, ight weigrt and 
Te contiqurahion allow this gonerarc’ to be 

siud mountes in any Sesilion, aroectly to structures 
WiNOL externa: support or Critical salt airgrment 
proble rs. 


The Model F538 Electromagnetic Snake: :s a cy.17- 
sricat pesmarert magnet shaxer. Tre magnet is in 
rgic cortact with the outer case. A mowmng col 
wound on ar alum num bohb:n surrounds ine mag: 
re “his cor and pobbia « suspendes from two fat, 
c.rcular retal surngs ensurng pure axiai mown A 
law certer ot Cravity MItUZOS (OlALON axcitation by 
the shaser. Tne added dynam c wcisht s lov since 
ihe suscerded weight does not eac: tetationai o 
dual nena above ts resonarce. “hs porevonts 
neorsiscenecs offen erecuntored wrerever tre 
fctianona. mpedarce of sireciures s low comaures 
te the axia impedance. 


The Mccel FSB is desanad {or operaton over 3 very 
w de range Gfaucie feaquencies. Itcezn be supslies 
wath a sensing VANGAUCEL COMANING AN ACnMeroM- 
cleranda force gage (Medel Z11 Ilr pedanee Head}. 

The shas<er drives iP) wasted Sstrumiuré Ihrounh fe 
mpadance “ead 


vhe Madat 411 Impacartce reac § a cy indrica’ 
Structure cottaining a piezoodlediia anceleromeler 
aadap-e7oelectriclorce gage The tracsduces base 
canbe uses lomecasuro aspked fe:ce and structure 
metion. From these measuremerts mochanica: im- 
padacce can ce obfamed 


Troe high meoesance cnatge signals tron the pre- 
ZOCPCING ‘orce gage and Ecce eromete’ srculd be 
cordit-annd USING chess y ampiiers or charon cen: 
vere (CO757}, and tratiompedarce voilage an. 
pl.igss (ALAS) 


Tee Mode! Z11 In.pedance head has a scecn en 
Ssontant Guavrneter of O75 ncres preventing oxces: 
sero st ffenrs cy mpedance haad attachment The 
vee Ow rinass below the force gage (1 grams Takes 
i PCss bie to Take MeassromMents Or relatively Ligh: 
SUGiules. SLE ay arf-ames. Modes ara light 
THAC™ ING. 


NOTES: 


Usanie Frequency a 
Bloexad -orre Oute 

Maximum ie ous Current .. 
Nominal Electrical tmoedance 
CC Elcerneal Rasistarce ...... 
Resonance Frequercy bocx<cd 
Connector 299 ww. 


Accelerometer Nominal Vatues 


Cnarge Sensitmity 

Voltage Sensitivity’... eae 
GAGRC TAMER” 5. Jcse0::.(. aera, erases) esau 
Froquency Range, -3d3 .. 

Canrector 

Gagic ...... 


Force Gage Nominal Values 


Charge Sensitwty oo... 

Vohage Sensitrwey . >. 
CRDACHANGE- oo nus ee 
Frequency ae 308. 
Connecter’ . 

Cahie ...... 


Mass Belov: Farce Gage 
fforge gage attached %¢ specinen) . 
Wass Below Fore Gage 
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SPECIFICATIONS OF MODEL 211 / PAPEDANCE HEAD 
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Model F3B/Z1) 


electromagnetic Shaker System 





SPECIFICATIONS OF MODEL F58 ELECTROMAGNETIC VIBRATION GENERATOR 
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APPENDIX B2. MANUFACTURER SPECIFICATIONS FOR TRI-AXIAL 
ACCELEROMETER 


Dimensans shown en esches exceot where AoKx 





Te 
} 37 Cube 
ky 
y?P 
| a 
ft, 
10-32 Thd Mig Hole 
“apecteanem 
— 
: > f ' 
Coos an 
10 t integras 6 Conductor 
Cable with 





3) 10-32 Tnd Somer 
ernmunal Conn. Adaptors 





338801 & 339B10 354817 354B33 

















GENERAL. PURPOSE 


sei: 
rent an A EE 
10 e700 4000 [0.705000 | 0707000 | 07104000 | 25000 
| “Hosonant Frequency sd Cie tems Cte Te 
SS a eT a a es Cl 
[Hesouion (broadband) pe | oon oot 9s Ps 0008 
[Mechanical Shock Limis [gpk {eam [000 | soon | s000 | e000 
Temperature Range st ~*d a 00 880 | 45004250 | 05004080 | 250-250 | esto «750_| 
“Temperature Coettinant” «dy a] eoes cous | ooo | <00s | <008 
Se Se SS TS ae ay TE eS 
"Transverse Sensitivity | % [ss | <5 | ss—S—=«dY(C ‘Sd CSS 
“Base Strain Sensitviy | ge | s0.00s | s 0000s | some | som | coona 
“Excitation Voltage _——S*d SVC | 20090 | 18090 | 1680 | tow | 20100 | 
“Constant Current Supply | mA | 2to20 | 2020 | 2020 | 220 | 2w20 
Output impedance SS OS 
































orcas 

a a 
[Element Configuration | structure | sheer | sheer] sheer] sheer | anew | 
[Sealing SS tye | wldecthermete | waldedinermetic | weldherm | woldherm | waldedtrermeti | 
[—wewnang Trea ae [rose [tose [sao | roseenveet [vee 
a 


9 a ee) sess | se | een | seas | 
LECH) ROR EEED PUPCRTRAL SDE EN OSA TERE SOCES PLOR ESET, 


fanosive Potro Wax CO 
Mounting Stu (6 es a 
iounting Base (6) Se | 
feiss Sees < ROSE) SSSR BES CIREEDELE, a 
Magnove Mounting Base ee a 
(2 HC SERRE LOSES 2 EOE SEELECIOT WEEORERES CPE TS TERE 
Sead eecomeaee ete ee eee cae 
Recommended Cable [series [ot Aten) 


NOTES: , 
4. NIST traceadts calibrabon Cerlificale from 10 Hz to +5% frequancy poirt 


2 Models 354817 arxd 339811 obtain ground isolation when used with supplied 6. Model 061805 mounting stud is 10-32 to 10-32. Model 061A27 fmounting stud 











mounting base. Models 356A27 and 354833 are case rotated.  inveaded 5-40 to $40. Model 081A52 is 4 0-80 tnreacted cap screw. Mode! 
3 See cables/accessones section beginnang on page 78 for 3 complete eee ga ma threaded cap screw. Models 081A54 and OB1A55 ars 
Cescription of the supplied and optonal acoessonas. throaded cap Scrows. 


6. Model 086433 ard OB0A61 are olactricalty nsulated mounting washer pads. 
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APPENDIX B3. MANUFACTURER SPECIFICATIONS FOR BASS SHAKERS 












“4 wer 
4% . ha 4 
Ay Bec are <“ 
Lm . , dg 
oe nT, 4 a * 
e™* WN ee ge nt at eo i 
age . oh Ait; ¥ ALS fae ne ie ¥ 
an", 4 ee * ka Sips te eX 


Bass Shakers” Let You Ay ; Mir. Sugg. Retail $199. é 
Feel All the Bass without . Ga oO: 


Breaking the Sound Barrier 
The 8ass Shaker 's a transducer that generates the 
sensation cf sound by vibration, no: by moving air. The result is a biq bass effect without 
a high pressure ievel t1at could distort sound or alow your speakers. They are ziso great 
for adding bass tc vehicles with coor sourd insulation, such as trucks. Mourt these Bass 
Shakers under your front seats ard you will feel the ounch of the drums anc the kick of the 
percussio’, as if you were Gn Stage or in tha stucio wren the music was recorcec! Hooks 
iy eG io any Standard audio ampit‘ier output channel. and can be used with existing 
jowoolers for an extra kicx. Each Bass Shaker requires only ZoW RMS of power. 
M ode! AST1B4. Item No. B-40070-492383 S/H $9.95-4°"° h/t 2 


AO He NSomam 





._— 





Specifications: AST-1B-4 Bass Shaker 
Fraime S1Ze: ..... eee. 4.75” dia, X 2.35” h 
Magnet [ypen ene Ceramic 

POWEr Kaui: 22a eres Nominal 25 W 

Force, Nominal: 3-2) 10 lbf (44.5 N) 
WEIN: ois. nee 2,5 Ibs cach 


Resonance Frequency (fo): ... 42 Hz 
Usable Frequency Range: ...... 20-100 Hz 


ye) 


Configuring Shakers. 


The Aura Bass Shaker™ has a nominal impedance of 4 ohms. 
Depending on the amplifier dnving them, they can be hooked up 
in stereo or bridged/mono configuration, just as regular 
Subwoofers. Please consult the amplifier manual for connection 
details. 


Connecting Multiple Shakers 


Two Shakers can be wired together to increase the delivered 
power . This can be accomplished by connecting them in either 
series or in parallel. Connecting in series increases the impedance 
to 8 ohms (see figure A), whereas, a parallel connection (figure 
B) reduces the impedance down to 2 ohms. Consult the amplifier 
manual to figure out the delivered power at various impedances. 
In either case do not exceed the rated 25 Watts RMS per 
Shaker (50Watts RMS, total). 





Page 3 
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APPENDIX B4. MANUFACTURER SPECIFICATIONS FOR UNI-AXIAL 
ACCELEROMETER 














"et TB. BERS 5 sf nts 
Vollage Sensitivity (1) | mvig | 
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Transverse Sensiiv 
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: 
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TAgheswe Mount | 
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"Water Resistant Connection |__| 
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Mounting Stud (9) _——s 
Adhesive Mounting Base eae 


Magnetic Mounting Base poy Ge) 


re as Soe ees eens 




























NOTES: 


4. 


« 
wo: 


1. Sensitivity tolerance is t 5%; bohter tolerances (¢ 26) avauiable upon request. 
a: 
3. Models 353813, 353817 and 353867 suppiied with 10 ft. ndbon cable 


Transverse sensitivity on shear structured acogierometers @ typically <3. 


termnating in a 10-32 coaxial piug. Cable is user repairabie. 
See page B for a desonption of tne optional modnis. (Speahcabons may 
d.tler stignity Consut factory before ordenng. ) 


Optionat high sensfivity Modets (353865, 353866, 353867 and 353868) not 


available with low output bias ootion. Output beas voltage range: 12-17V. 


Dimensions sown in inches exceot where noted 


10-322 Thd Coanial 

Connector 
\ a ae 50 
Ora 


ae 


] \ 


a 


t 


ol 





172 % 10-32 The 
Hex Mig Hole 


353802 & BOS 


GENERAL PURPOSE —___ 


Ty Rw «4 


Oy « Le Rt en} | = | Ae VL” oe ee 
By rier hte § Fi See ole | oe ate otllg By Aacasaa BS 























Ss ES BS 


[0 [500] 280 [500] 
[anes | _o0r | 0.005 [001 
[E59 +250] G5 to «280| 6510-250] 05 10 +250) 
[~s003 | 003 | soos | 2003 | 
ee ee ae | 
(eee ee cs ee | 
[181930 | 16090 | i890 | 1810390 | 
[220 | 21020 | 2120 | 20020 

[$100 






















770-32 coax | 10-32 oomx | 10-32 comx | 10a2 coax] 


19-32 







191038) | 101035) | 103s | 1010.35) | 
VIII EAE BARRE 
PNA [NA [NA | _NA | 


TBeseor | ess603 | ossecz | as9608 | 
“assapor | As3003 | ~A83602 | AIEI004 | 
Pwasse01 | was60s | wasoec2 | 

IDSA DEOL NCEEEL. 
[080A26 | 080A26 | 080Az« | OB0AZe | 
EG A SIEER DTe Git 
LILLE ELA DE CPR ETAL 
Cn CC 





















E Nigh frequency/iow mass sensors will have 2 5 sec. crecharge lime constant 
{typrcaity B seconds) with low frequency (‘CO’) option. 

7. See cables/accessores section beginning on page 78 for a complete 
description of the supplied and optional accessones. 

8. NIST taceadle caliiration certiicate su0plied with each accelerometer 

8. Supplied mounting shud is threaded 10-32 to 10-32. 

10 Models 0353643, 0353844, 0353851, O3538S2 have a maxmum 
tormerature of 1BO°F (82°C). 
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APPENDIX B5. MANUFACTURER SPECIFICATIONS FOR WATERPROOF 
TRI-AXIAL ACCELEROMETER 
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APPENDIX Bé. MANUFACTURER SPECIFICATIONS FOR BOWLING BALL 
UNI-AXIAL ACCELEROMETER 


WILCOXON 
RESEARCH 


Models 352/352T 


(seneral Purpose, Charge Output Accelerometers 








INOUSTRISS Sty PAROD 
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© Uuehtweigh! +3. Ste mies teen ste *O- 15,000 Fiz 
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’ 1 ELECTRICAL 
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APPENDIX C. DISCRETE MODE SOURCE TECHNICAL DRAWINGS 


This appendix contains the technical drawings for the various source-to-ground 


coupling mechanisms designed during the course of the research. 
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CONCEPT SOURCE TECHNICAL DRAWINGS 


APPENDIX Cl. 





ssununopy Joxeys 


te 


wOO'E SLO 
‘a 
deL .Z€/01 


jURO|W 
JN}OWIOIBFAOOY [eIXeV], “MOA GOL 





. \ ssununo|py 


JOYEYS 





del .Z€/01 dey .Z€/01 
JONOP] JoaUIOIJOPOOOY [eIXeN |, JUNOP JoIWOIIJaIYy exe] 


SMALA OPIS ‘MATA WOL 


106 


ee 
2 yaar > 
VINE Dt 


b aus 
od 
in} 


tee 
% 
As 3 


be 


+t 2 
& 


¥, 4 By om . % Z ‘ y at wi a e @ 
Utara & Oe , Seiten Ware Saale eel Ae e f “wh ff,” te oa 
a 

& 


ee © Sie 
fee, UT Tee = 5 
EOWA Lt S er 


“ 7 - > K we Dae % 
— ~ * ee ? “~~ . vce re Oe 


Lee tee Ecole, inate 





APPENDIX C2. FIELD TEST SOURCE — MOD I TECHNICAL DRAWINGS 
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Side View: 
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Top View: 


Top Piece: 





Triaxial Acclerometer Mount 
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APPENDIX C3. FIELD TEST SOURCE — MOD II TECHNICAL DRAWINGS 
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APPENDIX C4. FIELD TEST SOURCE — MOD IV TECHNICAL DRAWINGS 
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Side View: 
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APPENDIX D. TEST EQUIPMENT WIRING DIAGRAMS AND 
EXPERIMENTAL SETUPS 


This appendix contains the test equipment wiring schematics and diagrams of the 


experimental setups utilized during source experimentation. 


iy 


APPENDIX D1. CONCEPT SOURCE LABORATORY TEST EQUIPMENT AND 
SETUP 


The following diagrams provide the wiring schematics and tank layout for the 
laboratory test equipment used for concept source experimentation. The equipment setup 
consisted of two units, a signal generating rack, and a signal receiving rack. The signal 
generating rack consisted of two phase-locked function generators, with the master 
function generator controlling the left Mini-Shaker and the slave function generator 
controlling the right. The output from the function generators was sent through a voltage 
amplifier then to the two vibration sources, with the resulting signals displayed on the 
oscilloscope. The pre-amp provided power for the source mounted tri-axial 
accelerometer with only the longitudinal and vertical components as outputs to the 


oscilloscope. 


Signal Generating Rack: 


HP3314A HP3314A 
Master Function Slave Function 
Generator X O | Generator 







Techron 7520 Left Mini-Shaker 
ike Tri-Axial Accelerometer 
Amplifier Right Mini-Shaker 
PCB 482A17 
Accelerometer 
Philips ae ee 
Channels PH3384 Channels | 2%¢ Fre-Amp 
l 2s a4 Digital O-1-x 


x xX 





Oscilloscope Ow xX ee 
eee 


03x 
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The signal receiving rack consisted of a unit built by ARL:UT, which provided 
power to the two tri-axial geophones and output connections for each of the three 
components of motion. The three oscilloscopes displayed the received signals from the 


two geophones. 


Signal Receiving Rack: 


Ghanrcle Oscilloscope 
ix tome *x’ Component 










Channels | Oscilloscope Geophone 
ix 2x | Y Component Receiver 
CE as 
Channels | Oscilloscope Geophone #1 
i coe *z’ Component 
Geophone #2 
a 


The source and geophones were positioned in the test tank as shown below. 


[P 


Tank 
Walk-way 
[- 60 cm c= 80 cm at : 


/ RSV 


Sand Tank Setup: 


Sand Filled 
Redwood Tank 


Support 
Brace 







Concept 
Source 
Geophone #1 


Signal Receiving 


Signal Generating 
Rack Rack 


APPENDIX D2. FIELD TEST SOURCE —- MOD I LABORATORY TEST 
EQUIPMENT AND SETUP 
The following diagrams provide the wiring schematics and tank layout for the 
laboratory test equipment used while conducting source experimentation with the Mod | 
Field Test Source. The equipment setup and tank layout shown below are the same as the 
setups described in Appendix D1, except for the addition of two notch filters in between 


the function generators and the voltage amplifier in the signal generating rack. 
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APPENDIX D3. FIELD TEST SOURCE — MOD I SOFTBALL FIELD TEST 
EQUIPMENT AND SETUP 


The following diagrams describe the equipment and experimental layout utilized 
for the field experimentation of the Mod I Field Test Source. The equipment setup is 
nearly identical to the setup as described in Appendix D2, except that the two Bass 


Shakers were now operated independently. 


Signal Generating Rack: 
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APPENDIX D4. BASS SHAKER FAILURE TEST EQUIPMENT AND SETUP 


The test equipment used in investigating the operating limits of the Bass Shakers 
is shown below. For the necessary testing, only one side of the signal generating rack, as 
described in Appendix D3, was required. The HP35665A Signal Analyzer was used to 


monitor the output signa! of the case mounted uni-axial accelerometer. 


Bass Shaker Failure Test Equipment: 
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APPENDIX D5. FIELD TEST SOURCE — MOD II SOFTBALL FIELD TEST 
EQUIPMENT AND SETUP 


The equipment diagram shown below is identical to the setup shown in Appendix 
D3, except for the addition of the SD390 Signal Analyzer. Incorporation of the SD390, 
allowed for acquisition of eight data channels. The dotted lines in the diagram are used 


only for clarity of the test equipment setup. 
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APPENDIX D6. FIELD TEST SOURCE — MOD ITI FIELD SETUP 


The diagram shown below describes the experimental layout for the final ground 


phase testing of the Mod III Field Test Source. 
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APPENDIX D7. GROUND TRUTHING TEST EQUIPMENT AND SETUP 


Ground Truthing Equipment Diagram: 
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APPENDIX D8. FIELD TEST SOURCE —- MOD IV BEACH SETUP 


Beach Equipment Diagram: 
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APPENDIX E. CONCEPT SOURCE EXPERIMENTAL RESULTS 


This appendix contains the qualitative results for the effects of frequency and 
phase on the axial surface wave field for the concept source. The effects were 
characterized by the amplitude of the longitudinal and vertical signals of the geophone 


farthest from the source compared to those of the geophone closest to the source. as 


observed on the oscilloscopes. 


Frequency: 100 Hz 


Case 1 
Radial: No Noticeable Difference 
Vertical: Visible Reduction 


Case 3 
Radial: Visible Reduction 
Vertical: No Noticeable Difference 


Case 5 
Radial: Visible Reduction 
Vertical: Visible Reduction 


Frequency: 150 Hz 


Case | 
Radial: Visible Reduction 
Vertical: No Noticeable Difference 


Case 3 
Radial: Visible Reduction 
Vertical: Visible Reduction 


Case 5 
Radial: No Noticeable Difference 
Vertical: No Noticeable Difference 


Case 2 
Radial: No Noticeable Difference 
Vertical: Visible Reduction 


Case 4 
Radial: Visible Reduction 
Vertical: No Noticeable Difference 


Case 6 
Radial: No Noticeable Difference 
Vertical: Visible Reduction 


Case 2 
Radial: Visible Reduction 
Vertical: Visible Reduction 


Case 4 
Radial: Visible Reduction 
Vertical: No Noticeable Difference 


Case 6 
Radial: Visible Increase 
Vertical: No Noticeable Difference 


Frequency: 200 Hz 


Case | 
Radial: Visible Increase 
Vertical: No Noticeable Difference 


Case 3 
Radial: No Noticeable Difference 
Vertical: No Noticeable Difference 
Case 5 
Radial: No Noticeable Difference 


Vertical: No Noticeable Difference 


Frequency: 250 Hz 


Case | 
Radial: Visible Reduction 
Vertical: No Noticeable Difference 


Case 3 
Radial: Visible Reduction 
Vertical: No Noticeable Difference 
Case 5 
Radial: No Noticeable Difference 


Vertical: No Noticeable Difference 


Frequency: 300 Hz 


Case 1 
Radial: Visible Reduction 
Vertical: Visible Reduction 


Case 3 
Radial: Visible Reduction 
Vertical: No Noticeable Difference 
Case 5 
Radial: No Noticeable Difference 


Vertical: Visible Reduction 
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Case 2 
Radial: 
Vertical: 


Case 4 
Radial: 
Vertical: 


Case 6 
Radial: 
Vertical: 


Case 2 
Radial: 
Vertical: 


Case 4 
Radial: 
Vertical: 


Case 6 
Radial: 
Vertical: 


Case 2 
Radial: 
Vertical: 


Case 4 
Radial: 
Vertical: 


Case 6 
Radial: 
Vertical: 


Visible Reduction 
No Noticeable Difference 


Visible Reduction 
No Noticeable Difference 


Visible Reduction 
No Noticeable Difference 


Visible Reduction 
Visible Reduction 


Visible Reduction 
Visible Reduction 


Visible Reduction 
No Noticeable Difference 


Visible Reduction 
No Noticeable Difference 


Visible Reduction 
Visible Reduction 


No Noticeable Difference 
Visible Reduction 


Frequency: 500 Hz 


Case ] 
Radial: No Noticeable Difference 
Vertical: Visible Reduction 


Case 3 
Radial: Visible Reduction 
Vertical: Visible Reduction 

ase 5 
Radial: No Noticeable Difference 


Vertical: Visible Reduction 
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Case 2 
Radial: Visible Reduction 
Vertical: Visible Reduction 


Case 4 
Radial: No Noticeable Difference 
Vertical: Visible Reduction 

Case 6 
Radial: No Noticeable Difference 


Vertical: Visible Reduction 





APPENDIX F. AMBIENT BACKGROUND NOISE ANALYSIS 


This appendix contains the MATLAB program and output figures used to conduct 


the spectral analysis of the ambient background noise at the beach test site. 


oo 


Filename: noise2geos.m 

Written by: F.E. Gaghan 

Date Last Modified: 27 January 1998 

Purpose: This program is used to calculate the frequency spectrum and 
power spectral density of the background noise from one of 
two geophone signals for the various test sites. The user 
has the option to perform the analysis for a single data 
file or to perform an incoherent average for a number of 
data files. 


A of oP oO of oD oO 


fe) \o) 


Initialize user interface 

average=input('Do you want to average the noise plots 
Keor 1s) 7 Se 

directory=input (‘Input the noise file directory: ','sS'); 

Pitemed HeCEOLy— | clnoeeery |); 


ae 


if average=='y' 

mum averages; 
else 

putea Sorage— 1, 
end; %Send-if average 


$ Assign noise Signal component to appropriate data channel 
geophone=input ('Do you want the data from geophone 1 or 2: '); 
if geophone== 
revr_ horiz column=4; 
revr vert column=6; 
else 
GevE shore. columni—o; 
‘Cvmevert yecolumn—7/; 
end; send-if geophone 
sampling rate=5120; 


$ Loop structure used for signal averaging 
Ct. teers ton mum ravemagc 


chise e'); 
if iteration== 

file=input('Input the noise filename: ','s'); 
else 

file=input('Input the next noise filename: ','s'); 


end; %end-if iteration 
filename=(file]; 
évald(* load <=" f£ilexdireceory, \ 7titename |, 


heorizyne1se—-yres90 (+, eevr orl 2. colunn = 
Mean (<res 90 rey rE shonlzZ.COLluma yy, 

vert noise=xrc390(:,rcevr_vert_column)- 
mean (xrc390(:,revr vert column)); 
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$ Calculate frequency response 
Eee lengen—Veng en ine wt emo Use i, 
freqs HzZ=sampling vate* (0: (fft lengen/7) ye aleacien 
f£t_horiz noise-abs(fit (horiz noise, Sit lemeiin) 7 e iesbom opens 
fic horiz_ noise=fft horiz _noise/max(fft_ horiz noise); 
fit veripnoise. ie Gaaencoce noise, fft _length))/ff£t_ length; 
fft vert no1se= fft_ vert noise/max(fft vert noise); 


$ Calculate power spectral density 
psd horiz moise—psd(horiz noise, fit Vength sanpieigg oe) 
“(2/7 samp lange rate ) 7 
eo horiz noise= ai acchny ee horiz noise); 
horiz density amp= IO 5 el Sheksro | | ate neak 74 _noise/10); 


pedgver: MNO1se—psd (vert nolse, & Ee arengcen, samo timepiece) 
*(2/sampling rate); 

dBpsdever’e moise—-10” logl0tpsd verermomse); 

vert density amp=10.%*(dBpsd_vert_ noise/10); 


oo 


Perform incoherent average of frequency spectrum and power 
spectral density 
if iteration== 
hOriz Spectrum=zerces (si ze (tit Were z ene tocren 
horiz SpeGtrum=horizrspecr rum: ee noneze else. 7m, 
Vertn spec’ hUM=Z76res (seca eerie tse) ye 
Difeyaye cj of=lele te Wl if=iein us) ola(eene blues (gic i eieie elo able 3 |) 
NOrtiz dens fry Zeros (Siew nomi _density_ amp) ); 
horiz density=horiz densitythoriz density amp; 
vert density=zeros(size(vert density amp) ); 
VErepdensSilt y-Vere Oeislaytyere Cote Teyeaine: 
else 
hori z SpecCerum—NordZ spec er UM (a me mer gem Cdn mee, 
VELENSPCCErUM{=Vere Sspecerannn ( Phe. Jomtien iewse racy, 
NOrEZGensliEy=norizg —CenSsiey hors Zeceno ie mame: 
VGlLe Gemsley—-Vere Cele tev ome pichHo ley mater 
end; %end-if iteration 
end, »t2°enud=ter 1ceration 


fa) a 


Identify maximum frequency of interest 
freq marker=0; 
for i=l:length(freqs_ Hz) 
if freqs Hz(i)<=s500, treq manker—irequmarken aed 
end; tend ior: a 


Normalize frequency spectrum 

avg NOriZ SpPeCUruUmM—Sdme (sori 2Z spectrum: 

avg horiz spectrum=avg horiz spectrum/max(avg_horiz spectrum) ; 
avg Vert, Spectrum=sqrt (Verte spectrum) ; 

avg vert spectrum=avg vert spectrum/max(avg_vert_ spectrum); 


if iteration== 
avg NOrlZ demsit y=ner i Zacdems mia; 
aVGLVErrECcenest y—-Vert Cenmetun] 
else 
avg horiz density=horiz density/3; 
avg vert density=vert density/3; 
end; %end-if iteration 
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fo\\°) 


ae 


oie 


ave 


JBa7 Gone zecdensit.7—10* log l0 (avquhourmzedecisiey). 
teaver cecensily—-l0*loglO0(avg vere density); 


Plot normalized longitudinal frequency spectrum and power spectral 
density 

figore (1); 

suep lot (211); 

PICEterecs EZ(2.treq marker) ,avguhort az spectrum(2: freq maken. 
title('Average Longitudinal Background Noise'); 

xlabel('Frequency [Hz]'); ylabel('Relative Magnitude'); 
Subp ikee (2 2) 

plete theqemiz(2 treq Marver), dBav gue ere eelistEyi 2 freq Marker) ) 
xlabel ("Frequency [Hz]"); ylabel Spectral Density [dB] *); 


Plot normalized vertical frequency spectrum and power spectral 
density | 

jew ebb a — 5 2) er 

Subei oe (211); 

plot (ireas izi2. freq Marker) ,avg Vert Speeenume are coma ker i, 
title('Average Vertical Background Noise'); 

xlabel('Frequency [Hz]'); ylabel('Relative Magnitude'); 
subplomiz 12); 

plouyErecseiae steeq marker) ,dBavg vert density(2.2req mar ker), 
xlabel('Frequency [{Hz]'); ylabel('Spectral Density [dB]}'); 
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APPENDIX G. GROUND TRUTHING 


This appendix contains the various MATLAB programs written to analyze the 


ground truthing data collected. 


HW WP WO WO 


oP Wf oO 


foe) 


are 


ae 


Filename: truthingdisplay.m 

Written by: F.E. Gaghan 

Date Last Modified: 14 February 1998 

Purpose: This program is used to graphically display the evolution of 
the wave packet generated by the dropping the bowling ball 
in both the longitudinal and vertical directions as a 
function of time and range. The Rayleigh wave speed can be 
estimated from the graph. 


Establish user interface to control data input 
directory=inputi “npucy the data gtile diveceouy. 7 9s); 
file directory=[directory}; 

freq span=input (‘Input the frequency span in kH2: °); 

Sane Singm sake rtd spam. OU 2. a6, 

Saniple adelay—i1npur ( Input Che nunber sor edelay eo otploo.. 7G 
Samp horde tay—-abe(samplosdelay 


Assign input channels to correct data columns 
Eemaqer Teor — 1, 

GSOr pPomgneOl=4; qeolevereyco!=6; 

geo2 long col=5; geo2z vert col=/; 


While loop used to control data file input sequence 
iteration=1; 
continue='y'; 
while continue=='y' 
% Input and load data file 
2 hod a Er 
if iteration== 
file=input (‘Input the data filename: ','s'); 
else , 
file=input ('Input the next data filename: ','s'); 
end; % end-if iteration 
filename=[file]; 
evalt |’ loed)* ftilerdirectory, \*, £1 lemme i), 


range ft=input(" Input. che range<in ft to tthe fresergecpione =, 
range m=range ft/3.2808; 


% Select desired geophone 

geophone=input ('Input the geophone number for that range: 
ae ee, ce 

if geophone== 
Geelong col-geol Tong col; 
geo vert col=geol vert col; 

else . 
geo long col—geo2 long col; 
deo _verescol—geoZ vert col, 

end; %end-if geophone 
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o 


~ 


o 


oo 


oo 


Adjust geophone timeseries in accordance with trigger signal 
CEigge? Lameserics -xres0(:, Lriggerncolan 
[Eragger: value, trrgqger, pt) — 

mMmin{trigger timeseries(|]:sample delayio00j 9. 
Svame PlC-trlgge mares, 
record length=length(trigger timeseries) ; 


time scale=[0:(1/sampling_ rate): ((record length- 
Stapemet)/Ssampling rate) ]', 


geo stongle  meser 6s=xrc3 90 (Start mil weeco rem homo s 
Geog loncgacel)'; 

GeGsvere timeseries=xrc3904 (stamegee- secora Lengua, 
Been verea col, 


Filter and scale geophone signals 

b=ones(1,10)/10; a=1; 

Gea leng ti ltered=tiltii leo, an g-om long eiieeseGles .70lur 
olsreN ASig egees ee a arel eat le igiibies (ery Cision cpeie. eesicteibacy) 24010) 2 
ERIgG@ger EaMme@serles= (rigger ii ie seme slr 


Plot the timeseries of the filtered geophone signals linearly 
scaled by the range to the geophone 
adjust yaxis=range m; 
range—sprintt(°%2.28 7 ranger, 
if iteration== 
hold on; 
frgure (1); 
plot(time scale, trigger Teimeser wes Serene. meeore elem aina © 
time scale, geo) lonomerreered | deg ta em yar irsys, 
text (0.7,geo long filtered(recorad Tength) tadjuse yaxise le, 


iwange, mm’ lj] 
xlabel('time [sec]'); ylabel('range [m]'); 
holeaion; 


figure ( 2) ; 

plot (time scale, trigger uimeser ies (Sstareuape> record bengunyE 
epime scale; geo Vert fl lteredtadsuceaya ts 

text (0.7,ge0 vert Ttiltered(record length) tadgustuyaxis.024a, 


[range,' m']); 
xlabel('time [sec]'); ylabel('range [m]'); 
else 
igtenkelRoiay: 


Eraqure (1); 

Pleetrime scate geo 7long) PiMitereds demise cher e, 

text (0.7,ge0 long filtered(record length) +adjust. yaxic 0243, 
[range,' m']); 

hola on. 

figure(2); 

plot (time seave, Geotvney i iirered cust yaxts)y 

text (0.7,ge0 vert filtered (record Wengi ln) fadsils tiny acer a4 or, 
[Pange.,- Sim ju; 

end; tend-if iteration 
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opemeck tor additional data file inputs 
iteration=iterationtl; 


GES rg eet, 
continue=input('Do you want to input another data file? 
(y or n) ','S'); 


end; tend-while continue 


Rayleigh Wave Progression as a Function of Time and Distance 


range [m] 





0 0.1 Gr2 0.3 0.4 0.5 0.6 0.7 
time [sec] 
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Filename: gabormain2.m 

Written by: P.Ee. Cagnan 

Date Last Modified: 10 February 1998 

Purpose: This program is used to calculate the time versus frequency 
components from one of two geophone Signals generated by a 
dropped bowling ball in order to analyze the wave speed 
characteristics of a propagating medium. The user has the 
ability to chose from three different processing methods. 
The results are displayed in the form of a gabor matrix 
(Conecouree lon) . 


Declare global variables 
global slowness velocity freqs Hz freq marker 


Establish user interface to control data input 
direstony—1eut( input the data virile directory. ", Ss ); 
filemctieecrouy— (directory); 


Assiqgm input channels to correct data Columns 
triggerscce!—|], 

gecleleng col—4> geol vert colo, 

GeoZeleng  col—>, geez. vert (cera. 


Enable user to incoherently average the gabor matrices 

ai soto 

average=input('Do you want to average the gabor matrices? 
[yacrey) 8 ys) 7 

if average=='y' 
num average=3; 

else 
NUN verage—l; 

end; tend-if average 


Select geophone signal to process 
GUsSe ee Oe 
geophone=input ('Select geophone number: 1-#1, 2-#2 '); 
if geophone== 
geo long col=geol long_col; 
gG¢o vere col-=gcolwvere col, 
else 
geo Jeng col-geoZ long col; 
geo vert col-geo2 vert col; 
end; %end-if geophone 


For loop used to control data file input sequence 
[Or 1Teeratton— ls nuteaverage 


qdiset > “je 
if iteration== 
file=input ('Input the data filename: ','s'); 
else 
file=input (‘Input the next data filename: ‘,’s'); 


end; %end-if iteration 
filename=[file]; 
eval (| load *viile directory, \"; 2blenamne li), 


150 


ne  neerati1on=—— 


range ft=input(['Input the range to geophone ', 
HDUMZser(gqeophnone)," in ft. We 

range m=range ft/3.2808; 

freq span=input ('Input the frequency span in kHz: '); 

samp linge sate—-freq Span i000" 254. 

sample delay=input('Input the number of delay samples: '); 

sample Gdelay-abs (sample delay); 


2) 


ena. s end-1rr ererat Lon 


% 


Adjust geophone timeseries in accordance with trigger signal 

EMigger Cimeseries—xres90 (i, eae mee a 

[trigger value,trigger pt]= 
mMin(trigger tameseries(]:samplemdela, 7590) 7; 

Slate gpe te eigger oe-7, 

record lengen—lengen (trigger timeseur—s) 

time. scale=|0: (l/sampling rate). | tecordmlengiis 
Seartgee) samp ling, eaba) |, 


long -timeserres—-xres90( (start, pt; recora  lengen yay Oe Pongncol), 
Vert, CLimeseries=xres90 (| (start pesncecord engen geomyert col), 


Filter geophone signals 

b=ones(1,10)/10; a=1; 

Tongue iVeewed=1ile tilt (b>, 42, Lond geese nies). 
VOGE er ebored— tattle’ (>,a, Vere ees omnes 


Plot longitudinal and vertical timeseries of geophone signals 
fagure (1); 

Spe toe( 211); 

plot (time scale, long timeseries) 


title(['Longitudinal Time Series of Geophone Signal for File: 
Pade je; 
xlabel(*time [sec] *);-=vlabeli amplitude Jil: 


Subp LOECZIZ 7 
pletltime scale, vert timeserics); 


title(['Vertical Time Series of Geopnone Signal tor File: "| 
file) j)> 
xlabel(’tamenisec]*)> ylabelt amplitudes (ea, 


Plot hodogram of geophone signals 
EVquire (2)y- 

ploe longa t witered vert of Mitened |, 
xlabel('longitudinal amplitude [V]'); 
ylabel('vertical amplitude [V]'); 
axis equal; 


Display geophone signals in a hankel plot 
RVZcexts—-zeros(lengeha(time scalej, 1), 
y_axis=linspace(min(long filtered),max(long_ filtered), 
engin (tamenscale):); 
Zz axis=linspace(min(vert filtered),max(vert_filtered), 
length(time scale)); 
figure(3); 
plets{cime jcale, long ii blered, vere filtered: 
time scale, xyz _ axis, xyz _axis,'k',... 
Dz e ey os ony re ey ee 
HVZ axl S, Ky¥Z aks) 2 ex Sy kK), 
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xlabel('time [sec]'); 

ylabel ('longitudinal amplitude [V]'); 
zlabel ("verticalvanplitude [Vv]")> 
VLEWw (lS Asus 


Enable the user the ability to display selected portions of the 
time record 


disp San, 

isolate=input('Do you want to isolate the Rayleigh wave? 
CY sOie Oi) eerie 

if isolate=='y' 


cont mue— vy"; 


while continue~='g' 
timel=input ('Input the start time of the hodogram: '); 
pointl=floor(timel*sampling rate); 
time2=input ('Input the stop time of the hodogram: '); 
point2=floor(time2*sampling rate); 
6 Plot hodogram of user selected portions of timeseries 
Ruy e (45) 3 
pieced long. Ell tered (pom l peiney., 
VEre Pil tered) pened ec rmine gy am 
long filtered (point!) yvere, fi ltereciieeiae ae. 
long filtered (pointe?) vert tileterecdhpeineZ) 7. =< e 
xlabel ("longitudinal amplizcude™ [¥]')- 
ylabel (‘vertical amplitude [V]'); 
axis equal; 


oe 


Display hankel plot of user selected portions of timeseries 
x axis—ones (pointZ-pointe byl) timescale (portant ae 
y axis—lainspace (min (long ti ltered(poimrel pointe ae 
max (long frltered(peintl: peint2) |; perme? [ee tai ae e 
Z axis—linspace (mim (vere faitered) pore pe ae ar 
max(vert filtered(pormel perme), pommes eo li las. 
freure (os); 
Pligeo (Gime sce (ecint Lape iie 
long filtered(pointl:point2), 
Vere i Pil wered (Omi emi 2 yee. 
time, Seale (point!) longgrrltcred (porn iy 
dete ee eat dbictsagi=(ejeleatoue Ih) ye Gig ae 
tame Seale (poine2);) tong git tered poraec | 
Veto rie heme d (Om TA) yu sae ee 
time scale (point! point2),xyzeexis (peimel: pemmezas 
MY 2, xo Pomme: Ome 2) ee 
x axis) ders, Hy 2 axis point porme2) 76 ee: 
x AXIS, V2, ax uo perme ll: peinE 2 GZ soe 
Xxlabel('time [{sec]'); 
ylabel('longitudinal amplitude [V]'); 
Zlabel(‘'vertical amplitude [V]‘'); 
axis ( [time scale(peintl) time Scale (pean. 
mim(long filtercd) max (long. tw Porcar 
min(vert filtered),max(vert filtered) ]); 
view (15, 30); 
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end; 


oO 


Graphically display ratio=ot longitudinal tojverticet 

Signal amplitudes 

ell pse-=1nput (' DO yor wamesto calculate ratio Gb Minorca! 
axes<e (Y OF n)t. 'see oe 

if ellipse=='y' 
long _hilbert=abs (hilbert(long filtered (point1:point2) ))}; 
vert hilbert=-abs(hilbert(vert filtered (pointl:point2))); 
amp ratio=long hilbert./vert hilbert; 


fo)\eJ 


Pi giiee GC); 
Suppo lot (2 Li) 
PEOECeEime Scale (poimel= sly aeotmes fo 7 

Long filtered (peinul= sie. sointZrol2). oe 

time scale (pointl—olZscetartZ oi) 

Ver fa Liered (poimel ol worl 22002) 4. ==) 
xlabel('time [sec]'); ylabel('amplitude [V]'); 
Vegend (* longitudinal”; veretesim aa, 

Subp lot (202); 

eB OE(rame scale (poine |2 poime a romp amce to), 

xlabel('time [sec]'); ylabel('longitudinal/vertical'); 
end; send-if ellipse 


continue=input('Press any key to continue isolation or "gq" to 
Gitar es ) 4 
end; %end-while 
end; tend-if isolate 


Select data processing method 
if iteration== 
arse ( 9 )a 
processing method=input ('Select processing method: 1-slowness; 
2 elOCLey; S2o eranioroums. 
end; tend-if iteration 


BLE Jenmielo ssyenole; julie) ec 
gabor matrix=gaborslowness(file,long filtered,sampling rate, 
Bene eent)=, 
elseif processing method== 
gabor matrix=gaborvelocity(file,long filtered,sampling rate, 
range M); 
else 
S matrix=stransform(file,long filtered,sampling rate,range m); 
end; tend-if processing method 


Incoherently average the gabor matrices for three data files 
amplitude matrix=10.*(gabor matrix/10) ; 
if iteration== 
amplitude average=zeros (size(gabor matrix) ); 
amplitude Vaverage-ampli tude average+ampi 1 tude macmix, 
else 
amplitude average=amplitude average+amplitude matrix; 
end; %end-if iteration 
Send-for iteration 


Q 


fo} 


Plot incoherently averaged gabor matrix 

if average=='y' 
amplitude average=amplitude average/3; 
gabor average=10*1logl0 (amplitude average); 
dbmax=max (max (gabor average) ); 
contour “levelsaiaena. |: dbmax—6ie 
figurei7); 
CONNEC Eregouna2 hed Marken) silewmess, 


GabGr Matrix ereq Marker, :))) Comecutegeas 7 -weorm 
title('Average Propagation Speed as a Function of Frequency'); 
xlabel('frequency [Hz]'); ylabel('slowness [sec/m]'); 


end; tend-if average 


Enable user to plot a fitted group velocity curve to gabor matrix 
Gisp (aa 
group=input('Do you want a group veolcity curve for the gabor plot 


(Y OL jr 7. s~ |); 
if group=='y' 
gaborgroup 


end; %end-if group 


Enable user to cross correlate the geophone signals 

gasp.) 7 

correlate=input('Do you want to correlate the geophone signals? 
CV Or |) ys ye 

if correlate=='y' 
Toric Pe mes eigkes =a) 


({StCart pt. record Tengen) geo Wenge es 
iemgZ2 se eieser tes me] Oy 

( 

( 


Start pt: necerd Tengen), doo 2s lende-a 
start pt: record lengrh),;geet™ rere lr 
Start pts record lenge) -dqece7 ceeace ss. 


Wésicic IL Veabuiselaicaket cic 8210) 
vert2 timeseries=xrc390 


ee 


b=ones(1,10)/10; a=1; 

Vong! faltered sti i tilt oa eng eee imac cmre se 

long2 filtered=f1iltfilt(b,a,long2 timeseries) ; 

Tongl2 “ceonr—xcorn( long! filtered, long2 pe tii@eree, eee: 
longl2 hilbert=abs (hilbert(longl2 corr)); 

vert] filtered=filtfilt(b,a, vertl timeseries) ; 

vert2 filtered=filtfilt(b,a,vert2 timeseries) ; 

vertl2 corr=xcorr(vertl filtered,vert2 filtered, 'coeff"); 
vertl2 hilbert=abs(hilbert(vertl2 corr)); 


$ Plot the cross correlation of the longitudinal geophone signals 
figure (8); 
plot (long l2 sieeise nt. 7, 
title('Longitudinal Cross Correlation'); 
axis ([length (long! filtered), fongtn Peng! 2a Voenre 
O,max(longl2 hilbert)]); 


ac 


Plot the cross correlation of the vertical geophone signals 

LIguUre ioe, 

plot vertl2 hilbert)-; 

title ("Vertical Cross. Cerrelaticnu. - 

axis([length(vertl filtered), length(vert12 hilbert), 
O,max(vert12 hilbert)]); 

end; %end-if correlate 
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Filename: gaborslowness.m 

Weweeen soya Foe. Gaghan 

Date Last Modified: 10 February 1998 

Purpose: This program is a subroutine used to calculate the time 
versus frequency components of a geophone signal using the 
element location within the timeseries record to determine 
the time of arrival which is converted to slowness (inverse 
speed). The subroutine returns the gabor matrix to 
gabormain. 


function gabor matrix=gaborslowness (datafile, timeseries, 
Sapling 2Leqd, separation) 


Declare global variables 
global slowness freqs Hz freq marker 


Initialize control variables and storage vectors 
Eft length=1024; 

ese cuptoles<= Ile ie ele e irate ee) 

(Pee licde x tase nde x— (htt bengtni- i), 

fft overlap=fft_length/16; 

gasor matrex—[j7 

slowness=[]J; 

Ps 


Conduct time versus frequency decomposition of geophone signal 
tem incdex—litiproverlapstft index=] 

blockimidpoint=index+(frt Lengtn/2)—1; 

arrival time=block midpoint/sampling freq; 

slow=arrival time/separation; 

slowness=[slowness, slow}; 


Pep bloeck—|(tameseries( i idex tidex a tim ron mem 
bescdcetr block, treqs (na) —psdi tit plesk) theplengtih, sampling gh ecyT, 
Gen tee olock— LO loglu(csdyEertolock); 
gabor matrix=[gabor matrix,dB fft_ block]; 
By bap alles 
end; end for index 


Determine element addresses for major frequency components 
freq marker=0; 
for i-l:tengeh (t@regs Hz) 
if freqselz(i)<—100). freq marker= req markers] wend 
end; %end for 1 


Plot the contours of the gabor matrix 
dBmax=max (max (gabor matrix) ); 
contour levels=(dBmax:-1:dBmax-6] ; 


figure (7); 
Gontour(freqs Hz(2:freq marker), slowness, 
Gabor Mates (2: Preqamarker,:) contour levels, 
title(['Propagation Speed as a Function of Frequency for File: ', 
datafile]); 
xlabel('frequency [Hz]'); ylabel('slowness [sec/m]'); 


return: %tend-function 
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Filename: gaborvelocity.m 

Written by: F.E. Gaghan 

Date Last Modified: 10 February 1998 

Purpose: This program is a subroutine used to calculate the time 
versus frequency components of a geophone signal using wave 
speed to identify the associated element location within the 
timeseries. The subroutine returns the gabor matrix to 
gabormain. 


function gabor matrix=gaborvelocity (datafile, timeseries, 
Sampling freq, separation) 


Declare global variables 
global velocity freqs Hz freq marker 


Initialize control variables and storage vectors 
fft length=1024; 
max velocity=200; 
min velocity=50; 
Vel@cieye lac remene=4. 2; 
Puno Behe renenus— 
floor( (max velocity-min_velocity) /velocity increment) ; 
Gabor wmMatrix—([1, 


Conduct time versus frequency decomposition of geophone signal 
for index=0:number increments 
velocity (index+1)=max_ velocity-(index*velocity increment); 
block midpoint=floor((separation/velocerty (index!) “Sanp lang sea ape 


if block midpoint < (fft_length/2) 
padding=(fft_ length/2)-block midpoint; 
Zeno padding —-zeLos(padcinc |, 
block vend-block miGgpoinut(Eruelenguay. i i 
tft bleock—|zero padding; timesewicea Vo leec send: 
else 
block Start-bilock midpoty es (paeehencirn 7) ne, 
block _end=block midpoint+(fft_length/2) ; 
if block end > length(timeseries) 
padding=block end-length (timeseries) ; 
zero padding=zeros (padding, 1); 
occa jo lise) <— 
[timeseries (block start: Tengen (timeserites) ))) zeroueadar ie 
else 
Le Dloek—eimescrl es (block Usbavmre heel somal |; 
end; send-it ebock ond 
end; send-if block midpoint 


[psd fft block, freqs Hz)=psd(fft block, ttt Tengen, cane iar cears 
dB fft bleck-10" loglo\psd tte oock iy 
gabor matrix=(gabor matrix,daB fit block]; 

end; tend-for index 


Determine element addresses for major frequency components 
freq marker=0; 
Pons leng en ( freqs jz) 
if freqs Hz(i)<=100, freq marker=freq marker+1l; end 
end; ~send=-tor i 
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Plot the contours of the gabor matrix 
dBmax=max (max (gabor_matrix)); 
COneeurelevelo—{CRma se cbia ae: 
figure(7); 
Gemm-our(treqs Hz(2:freéq marker), 
mel Oca ly, Gabor Maur ix (2: (eequmamera.) ,cConecur levels) 
title(['Propagation Speed as a Function of Frequency for file: ' 
Cia Gann oie; 
xlabel('frequency [Hz]'); ylabel('velocity [m/sec] ''); 


f 


return; %end-function 


Filename: gaborgroup.m 

Written by: F.E. Gaghan 

Date Last Modified: 09 February 1998 

Purpose: This program is a subroutine used to fit a group velocity 
curve to the gabor matrix plot created by gabormain. 


function gaborgroup 


Declare global variables 
gl@ea Pwrrequspeed Eype fit 


Establish minimum user interface 

libs yer) 

disp('Pick off the group speed points from the gabor plot with 
the mouse'); 

disp('Press return when done'); 

{ freq, speed] =ginput; 


Determine desired curve-fitting routine 

type: Tt inpur( Pick, Ehe«type. Orsrit sel linear Jorn oeee, 
fit steps=100; 

ewer tred—rued (1 ).; 

Upp ee@eered——req (lenge (i recd).)., 

freq step=(upper freq-lower freq)/fit_steps; 
frequencies=lower freq:freq step:upper freq; 

fie Svepse=lenguh( ireq step); 


Calculate fitted curve values 

Pie peeeit=—) 
speeds=interpl (freq, speed, frequencies) ; 

else 
speeds=spline(freg, speed, frequencies) ; 

ends) send-it tCypereic 


Plot fitted curve 
hola on; 


plot (frequencies, speeds, '--k'); 


return: ¢end-—function 


ey, 





APPENDIX H. DISCRETE-MODE ANALYSIS 


This appendix contains the MATLAB programs written to analyze the data 
collected during the beach phase experimentation with the MOD IV source. 


Filename: filtereddisplay.m 

Written by: F.E. Gaghan 

Date Last Modified: 09 February 1998 

Purpose: This program is used to graphically display the evolution of 
the wave packet generated by the field test source in both 
the longitudinal and vertical directions as a function of 
time and range. The Rayleigh wave speed can be estimated 
from the graph. It also displays the associated hodograms of 
each Rayleigh wave. 


oP oP oN NW oP GO oo OP 


ro)\% 


Establish user interface to control “datagmmput 


fo)\9) 


directory=input (‘Input the data file divéeéctory: ",'s'); 
file directory=[directory]; 
drive freq=input ('Input the shaker drive frequency in Hz: '); 


freq Spaneinpuc(’ Input the frequency seanvtnecde. . 
Sampling (sate=freq span~1000~2. 56; 
sample delay=input ("Input the number of delay samples, )- 


record length=8192; 
time scale=[(sample delay/sampling rate):(l/sampling rate): 
((recora, lengtun- sample dellay=1)/ samp linagg@e es ee, 


oo 


Assign input channels to correct data columns 
acce long coll; acco vert eol—Z- 
geol long@meet—2; “geo! vertycol—7, 
geo2 longqueol-—6,) geo2Z vert aeol—a; 


oo 


Create band pass filter centered at the drive frequency 
pass freq=drive freq/(sampling rate/2); 
[DAleghpass,a Nilchpass | pultten (2, pass areq, sian. 7, 

ib Nowoass a lowoass |= oueeer(Z pass (req): 

fre length 02d, 


oP? 


While loop used to control data file input sequence 
iteration=1; 
continue='y'; 


while continue=='y' 
% Input and load data file 
Cats p65 Ta; 


if iteration== 
file=input ('Input the data filename: ','s'); 
else 
fide=1 np Inputs the mMext data filename mae. an 
end; % end-if iteration 
filename=[file]; 
eval { load.” ,fule directory,” \ .eisenamel), 


Gangestt—inpue ll npub ene rangs in@ft to Ehe eirst gecpnoge. 7, 
range m=range ft/3.2808; 


is? 


3 


= 


° 
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Select desired geophone 
geophone=input (‘Input the geophone number for that range: 
See 2 ea 
if geophone==1 
gee long Jeo Jeo peng col; 
offebe es ice IE cpelo i chonmemere den 
else 
Géo signgmeo—Jeo7 bong col; 
geo vert col=geoZ vert col; 
end; tend-if geophone 


Filter accelerometer signals 

ae ecetmrer ye ameserles—=xr e390 |. acee Pm eicmec as, 

jeCeleverEwrimeserires— HES Olln patecell Vere col, 

accel long. filtered=filtfilt(b Ciagnpass, a_ highpass, 
eecee long trmeserives i; 

acee ly vyert filtered=f[1ltiilt( Ss nignpacs aunties, 
gece evVert selmeser les), 

accel leng fi ltered-1iltiiie (bp teweass oho ecoc, 
accel Vong i yiterec)) 

accel vert filtered=filtfilt(b lowpass,a lowpass, 
accelovert: ti lvered)? 


[power =eeags ni) 
psedtaccel Tong filtered, record engin, sane eran ie 


Determine element address for drive frequency 
Woe ee, ieugldene 9); 
for et) bengen (ee aemizy 
if freqs Hz(i)==drive freq, freq marker=1i1; end; 
end; Send for i 


Normalize filtered accelerometer signals by magnitude of 
accelerometer fft at the drive frequency 
accel long fft=abs(fit(accel long filtered)/fftt length); 
accel vert fft=abs(fft (accel vert filtered)/fft_ length); 
accel long filtered= 

accel long fi ltered/accel lenge f(r ecm ma icc ise, 
accel vere rritered= 

accel: vert filtered/accel vert lrrtt meq mam cers, 


Filter the geophone signals 
geo long timeseries=xrc390(:,geo long col); 
geo vert, tameseries=xre570 (|. gecuven geo, 
geo long filtered= 
Filtfilt(b highpass,a highpass, geo Tong geimeser tes? 
dec vere ti itenea 
filtftilt(b highpass,a highpass, geo vert seumeser les), 
gecn long tileerca. 
filtTilt(b lowpass, 4a lowpass, geo ulong iiitered, 
Geo Vert si liered— 
Filtftilt (5 lowsass;a lowpass, geovvert ii ltemecdi:, 
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Normalize the filtered geophone signals by magnitude of 
accelerometer fft at the drive frequency 
geo long ttl tencd= 

l.79* (geo long tiltered/accel Vong fit(freq marker) ); 
geo vert filtered= 

l /o "(Geol vert ill tered ae-c Wa ee tit (ireq marker); 
Geos lenge ti sere 

I 13 (geo long, EVitered/ Tai cmolcmt Tere): \; 
Geos ver ine ee rce.. 

ee (OSOnvert el tered) Mara igeem ele tered) ) 7; 


Plot the timeseries of the filtered geophone signals linearly 
scaled by the range to the geophone 
adjust yaxis=range m; 
range—speime kl (+2226) Langemm) s 
if iteration== 

holo (on, 

frqure (1); 

plot(time scale,accel long filtered timeyseale, 

geo lend. filteredtadjust _yaxis); 
text(O-7,gc0 long filtered (recone Wendel nods ioemy a mono, 


[frange, ° om | )7 
xlabel('time [sec]'); ylabel('’range [m]"); 
hole jon, 
figure(2); 


Plot (timegscale, acee! Verte (tered tances eale, 
deonvert: _filtered+tadjust ~Vaxts)7 
Cext( 02/7 Geo vere, Filtered pece tel Siig ei mami oie le eee 
[wange,. m1); 
xlabel('time [sec]'); ylabel('range [m]'); 
else 
hee ron; 
Paveniics (1); 
plee (time scale,geo long) filtered ragga 15), 
terw(0.7,ge0 long filtered (record yanmar sta yas 10-47 
[range,’ m']) 
he fason; 
figure (2); 
ploe(eime scale, cco vere Ellterec™ scg Mees 
Eext4O0.7,0ge0 Vert qtilterea records tengery a ccmis omy a els or og 
[range,' m']) 
end; %end-if ecation 


Plot hodograms of filtered signals 
Of fset.o ls; 
if iteration== 
Eaume dt) sa= 
input (‘Input the start time of accelerometer hocogram: ‘); 
pointe ya loor(eime! parsanp ling Gate) 7oetocr, 
tine 2 a= 
input ( Input, the Stop time of accelerometer noaegran.= |); 
point2 a=floor(time2_a*sampling rate) toffset; 


timel g=input('Input the start time of geophone hodogram: '); 
pointl gG-tlGor(timel G=sampeling rate) +ortser; 
time2 g=input ('Input the stop time of geophone hodogram: '); 
pe nee o- tleer (lime 2. same ling share} Orrsot, 


16] 


range {[m] 


—* 
(om 7 


heaven; 
figure (3)" 
plot (accel longer tileered (point! esperar 
accel Vere fa tered (poime lea: Celta aa 
geo Jongm ri leered (polit lo perme saa. 
geo_vert silbtered (point! ¢:poimt2 Oj 2cn uous a 
xlabel('longitudinal displacement"); 7 
ylabel('vertical displacement'); 
axis equal; 
else 
Eimel G-inpul, input the start Lime of geophone. hecogc am, mee 
pointl g=floor(timel g*sampling rate)+offset; 
Eime7eg— 0pue( Input the “step time cfegecesnoene Modegtan ma. 
point2 g=floor(time2 g*sampling rate)+toffset; 
holdson- 
Pure (3); 
ploeigeo Jong fritered\ pointl . poantZ soi 
geo vert filtered (point Wieq {pein ? iG) Soe cemya se 
xlabel('longitudinal displacement'); 
ylabel('vertical displacement'); 
axils equal; 
end; %end-if iteration 


oo 


Check for additional data file inputs 
iteration=iterationtl; 


ES); 
continue=input('Do you want to input another data file? 
(Ve Cae) S14 


end; %tend-while continue 


Rayleigh Wave Progression as a Function of Time and Distance Hodograms 
30 


30 
25 


20 


—_— 
qn 





-0.1 0 0.1 2 0.3 0.4 05 0.6 0.7 0.8 
time [sec] 
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Filename: sourcecorr.m 

Wertten by: F.E.) Gaghan 

Date Last Modified: 10 February 1998 

Purpose: This program is used to calculate the cross correlation 
between the input drive signal to the bass shakers and the 
received accelerometer and geophone signals. From the output 
of the cross correlation estimates for P-wave and R-wave 
speeds are labeled. 


Establish user interface to control data input 


diréctory=input(’ Input the data file directory: *,'’st); 
file directory=[directory]; 
file=input(*Input the data filename: 7s"); 


filename=[file]; 
eval(['load ' file directory,'\', filename] ); 


drive freq=input (‘Input the shaker drive frequency in Hz: '); 
rangel ft=input([’Input the range to the geophone #1 in ft: ']); 
rangel m=rangel ft/3.2808; 

range2 ft=-input( [> Input the range: vo ene geophone 12 in ft. 9 })-; 
range2 m=range2 ft/3.2808; 


freq Span=input (‘Input the frequency span in kHz: "); 
Sampling rate—t req span 1000* 2.56; 
Sauple delay—iInpul( Input the number wor delay canples a): 


Beeonudelength—oi 72, 
time scale=((sample delay/sampling rate): (1/sampling rate): 
((record length+sample delay-1)/sampling rate)]'; 


Assign input channels to correct data columns 
a@eel Long Jeel—il,;  aceel iver emeo!— 2, 
gece IM long col), gee ertgce l=, 
gece Vong 7eol—G;— Gee2 verre cel, 


Assign storage locations to data file columns 


accel Menger tmeSentes es 10 (27 eece eng gee, 
accel Vere tC imeser tes=xKrc 390: acce liver egeo 
geo Vong  etmeserires—xre3 90. geolelonggee! 


geol vert timeseries=xrc390(:,geol vert _ col); 
geo2 long timeseries=xrc390(:,geo2 long col); 
Gee Vere tameooen tes Krell) (1, geese vere col) 


s 


Generate the reference input drive signal 

N=sampling rate/drive freq; 

n=(0.20* Nii; 

drive signal=sin(2*pi*n/N) ; 

Ereme padding=Zeros(abs{sample delay), 1); 

padding-record length-{length(front padding) +length(drive signal) ); 
rear padding=zeros (padding, 1); 

reference signal=[front_padding;drive signal';rear padding] ; 


Create band pass filter centered at the drive frequency 
pass freq=drive freq/(sampling_rate/2) ; 

[b highpass,a_ highpass]=butter(2,pass_ freq, 'high'); 

[pe lowpess,; aulowpass | =buEter (2, pass freq); 
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Filter the accelerometer and geophone signals 
reference filtered= 

Filtfilt(b highpass, a highpass, reterencesoiguains, 
accel long filtered= 

filtfilt (> highpass; a highpass, accel longeeimescriacs 
accel vert filtered= 

filtfilt(b highpass,a highpass,accel vert timeseries); 
geol long filtered= 

Filet lee ianpass, ajhighpass, geoltlongurimeser ic... 
qecle verte erence — 

filtfivlet> nignpass,a highpass, geol vertytinessmn = = 
geo2 long filtered= 

fi lEEPIE (Db highpass,a highpass; geoZplong emeser ies): 
efaver’, Siccncis 16 aL, |Listshetstols 

filttvilc(b highpass,a highpass, geoZ (vert timeserices) 
reference filtered=filtfilt(b lowpass,a lowpass, reference signal), 
accel long filtered= 

Taidtiilt (se lowpass, 4a lowpass, acco ongeitilterca) ; 
eye eel ence ti71 betesacie i 

file tebeyoe lowpass, a loOwpassyaccch VemteG) Inpemec) 
geo long wirltered=filtiiie(s _lowpass, a _lowpass, Geol yong ste ee 
geol vert filtered=filtfilt (> lowpass,4a lowsass, geo) vere eilecmed 
geo2 _long_ filtered=filtfilt(b lowpass? a lowpass, geoZ _long_ filtered 
geo2 vert filtered=filtfilt(b_lowpass,a _ lowpass,geo2 vert filtered 


; 


; 


; 


; 


Calculate the cross correlation between the reference signal and the 

accelerometer signals 

[accel Plongecorn accel Longe lagc |= 
GORE (reterence Sadma ll, acco) Wlonyee Eber CG meee cilwaE, 

aecel long corr—abs (iiibere (accel 1 Ongpcoga ae, 

accel long lags=accel long lags./sanipling rate, 

Lqeeceiy Ver EReCorryacee ever rel aO ots 
RCOrr(reterence Signal, accel Vere srritencd, seoern |] 

accel Vert seOrr-abs (rbaere (aecel Verve mecorn) ) 

aceol vere lags—deeol vere lagse, Sali lincueduc, 


Plot the cross correlation as a function of time 

EULCiiise 1s) 

Sulop Poti); 

pilov (aceel longs lage radee mone ce oir) | 

title('Cross Correlation of Reference and Longitudinal Accelerometer 
Signals ye 

xlabel('lags [sec]'); 

Subp loc( Zia, 

plot (accel Vcmemlags,accel vertacorr |, 

title('Cross Correlation of Reference and Vertical Accelerometer 
Signals: 

xlabel('lags [sec]'); 


Bstablish the boundary criteria for the various wave speeds at 
geophone #1 

plower lag=rangel m/1000; pupper lag=rangel m/500; 

pwave lower=record length+floor(plower lag*sampling_ rate); 
pwave upper=record lengtht+floor(pupper lag*sampling rate); 
rlower lag=rangel_ m/100; rupper lag=rangel_ m/50; 

rwave lower=record length+floor(rlower lag*sampling rate); 
EWaVve.Wpper—recond *length+floor (ripper Vag sanerngeraee), 


164 


o\o 


Calculate the cross correlation between the reference signal and the 
longitudinal signal for geophone #1 
[fjeci eng teou.  geol tong lags|— 
eG (Sereronce signal, geoumlonguer tered, “cocit')), 
geo! long cerr—avs (hilbert(qceWbemgueeoerr |')/, 
[rmax long val,rmax_ long loc]= 
max (geol long corr(rwave lower:rwave upper) ); 
[long_val, long_loc]=max({(geol_ long_corr(pwave lower:pwave upper)); 
olepaue miaylivele wigs. 
1fplong val > gecl Jong cere, receqemuenagrh, 
pmax? Heng loc] long loc, 
else 
Cont ploe— 57 
Sjately | taJ9 Clercmdl cl gte, teu! 


reyes 


oe 


Calculate the cross correlation between the reference signal and the 
vertical signal for geophone #1 
[geol vert corr, geo! vert Jegs— 
RCOUE (PeTOreice Signal, gcol Vere fy eer Cc mecc uy 
geol vert corr=abs (hilbert (geol vert corr) ); 
[emax vert val, rmax vere lec|— 
max(geol vert corr(rwave lower:rwave upper) ); 
[vert val) vert lee)=max(geol vert scour (pwave fewer puavesUppe= i, 
Diievert, Valse Geele vere (eerr (recoucu lena any 
pMax vert locevert loc; 
else 
dont plot='y'; 
end; %end-if vert val 


oo 


$ Plot the longitudinal cross correlation as a function of slowness 
(inverse speed) 
slowness long=geol long _lags/(sampling rate*rangel m); 
slowness vert=geol vert lags/(sampling rate*rangel m); 
alos | eho) sue My ovlite 4s a 
Eee (2) 
Ssubplor (22.1), 
plee(slovness (lend geo! longmeenr a... 
slowness long(rwave_ lower+rmax long loc), 
geol, long cer: (rwave Vower+rmaxy tend grec, eo oe 
text (slowness, long(1)+0.0025, rmax long val—U.25, 
['rwave velocity = ', 
numzstr(17/sleowness Jong (twave. lower max lenguloe ee mig cee is, 
text (slowness long (1) +0 .00Z5, «smaxelong val=0 735; 
['rwave corr = ', 
numZstr(geol long _corr(rwave lowert+trmax long loc))]):; 
title('Cross Correlation of Reference and Geophone #1 Longitudinal 
Siomeuks:' )¢ 
xlabel('slowness [sec/m]‘); 
Subplot( 212); 
plot{slowness (vert, deGleverts Corr)... 
slowness vert(rwave_ lower+rmax vert loc), 
geo livieme .COre (rWave  lowertrmax vere, loa. G4, 
text(slowness vert {1)+0.0025,rmax vert val-0.25, 
['rwave velocity = ', 
num2str(l/slowness vert (rwave lower+trmax vert _loc)),' m/sec']); 
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Eext(Slowness vere!) 10.0025, rmax Vern cm gas Urecee 
('rwave corr = ', 
numZ2str(qeol vert core (rwave lower tena vememloc wale, 
title('Cross Correlation of Reference and Geophone #1 Vertical 
Siona is. )a 
xlabel('slowness [sec/m]'); 
else 
figuret Zi); 
SubpleetZ yer, 
plottslowrnescmong, geo! Jong corr, 22] 
slowness long(pwave_lowert+pmax_ long loc), 
Geclelcnagtcer: (pwave lower + pmax Slene os) row, 
slowness long(rwave_lower+rmax long loc), 
geol long corr (rwave_lower+trmax_long_loc),'o'); 
text (slowness long(1)+0.0025,rmax_long_ val-0.05, 
['pwave velocity = ', 
num2str(l1/slowness_ long(pwave lower+pmax long _loc)),' m/sec']); 
text (slowness long(1)+0.0025, rmax Ponce! tn, 
["pwave corr = ', 
numZ2str(geol long corr (pwave lower tomas mong aloc) |i); 
text (slowness long(1)+0.0025,rmax_long val-0.25, 
['rwave velocity = ', 
num2str(i/slowness long(rwave llowertrmax long tec))>' m/sec ae 
text (slowness long(1)+0.0025, rmax long _val-0.35, 
['rwave corr = ', 
numZstr(geol long corr(rwave lowertrmax long loc) )]); 
title('Cross Correlation of Reference and Geophone #1 Longitudinal 
Somes ) 7 
xlabel('slowness [sec/m]'); 
Supplet (202) ; 
Plot (slowmess Vert, geo! Vererceur 
slowness vert (pwave lower+pmax vert loc), 
geol vert _corr(pwave_lower+pmax vert_loc),'o', 
sicowness) vere (EWwave Vewer ile ioe 7 
Geol vere corr(rwave lower) (Maw vercEleoc), @ ); 
text (slowness vert(1)+0.0025, rmax vert val-0-035, 
['pwave velocity = ', 
num2str(1/slowness vert (pwave_lower+pmax vert _loc)),' m/sec']); 
texte (slowness Vere(!)7 050022) mnie ve rege! Oo 
[‘pwave corr = ', 
nuMzSEr(geel Vere Corre (pwaverlewertpma  longmice). |), 
text (slowness vert(1)+0.0025,rmax_ vert val-0.25, 
['rwave velocity = ', 
numZ2str(l1/slowness vert(rwave lower+rmax vert _loc)),' m/sec']); 
text (slowness. vemm() +0. 00257 emake eee ae ee, 
['rwave corr = ', 
numZ2str(qecl vert corr (rwave ewer? mace lencm oe) aly. 
title('Cross Correlation of Reference and Geophone #1 Vertical 
sploigigdlcny i 
xbabel (slowness [see/m|]*)> 
Siig Sene Fr "Aden. ple: 


Establish the boundary criteria for the various wave speeds at 
geophone #2 

plower lag=range2 m/1000; pupper lag=range2_m/500; 

pwave lower=record length+floor(plower lag*sampling rate); 


oP oO 
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pwave upper=record length+floor(pupper_ lag*sampling rate); 
rlower lag=range2 m/100; rupper lag=range2 m/50; 

rwave lower=record length+floor(rlower lag*sampling rate); 
rwave upper—-record length+floor(rupper lag*sampling rate); 


Calculate the cross correlation between the reference signal and the 
longitudinal signal for geophone #2 
rqeo2eTongueerr, geo2 long lags|— 
RCOMMLelCronece Ss. dnd W,dco2. leng@ee irene, -COCEE* ) ; 
geo? leongmeoerr—abs (hi lbert (geoZ ml oncmcecr a, 
[rmax  Hemguval, ymax slong eloe]/— 
Max (Gece wlong Corr (twave lover avavemnippel) 
[long val,long loc]=max(geo2 long _corr(pwave lower:pwave upper) ); 
done wel ou hee, 
If long val )> deo2 ond Gorn | recone mae) 
pmax long loc=long_loc; 
else 
dont plot='"y'; 
end; %end-if long val 


Calculate the cross correlation between the reference signal and the 
vertical signal for geophone #2 
[geo2 vert corr,geo2 vert lags]= 
xcorr (reference Signal,geo2 vert filtered, ‘coeff"'); 
Geo2 Vert Corr—aps (hi lbere(geo2e vere meemnm 
[Pd xavier e vol, (May avert aoe) = 
Max (G6OZ0Verte Gorm (mwave lower. evavesUpecm ai, 
[vert val,vert loc]=max(geoZ vert corr(pwave lower -pwave upper) ); 
Lf Vert. val > ge02 Vert corr (record lenge 
pmax verre noc —vert lec; 
else 
Pmaxeverce loc-imax vert loc, 
end wrect@eet. Vere vol 


Plot the longitudinal cross correlation as a function of slowness 
(Inverse speed) 
Slowness long=geo2 long lags/(sampling rate*range2 m); 
slowness vere=geo2 vert lags/( Samet migmrace canges m) > 
Pe ont plor—  ¥! 
1g Ak(e Deis wl EC N ies 
Subp loumi21 1); 
ploc(slowness long, geoz lonatecrn,..= 
slowness long(rwave lower+rmax_long_ loc), 
GeoZ long. Cor, (Yyweve over iam long elec) 7.16. )i: 
text (slowness long (1) 7070022) mex eng ia -0225, 
['rwave velocity = ', 
num2@str(l/slowness_long(rwave lower+rmax long loc)),* m/sec']); 
text (slowness Jong ()) 7070025, 2max long val-0. 35, 
['rwave corr = ', 
num2str(geo2_long_corr(rwave_lower+rmax long _loc))]); 
title('Cross Correlation of Reference and Geophone #2 Longitudinal 
Sigmals*) 4 
xlabel('slowness [sec/m]'); 
Subolor (212). 
plot (Slowness vert, g202  Verty Corr 4... 
Slowness vert (rwave_lowert+trmax vert loc), 
geo2Z vert corr(rwave lower+rmax vert loc),'o'); 
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text (slowness vere 1) 100025, rmax Veruny a0 ay 
['rwave velocity = ', 
numZ2str(il/slowness vert (rwave_ lowert+trmax vert _loc)),' m/sec']); 
text (slowness vert (1)+0. 0025, YMaxX Vere Veneaeeree, 
['rwave corr = ', 
num2str(geoz vert corr(rwave lower trmaxgiengutca) iia 
title('Cross Correlation of Reference and Geophone #2 Vertical 
Si otal au. 
xlabel('slowness [sec/m]'); 
else 
figure; 
SubplemeZ it ie, 
ploetstewnessn ong, geo2 Fong eecor™ 
slowness long(pwave_lowert+tpmax long loc), 
geee longveore (pwave lower 1 ena lenme mice) Oo yaa: 
slowness long(rwave_lower+rmax_long loc), 
geo2 Vong corr (rwave lower trmeaxs lengu lec). oF jy 
text (slowness long(1)+0. 0025, Imax “leng .Vval—02 05, 
["pwave velocity = ', 
num2str(1/slowness _iong(pwave_lowert+tpmax_long loc)),' m/sec: |e 
text (slowness long(1) +0. 0020, rmax “longer =O 75, 
[‘pwave corr = ', 
num2str(geo2 long corr (pwave_lowert+tpmax long loc))]); 
text (siowness long(1)10.0025, rma Ven vale =, 
['rwave velocity = ', 
num2str(i/slowness long(rwave_ lower+rmax kong loc) );" m/see wie, 
text (slowness long(1i)+0.0025, rmax Vong val—0335, 
['rwave corr = ', 
num2str(geo2 long corr(rwave lowert+trmax long _loc))]); 
title('Cross Correlation of Reference and Geophone #2 Longitudinal 
SLGmas 4); 
xlabelt “slowness [see/m]  }; 
subplot (212); 
plore (stewmess) Veur, Geol ay Cr uee oO iat ae 
slowness vert (owove lowertemaxgvert Loe), 


Geol Vert seer, (pwave Vewer ya) 7emiemel OG) o yen 
SlOWNESS Vere ( uWwave slower) miax | Verne roCc)., 
G€02Z VCrCe COrr  twaven lOwebrrmaxevcrue Oe) 7s © Ja 


text(slowness vert (])4+0. 0025; rmaxi vert val-0.05, 
['pwave velocity = ', 
num2str(l/slowness_ vert (pwave_ lower+pmax vert _loc)),' m/sec']); 
text (slowness vert (1)+0. 00Zo, amas Vere nvala 0, eS, 
['pwave corr = ', 
num2str(geo2 vert _corr(pwave lower+pmax long loc))]); 
text (slowness vere (1) +0]0025 ma vert va Ue, 
['rwave velocity = '', 
num2str(i/slowness vert (rwave lower+rmax vert loc)),' m/sec']); 
Lext(slowness vert (1)0, 0025 cere een le oe, 
[*rwave corr = ", 
num2str(geo2 vert _corr(rwave lower+rmax long loc))]); 
title('Cross Correlation of Reference and Geophone #2 Vertical 
Sa Onia ls) - 
xlabel('slowness [sec/m]'); 
ene send Pedont (plo: 
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Cross Correlation of Reference and Geophone #1 Longitudinal Signals 


rwave velocity = 88.7542 m/sec 
rwave corr = 0.89295 





O 
-0.03 -0.02 -0.01 0 0.01 0.02 0:03 
slowness [sec/m] 
Cross Correlation of Reference and Geophone #1 Vertical Signals 


rwave velocity = 83.2317 m/sec 
rwave corr = 0.94745 





0 
-0.03 -0.02 -0.01 0 Ogi 0.02 0.03 
slowness [sec/m] 


Cross Correlation of Reference and Geophone #2 Longitudinal Signals 


rwave velocity = 87.022 m/sec 
rwave corr = 0.84913 





0 
-0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 
slowness [sec/m] 


Cross Correlation of Reference and Geophone #2 Vertical Signals 


rwave velocity = 85.2163 m/sec 
rwave corr = 0.95461 





O = 
-0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 
slowness [sec/m] 
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Filename: sourcesignals.m 

Written by: F.E. Gaghan 

Date Last Modified: 10 February 1998 

Purpose: This program allows the user to display the timeseries of 
the filtered and unfiltered received accelerometer and 
geophone signals. The user can also display the hodograms 
and hankel plots of each of the signals. It also enables the 
user to isolate selected portions of the geophone signals 
For ease of display. 


Establish user interface to control data input 
directory=input ('Input the data file directory: ','s'); 
file directory=[directory]; 


file=inpub (input the data filename: "7 %s"); 
filename=[{file]; 
eval({'load ' file directory, '\', filename] ); 


drive freq=input({'Input the shaker drivegercquctiGugin Hz: ~ } 

range ft=input({'Input the range to the geophone #1 in ft: ' 

range m=range ft/3.2808; 

Eregueerpai—imnpue( Input the frequency span ema kha: =); 

sampling rate=freq span*1000*2.56; 

sample delay-input (’Input the humber ot delay samples: ~},; 

record length=3192; 

time scale=[(sample delay/sampling rate): (1/sampling_ rate): 
((record lengthtsample delay-1)/sampling rate) ]'; 


Ds 


Assign input channels to correct data columns 
acce Rl ongmeot. | .aeeel  Verecol—Z, 
geol Mong teel—s; Sgeco verre neol= |; 
geo2 long col=6; geo2 vert col=3; 


Assign storage locations to data file columns 
accel Vong wet@mesertes sre] 30 acee Te ongmec ls, 
accel vert eameseries-xre 370 (>, deceiver coing, 
geol long (elmeser tes nee l [0 | gece longeo ey, 
geol verte Eameser tres-—xre 730 geOln vere ecols), 
geo2 long timeseries=xres90(:, geez long col); 
geo2 vert timeseries=xrc390(:,geo2 vert col); 


Plot unfiltered timeseries of accelerometer Signals 
frgure i): 
Subplor. 2), 
plot(time seate, accel lengreimeseries): 
axis ((Eimenscale (ij timepeea l= (eceo- cm lengen 
min(accel long timeseries)-0.2,max(accel long timeseries)+0.2]); 


title({'Longitudinal Time Series of Accelerometer Signal for File: ', 
Fert): > 
xlabel('time [{sec]}'); ylabel('amplitude [V]'); 


Supp ler. 21k 
plot (time scale, accel vert timeseries); 
axis({time scale(1l),time scale(record length), 
mainiaece] long timesertes) —0.2,max(accele tong timeserics ae 17 


title({'Vertical Time Series of Accelerometer Signal for File: ', 
Eile); 
xlabel( timely (sec) );vylabel (vanplitrde Vivi je, 
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Create band pass filter centered at the drive frequency 
pass freq=drive freq/(sampling rate/2); 

femivignipace scm tghipass|=butter(2, pass freq, “high’); 

[pe lowpass, a) Owpass ||—butter(27 pass tired) ; 


Filter the accelerometer and geophone signals 
aG@e | PEOnG Man eerea— 

Piet tet ben gihpass, a highpass, accel long timeseries); 
gece ly eee tite! ec 

fr Tet eee hgipass, ato nbass7ceec myo me LIMeser Tes). 
eteiolls Ieiele | seal ee ie eyo= 

frlttile(s tnghpess,a highpass, gesclmongeeimesertss); 
Pool Vert ee eered= 

filet be pep echeacs, a highpass dqeOlsveu timeserles),; 
Seo? onget eit ered— 

PLC t Peon ore tpacs, ae iighpdss,QcO-n bom metnicscer logy, 
geo2 vert filtered= 

ERlChi te (So hwgapass, a ieipass, geo? Veuemeimeceries) ; 
eee long ey fhe rea 

PileLile (eelowpass,a lowpass , accel ploncge el ber ea); 
accel vert filtered= 

EPEGTRE (a heOwpass, an OWpass, aceon vert yt eleemead) | 
geole rong -erkterea= 

valle stale (Jel lev ieyelerer el Weverctoc este ll llere; sen die cic, | 
geol vert filtered= 

filttilt(b> lowpass,a lowpass, geol vere filtered), 
jeoZeteng til EeCrea=— 

PibGr be QovPacc jon LOwodss, GeO4 sheng eg ia bechoe js, 
Geo2 vert tiltered= 

fi lerile( > lowpass, a lowpass, geo" Vert filtered); 


Plot the filtered and unfiltered timeseries signals for geophone #1 
figure (2); 

Sue lori 2 iy) - 

PLO timescale ceo Pe rongsEImeser les), 


title(['Longitudinal Time Series of Geophone #1 Signal for File: ', 
fake | 
xlabel('time [sec]'); ylabel('amplitude [V]'); 


Supp home) 

PISCE (eime (Seale, geo! Mong pr beered ); 

title(['Vertical Time Series of Geophone #1 Signal for File: ', 
file] y; 

xlabel('time [sec]'); ylabel('amplitude [V]'); 


Plot the filtered and unfiltered timeseries signals for geophone #2 

Pigunets |: 

Subploet2 lr); 

Pigs Lime seats, deo7 Vong -imlescrics), 

tittie([*Longrtudinal Time Series of Geophone #Z Signallter Fale: 
fied je 

xlabel('time [sec]'); ylabel('amplitude [V]'),; 

SUS plore 12) ] 

plot (time scale, geoZ long filtered); 

title(['Vertical Time Series of Geophone #2 Signal for File: ', 
eolkd Woe a Be 

xlabel('time [sec]'); ylabel('amplitude [V]'); 
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Plot the hodogram of the filtered and unfiltered timeseries signals 
of the accelerometer 

figure (4); 

subplot (211); 

plot (accel long timeseries,accel vert timeserieés); 
xlabel('longitudinal amplitude [V]') 
ylabel('vertical amplitude [V]'); 

axis equal; 

Suppletlziz yg 

plce (acce! Migngyemieered accel Vert nit ieewead):, 
xlabel('longitudinal amplitude [V]'); 
ylabel('vertical amplitude [V]') 

axils equal; 


Plot the hodogram of the filtered and unfiltered timeseries signals 
for geophone #1 

elie purara tc ar 

subplot (211); 
picel(cgeolelengetimeseries, deol vere elles coi: 
Riedel “Tlengieudinal amplitude ty iy; 
ylabel('vertical amplitude [V]'); 

axis equal; 

SUbproE (212) 

plct(geol, long fittered, geo leverie (1 tee wed); 
xlabel('longitudinal amplitude [V]'); 
ylabel('vertical amplitude [V]'); 

axis equal; 


Plot the hodogram of the filtered and unfiltered timeseries signals 
for geophone #2 

figure(6); 

Suppor (2c 

DICE GeCOZeleng Limeser ies gee7 Vet ume ies ctres,) |, 
xlabel('longitudinal amplitude [V]') 
ylabel('vertical amplitude [V]'); 

axis equal; 

SuCpleti(Zi2) 

plet (geez long Evitered, dec2gverestyrercad), 
xlabel (’ longitudinal ampli Eeuder (i ino. 
ylabel('vertical amplitude [V]') 

axis equal; 


Display hankel plot for geophone #1 
YZ axlS—zeros( Pengen (timescale) 


y_maxl=max (abs (min(geol long timeseries)),max(geol long timeseries) ); 
y axisl=linspace(-y maxi, yomaxl engin (cameuseale) i: 
z maxl=max(abs(min(geol vert timeseries)), max(geol vert timeseries) }> 


2 aos.) iispace (- z maxi, z maxl, Pengeh (time isco bei 
y_ maxlf= max (abs (min(geol long filtered)),max(geol long filtered) ); 
vy -axteli=linspace(—-y maxli, yo maxlt, lenge eamersea re a7, 
Zz maxlf=max (abs (min(geol vert _filtered)),max(geol vert filtered)); 
Z axis ilr—linspace(—z max lt, 2 meme lenge lemme ssaate sy, 
figure(7); 
SubpROE 2b) 
pPlole(time wcale,cgeol tong ttimesertes, Geet iver ret aneseries a 
TAME CSCalLe) xX V2 Axis, ay Zeanlo, 
BYZ axis; y axisl  eyZcaxts, ka, 


2 


MY Zlaxbo, XYZ axis,Z axisl,*k’ ); 
xlabel('time [sec]'); 
ylabel('longitudinal amplitude [{V]') 
zlabel('vertical amplitude [V]') 
view(15, 30); 
Sulgelot (212); 
plot3 (time scale,geol long filtered,geol vert filtered, 
time SGale,xy2 axis;xyz axis, ki; 
XYZ Gxkisye sexs PE xyclaxds § le, 
KYZ GRUS ey ceaxXlS,.2 Axisti 7 k ji; 
axis ([time Scake (212) time wecale( Vergem timerscale)), 
min(geol long filtered), max (geo Werong yeiimered), 
min (gqeol vere filtered) pmax (geolmvereyrrirered) |); 
xlabel('time [sec]'); 
ylabel('longitudinal amplitude [V]') 
zlabel('vertical amplitude [V]'); 
view(15,30); 


Display hankel plot for geophone #2 
Vomak2 tea a aes (Mem dee7 Ong uetmescmtes) i Maw GoO7 beng eeEmeser ies )), 
Yak TS2—1aneepace | > an2, eax, POngE ll elma sea vee): 
Z. _max2=max (abs (min (geo2_ VeGEREIMeSet tec) |) Maw GcOZmVCbE mL Imcoct Les) jG 
z axis2=linspace(-z maxZ,z _max2, length(time scale) ); 
y _max2 f=max (abs (min (geo2 Bitona PEC ered )))) Max ceo7 longs filtered )); 
Vy axkis2t—linspace (-y maxZt,y maxZ2t, lengeh (time seale))7 
Z MaxZ2i=maxtabs(min(geoZ vert filtered) ),max(gqeo2 vere filtered) ); 
z axis2f=linspace(-z max2f,z maxZ2f,length(time scale)); 
fagure (b) > 
subplot (Zi iis 
plot3(time scale, geo2 long timeseries,geo2 vert timeseries,... 
time scale,xyz axis,xyz axis,'k', 
XYZ axisS,y axis2,xyz_axis,'k', 
XYZ aKES XYZ axis, 2 axis7, i); 
xlabel('time [sec] ') 
ylabel (' loneueuaine amplitude [V]') 
zlabel('vertical amplitude [V]') 
View, 15230). 
siulejede nen GZ ae 
OES ears emma ace eee on ie 
time Scale,xyz axis,xXyz axis,'k', 
RY 25akiS, VY dxisZe eye axis, 
RV A2Gax1S, My oeakis, Zon sZt, ok - 
axis([time scale(512),time scale(length(time scale)), 
min(geoZ ong filtered);,max(geoZ Tong filtered), 
min(qeoZ vert filtered) ,max(geoZ vert filtered) ]j; 
Xxlabel('time [sec]'); 
ylabel('longitudinal amplitude [V]') 
Zlabel('vertical amplitude [V]') 
view(15, 30); 


oo 


ore 


Enable the user the ability to display selected portions of the time 


record 

Gas = a ye: 

isolate=input('Do you want to isolate the Rayleigh wave? 
(iO “ee roe 

if isolate=='y' 


continue— yy 


while continue~='gq' 
GaSe ae, 
tanec le i= 
input( "Input the start time of the hodogram foreirgune(y ee 
point! l=floor(timel 1*sampling rate); 
time2 1= 


input ('Input the stop time of the hodogram for figure(7): '); 


point2 I=floor(time2 1*sampling rave); 


means! —ones (point? 1-poineigis tli eome sea te (oon > ar 
Viceoeeimescr1es— 


Zz 


y 


zZ 


uy 


Za 


y 


Za 


o 


= 


{o) 


max(abs(min(geol long timeseries(pointl lspoint2 1))), 

max(geol long timeseries(pointl] l:point2 1))); 
Dia me Ser res= 

max (abs (min(geol vert timeseries(pointl l:point2 1))), 
max (geol vert timeseriecs pornr Wr pete Zl) ye 
Sans Wet imeserics— 


linspace(-y maxl timeseries,y maxl timeseries, 
poine2 l-peint soa, 

Seve reabuleerere les 
HinSpace (2 Maks ties. sles eer ed Moccia oar 
Point Ze) point) lie 


“Maxi tiltered= 
max (abs (min(geol longer tered poles een Zm. ae 
max(geol long filtered(pointl l:point2 1))); 
Mier es eiena scl 
mMax(abs(min(geol vert) Filtered (pointle pointZ a aay 
Man Geollvert tiibered (porn r We Cones a) r- 
Pep ciedl Maik weic sol 
linspace (-y_maxl filtered,y_maxl filtered,point2 l-pointl 1+1); 
waxveonee ease acel— 
linspace(-z maxl filtered, z maxi filtered; peint 2) spoine tile 


Plot hodogram of user selected portions of geophone #1 
faqune cs je 
plotigeol Vong trtlvered) polite eine ay 

geol verteitl tered (poli: Ir Cems ee 
xlabel('longitudinal amplitude [V}'); 
ylabel('vertical amplitude [Vj'); 
axis equal; 
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$ Display hankel plot of user selected portions of geophone #1 
saben agl a (ac 
Supe vor tzw.) 
ple oieimerscale( point ly pocmmes (a), 
Gee Cage Nese es (gerne folie |), 
Cees Ver emriIMese: les pommel pollu? 1 )7 
time scale(pointl 1),geol long ee eg Oey 
Gee vere _timeseries(point1_ (le ie aaa © a ene ee 
time wseale (oOIntZ) die geol_long timeseries(point2 1), 
Geol vere ne meson tes. (OO mie) eee 
Pime —seele (point ll) sore. bi. 
RY2 aAxUstee int Wl spomnrZ gay 
Syenexrreqeormel DT speimez Ws) ste a: 
Retro yeare sl ElMesemres, 
my Zea eou Orme. ls poune2 slices, 
Moraceae aS (point lel pomme2 wa) | 
Ze he beetles er Wes 7 kiln 
xlabel('time [sec]'); 
Vlabel (’ long@rudanal amplitudes! 7]. \; 
Zlabpel('verticaleamplitude [Vii i; 
axis (| Peime scale(poimel ))) eimeyscalies ponies ly), 
min(geol long timeseries),max(geol long timeseries), 
Min Geol Vere eimeser ies | ten (Goomene memes omits). \e 
VLewl 1S, 30) : 
Sub pTOE C212); 
DICE ls Veime scale (poimmt les peineas a, 
geole ong ee tered (pote he: pene 
efeSlk Neuere pectiliccsaslelversehe lL Jlkgjecuuigie 21 iL) 
time scale(pointl 1), gael lone Seeaeciisonmel_ i, 
qeol vere tf Pitered (poi ce eae wcue.a a 
time scale(point2 1),geol ae os Se 
GeOUNVernee lt rered (polnecuay ee 
time _scale(pointl_ Terie Zamey., 
xyZ axes (point) 1: point2 1), 
RV 2 rd aes Oa le Od me 2 § jy kya 
Roakitsi, yo ax Slee Geued, 
ra ebcalts\(Qofenl ghey a sjoye Weue f. Mljie 106 7 
Spake 1 xy 2 vexis (oCIme ly spermine ly, 
Cameos ae Ge ee hs de, 
xlabel( time [seel’)> 
ylabel('longitudinal amplitude [V]'); 
Zlabel('vertical amplitude [V]'); 
axis ([Cime Seate (pointe), winerscale (poms an, 
min(geol long filtered),max(geol long filtered), 
Min | Geol ever yeti tered) max (gee Vomege rl tcmcena.) 
View (lS. 3004 


das psa, 
timel 2= 

input Ilmpwt the stare timeror Ene Nodeograr toner loure( 3): i, 
Peici 2st loor (cimel 2 sanelimgs wae. 
time2 2= 

input( "Input thersstopo time of the hododran sor i1guveta)..  j 
poinczs42—-tloor(times 2" sampling eace 
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x 


\ a 


fi 


ies 


2 


se 


z 


ae 


Zz 


_axis2=ones (point2 2-pointl 241, Y= timemeeal— (penn. lezen 


max2 timeseries= 
max(abs(min(geo2 long timeseries(poime ly pointe ee 
max(qeo02 Tong timeserics (pointe ly 2 pela ae. 

_Max2_ timeseries= 

~max (abs (min(geo2_ Verte tLimeseries (botnet ene) am, 
max(geo2 vert _timeseries(pointl | Ze point2_ Ze: 

axis2 fimeseries= 

linspace (-y max2 timeseries,y max2 timeseries, 

POtne Aero ne Woz); 

axis2 timeseries= 

linspace(-z max2 timeseries, z mMax2 timeseries, 

Pome? ererne hoe); 


max2 _filtered= 
max (abs (min (geo2 long  falrered)(Pommelee  pOltmezm sa), 
max(geo2 long filtered(pointl 2: point2 2))); 


Ss ical Ih nienetete ee 


Maxtabs(mim(gqec2e vert filtercaipemmee omnes) oe 
tax (GeO? vert filtered point ley aee mies ee eo. 
a<us2 tii tercad— 

Mmspace(=y Max2 ti ltered? yy maxi oreqr 
POrne2. 4 polnel 2, 
_axis2 filtered= 

linspace(- ZaMaxcoutltered. 2 Maw t et Eonca, 

Pomme 22 -pe melas 


Plot hodogram of user selected portions of geophone #2 
figure (6); 
pPleoe(gqeo? Wong [1 iremed (polme lA poine Aueae 

GeO2 VErent i bencd (pore 2 pone ae, 
xlabel('longitudinal amplitude [V]'); 
ylabel('vertical amplitude [V]'); 
axis equal; 


Display hankel plot of user selected portions of geophone #2 
fLigune( LO) 
Supplor (2 1) 
pICtS(Timersea le (povnb 2 oO lice a, 
Seo7 Vongeeimeserles (Petar lee ieee 
deo2 Vere eImesertes ( pomne a2 perme a es 
time scale (point iy2) -ocozelencm timeseries (point1_ ae 


(eOZ Vere EIMeser ies (Seren rma) ero ue eee: 
time scale(point2 2),geo2 long So TMeSe Pe ecet ees Zoe 
GeO2e 7 er eet Miles mies (Ol ee ere 


Lime Scale ( se tari seo Ehren) a 
KY Zak Ps(pOlne Pe eOimes a. 
xyz axis(pointl 2:point2 2),' aa eee 
x axis2,y axis2 timeseries, 
KVZ axis (POlMei 2 polmezazye ik 
AXIS; X¥C axtepolnele. pomluzes a 
Z axis2 timeseries,'k"'); 
xlabel('time [sec]'); 
ylabel('longitudinal amplitude [V]') 
Zlabel('vertical amplitude [V]') 
axis([time scale(pointl 2),time scale(point2 2), 
min{geo2 long timeseries),max(geo2 long timeseries), 
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oe 


oo 


Mmim(geo2 vert timesseries | ;max(geo2.vert tiameseries) |); 
view 57 30); 
SUDD Lot U2). 
ploOestt imensecale (point yo 3permeZ 2), 
djeoz nr ong wav ltereca pommrein2  poant2 2), 
Jeo7 *Verenr bberod (pOmmmleZ spOlIne2 2), 
time scale(pointl 2), Belong fered per ntieey, 
GeO7 Verne em ttered Oosme VEZ iO ee. 
time scale(point2 2),¢ge02_ eke egg ee Ze, 
Geo2 Vere ii ened (pole 7 
time: scale (pointin2  pornes 2), 
RY Zax S Ome le OO] liam), 
RY 2oaxis (POInNe l= Ol nee ney ee 
pe Care omy at cva ee 1 Jie (iclere |) 
XYZ axtS (pOrmEel 2 POineZe 2) ks 
Mase seks (POlme lye point2_ 2), 
Z Pic Pim tered, kh; 
xlabel('time [sec]'); 
Virabeite longitudinal amplitudes Vi): 
Zlabel('vertical amplitude [V]'); 
axis({|[time scale(pointl 2),time scale({point2 2), 
iat (Jeez m long pri weercad) Max (Geog Olge biecrea) , 
min(geo2 vert filtered),max(geo2 vert filtered) ]); 
VMS We Me 0, 


Graphically display ratio of longitudinal to vertical signal 
amplitudes 
dispt’ “) 
ellipse= 
input ('Do you want to calculate ratio of elliptical axes 
(y Ofte 7 S }; 
if ellipse=='y' 
Jeol one sniilbert= 
deeper e (Geol, long i: ltered (pein: be = peime2 1)')}> 
geolgvert jailbert= 
abs (hilbert (geol_ VerGlr reread (point) Iypeint2el a 
geclwemp matito-geo! long snilbert. /geol_vert hilbert; 
£ Powee. (1 hie: 
subpRoEt 21 1); 
pleoUerncwseale (point l 1-256. peineZe a 256), 
geol@liong firltered( point) 1-256: point2 14256), 
i Memeee Le (POINEl Tl —>256: pOlmuZ #256), 
Ceo ewe weet) | tered (point bel] 2eG- pol me2 0) ame 
xlabel('time [sec]'); ylabel('amplitude [V]"'); 
beqencir wongrtudineal 5 verticel* ); 
SubpleGizi2) ; 
plOu Game toca le( point ll point? 1), Geo lwamee ra ete, 
Xbebeis@erame (sec) *); ylabel oi a ie 
axis (ibmme s scale (pointl yl) time scale (pe une? a) He ie ers 


Lae 


geo2 Vong ihaitent— 
abs (halbert(geo2 long filtered(pointl 2 vomit Zee 
geo2 vert hilbert= 
abs (hilbert (geo2_ vert filtered(pommule 2 peinez 2 ji; 
geo2_ amp ratio= geo2_ long_ hilbert./geo2 vert hilbert; 
Pigume( Zi", 
Subp Lome ii, 
PlOE(timerscale(point] 2-256: pointzeZ aco), 
G@eo2 Tong filtered (porn li 2225¢ polmec nao), 
time scale(pointi, 2-256: poimezes.2oo, 
geo2 vert filtered point le2-2aG. po rec 242 Oi sor 
xlabel('time [sec]'); ylabel('amplitude [V]'); 
fegena Wonga tudinal" ,“vertiealeys 
SUbpplOul 22); 
plet(time scale (pointl 2 perme a2 Wide o2 sano ecto, 7 
xlabel('time [sec]'); ylabel('longitudinal/vertical'); 
axis ( [time scale(porntl 2), Gime seater eoinrs 92770 i 


longitudinal/vertical 


am plitude [V] 
So 


end; %end-if ellipse 


WE oa ae any key to continue isolacuon Or "G4 te Gull. ey 
' S i] shes : 
end; %end-while 
end; tend-if isolate 
x 10° Time Series of Geophone Signal for File. sh7950p4 x 19 Hodogram 


— 





-1 __ longitudinal 
--- vertical 

-2 

02 025° 025 035 04 045 05 055° 06 ~ 065 
time [sec] 

Ratio of Particle Displacements for File: sh7950p4 

0:5 

0.4 

U3 

0.2 

0.35 0.4 0.45 0.5 0.55 
time [sec] 
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